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Abstract 

Two types of Gamma-ray bursts (GRBs) are observed: short duration hard spectrum GRBs and long duration soft 
spectrum GRBs. For many years long GRBs were the focus of intense research while the lack of observational data 
limited the study of short-hard GRBs (SHBs). In 2005 a breakthrough occurred following the first detections of 
SHB afterglows, longer wavelength emission that follows the burst of gamma-rays. Similarly to long GRBs, afterglow 
detections led to the identification of SHB host galaxies and measurement of their redshifts. These observations 
established that SHBs are cosmological relativistic sources that, unlike long GRBs, do not originate from the collapse 
of massive stars, and therefore constitute a distinct physical phenomenon. One viable model for SHB origin is the 
coalescence of compact binary systems (double neutron stars or a neutron star and a black hole), in which case SHBs 
are the electromagnetic counterparts of strong gravitational-wave sources. The theoretical and observational study 
of SHBs following the recent pivotal discoveries is reviewed, along with new theoretical results that are presented 
here for the first time. 
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1. Introduction 

Gamma-ray bursts (GRBs) are short intense flashes of soft (~MeV) 7-rays that are detectable once or 
twice a day. The observed bursts arrive from apparently random directions in the sky and they last between 
tens of milliseconds and thousands of seconds. Their physical origin remains a focal point of research and 
debate ever since they were first detected, more than three decades ago. 

More than twenty years ago hints that the GRB duration distribution is bimodal have emerged (Mazets 
et al. , 1981; Norris et al. , 1984) , suggesting that the GRB population might not be monolithic. The idea 
that GRBs are most likely composed out of two major distinctive sub-populations became the common view 
following the influential work of Kouveliotou et al. (1993). Based on 222 GRBs detected by the Burst And 
Transient Source Experiment {BATSE^ ) on board the Compton Gamma- Ray Observatory, Kouveliotou et 
al. (1993) confirmed that the GRB duration distribution is indeed bimodal, with a minimum around 2 s. 
They have further shown that bursts with durations shorter than 2 s are composed, on average, of higher 
energy (harder) photons than longer bursts. It took a dozen additional years before recent observations 
confirmed that the two sub-populations, defined in duration-hardness space, indeed represent two distinctive 
physical phenomena. Throughout this review I use the term short-hard GRBs (SHBs), or simply short GRBs, 
for bursts that share the physical properties of the population that constitutes the majority of the short- 
hard events in the duration- hardness space ^ . Bursts with physical properties common to the population 
that constitutes the majority of long-soft events are denoted simply as long GRBs. The term GRBs is used 
to describe the whole population of the busters, short and long. 

A major revolution in GRB research came with the launch of BATSE in 1991. BATSE detected bursts 
at a rate of one per day, a quarter of which were short. The isotropic distribution of bursts on the sky (both 
long and short) and the paucity of weak bursts indicated that the origin of GRBs is extra-Galactic and most 
likely cosmological. However, because of the inability to identify the distance of any specific burst based on 
the observed gamma-rays alone, the final confirmation of the cosmological origin of GRBs had to wait (for 
a review of BATSE results see Fishman & Meegan, 1995). The final validation that long GRBs arrive from 
remote locations in the Universe came following the discovery by the Dutch-Italian satellite BeppoSAX that 
the prompt 7-ray emission of long GRBs, lasting minutes or less, is followed by X-ray emission that can 
be detected for hours and days (Costa et al. , 1997). Observations of this X-ray emission enabled accurate 
localizations of bursts and led to the detection of associated optical (van Paradijs et al. , 1997) and radio 
(Frail et al. , 1997) emission, that can be observed, in some cases, for weeks and years, respectively. Longer 
wave-length emission that follows the prompt GRB is called the "afterglow" and it provides a wealth of 
information about the physics of the burst. Most importantly, sub-arcsecond localization of optical afterglows 
enabled secure identification of galaxies that host long GRBs, leading to the measurement of their redshifts, 
and finally establishing an unambiguous distance scale for these events (e.g., Kulkarni et al. , 1998). In some 
cases, the burst redshift was measured directly by detection of absorption lines in the spectra of optical 
afterglows (e.g., Metzger et al. , 1997). 

These observations revolutionized the study of long GRBs. The confirmation of their cosmological origin 
and the detailed observations of a growing sample of afterglows supported a physical picture according 
to which long GRBs arc produced by a catastrophic event involving a stellar-mass object or system that 
releases a vast amount of energy (> O.OOIMqc^) in a compact region (< 100 km) on time scales of seconds 
to minutes. This energy source, referred to as the "central engine" , accelerates an ultra-relativistic outflow 
to Lorentz factor > 100 and it is this outflow that generates the observed 7-ray prompt emission, and later 
the afterglow, at large distances from the source. This model is now generally accepted. 

The nature of the stellar progenitor of long GRBs remained a matter of debate for many years. Until 2003, 
when a detection of a long GRB at a relatively low redshift (z=0.168) by HETE-2^ led to the identification 



^ http:/ /www. batse.msfc.nasa.gov/batse/ 

^ Note that it is not necessary that the duration of a specific burst corresponds one-to-one to its physical properties (e.g., the 
nature of the progenitor). It is possible for example that there is a small fraction of bursts that are associated physically with 
long GRBs but their duration is short and vice versa. 

^ The High Energy Transient Explorer; http://heasarc.gsfc.nasa.gov/docs/hete2/hcte2.html 
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of a type Ic supernova (SN) spectrum superposed on the afterglow of this burst (Stanek et al. . 2003; Hjorth 
et al. , 2003). This observation confirmed previous suggestions that at least some long GRBs are associated 
with SNe (e.g. Galama et al. , 1998; Bloom et al. , 1999). Following this association and several additional 
evidence (e.g., Fruchter et al. , 2006) the consensus today is that most, and probably all, long GRBs are 
produced by the collapse of very massive stars (e.g., Woosley, 1993; Paczynski, 1998; MacFadyen & Woosley, 
1999). 

Comparable studies of SHBs were not conducted during this time. Being shorter and harder they eluded 
accurate localization and no SHB afterglow was detected despite the effort. The breakthrough in the study 
of SHBs occurred, finally, during the spring-summer of 2005 when Swift and HETE-2 succeeded in localizing 
several SHBs, leading to afterglow detections and to the determination of their redshifts. Following these 
discoveries, and fueled by a continuous flow of SHB detections by Swift, the study of SHBs progressed 
rapidly. The flrst conclusions from these observations were that SHBs are cosmological, but unlike long 
GRBs, their progenitors are not massive stars, thereby confirming that these two observationally defined 
classes are distinctive physical phenomena. On the other hand, the comparable luminosities and roughly 
similar afterglows of long and short GRBs, suggested that similar physical processes are involved in both 
types of explosions. 

The aim of this paper is to review and summarize the theoretical study of SHBs following the recent 
observations. These observations are described in §2 and aspects that have direct implications for the current 
theoretical understanding of SHBs are emphasized. Special attention is given to the observational challenge 
of identifying a burst as a SHB. 

Relativistic effects and the prompt emission theory model are discussed in §3 while the afterglow theory 
is discussed in §4. The theory of both the prompt and the afterglow emission is discussed in the framework 
of models in which the emission is produced by a relativistic outflow that flrst expands quasi-spherically 
(see definition in §3.1.2) and that later its interaction with the ambient medium leads to the formation of a 
blast wave that propagates into this medium. These models include the "classic" fireball model of radiation 
accelerated baryonic plasma as well as other types of relativistic outflows (e.g., Foynting-flux-dominated) * . 
The physics of GRB emission (prompt and afterglow) was intensively studied in the context of long GRBs, 
while there is little work speciflc to SHBs. As both the prompt emission and the afterglows of long and short 
GRBs are quite similar, long GRB models are typically applied to SHBs as well. The fireball and related 
models were comprehensively covered by several excellent reviews, and I refer the reader to these reviews for 
a detailed description of these models (Piran, 1999, 2000; Zhang & Meszaros, 2004; Piran, 2005b; Meszaros, 
2006; Lyutikov, 2006). In order for this review to be self contained, I include a brief overview of the basics of 
the fireball and related models, emphasizing aspects that are relevant to SHBs. I focus on the interpretation 
of SHB observations in the framework of these models. In several cases I extend existing models and re-derive 
predictions and constraints that are relevant to SHB parameter space. These derivations and some novel 
conclusions that follow appear for the first time in this review (see highlighted sections in Table 1). 

I go on to review proposed SHB progenitor models and the formation of the central engine §5 (for an 
excelent comprehansivc review of the topic see Lee & Ramirez-Ruiz , 2007). The leading progenitor candidate 
is the coalescence of a neutron star (NS) with another neutron star or with a black hole (BH). This idea 
was commented upon some twenty years ago (Blinnikov et al. , 1984; Paczynski, 1986; Goodman, 1986; 
Goodman, Dar & Nussinov, 1987) and was first explored in detail by Eichler et al. (1989). It successfully 
survived twenty years of observations and got some support from the properties of SHB afterglows and of 
the identified host galaxies. I therefore dedicate most of the discussion on SHB progenitors to the aspects 
of these mergers that are relevant to SHBs. Nevertheless, a merger origin for SHBs is not confirmed, and 
other progenitor models are still viable, and are discussed here as well. 

Mergers of NS-NS or NS-BH binaries are also the most promising sources of gravitational-waves that 
may be detected by ground-based observatories. This of course makes SHBs very interesting, as the possible 
electromagnetic counterparts of gravitational-wave signals that are expected to be detected within a decade. 
The prospects of gravitational wave detection from SHBs are discussed in §6. 



^ Other GRB models, such as the cannonball model (for review see Dar & de Riijula, 2004) and the precessing Jet model (e.g., 
Fargion, 1999; Fargion & Grossi, 2006), arc not discussed in detail in this review. 



4 



For convenience, Table 1 summarizes and compares observational properties of short and long GRBs, 
and their theoretical interpretations in the fireball model framework. This table can also serve as a quick 
reference guide for the different topics covered in this review. Through the review flm = 0.3, SIa = 0.7 and 
Ho = 70 kms^^ Mpc~^ cosmology is used. 

2. Observations 

Gamma- Ray burst observations are grossly divided into two main phases - the prompt gamma-ray emission 
and the afterglow. The prompt gamma-rays are intense bursts of ~ MeV 7-rays that are detected by gamma- 
ray space observatories and are localized by current instruments to within several arcminutcs. Afterglows 
are X-ray, optical and radio emission that follow the prompt gamma-rays and can be observed in some cases 
weeks, months and years after the bursts, respectively. 

BATSE detected the prompt emission from ^ 3000 GRBs, about 1/4 of which are SHBs. Additionally, 
prompt emission from more than a hundred SHBs were observed by Konus-Wind^ and other spacecraft 
that were, or still are, part of the Interplanetary Network^ (IPN). This is compared to a handful of SHB 
afterglows that were observed following accurate localization of the prompt emission by Swift and HETE-2. 
In some of these cases the sub-arcsecond localization of the afterglow led to an unambiguous identification 
of the burst redshift and host galaxy. These observations are briefly reviewed below. 

2.1. Prompt emission 

Most of the prompt emission properties are derived using BATSE bursts. Typically, the SHB sample is 
drawn out of the complete BATSE GRB sample using the criterion Tgo < 2 s (defined below). This is clearly 
a rough cut (see §2.7), but it is sufficient for the purpose of statistical studies which are weakly affected by 
a small contamination of the sample. 

2.1.1. Duration 

The true (intrinsic) duration distribution of SHBs is unknown. At the short end, the observed distribution 
is affected by the minimal BATSE trigger time (64 ms) while the long end the distribution blends into that 
of long GRBs. Figure 1 depicts the Tgg distribution of the entire BATSE GRB sample, where Tgg is the 
duration encompassing the 5'th to the 95 'th percentiles of the total counts in the energy range 20 — 2000 
keV. The duration distribution is bimodal with a minimum around 2 s (Kouveliotou et al. , 1993). For this 
reason the dividing line between short and long GRBs is usually drawn at Tgg = 2 s. This distinction suffices 
for the purpose of statistical population analysis, but one should bear in mind that there are short GRBs 
with Tgo > 2 s, and long GRBs with Tgo < 2 s. Horvath (2002) finds that this bimodal distribution can 
be decomposed into two lognormal distributions, as presented in figure 1. The SHB distribution peaks at 
Tgo ~ 0.8 s, and the full- width-half-maximum of the distribution is 1.4 dex. If this decomposition represents 
the physical one then 25%[3.5%] of BATSE SHBs last longer than 2[10] s. 

It was recently realized that at least in some SHBs the initial short and hard 7-ray prompt emission is 
followed by a much longer X-ray "tail" that lasts tens to hundreds of seconds. Typically, this tail is too faint 
to affect BATSE Tgo, although in some cases it probably does (Norris & Bonnell, 2006). The properties of 
this soft tail and its implications for the identification of SHBs are discussed later (§2.2.2 and §2.7). 

2.1.2. Temporal structure 

Only a small fraction of the BATSE SHB sample have sufficient signal to noise ratio (S/N) to conduct 
analysis of the fine temporal structure. Nakar & Piran (2002b) analyze a sample of 33 bright SHBs at a 
resolution of 2 ms, looking for variability (separate pulses) within the bursts. The low S/N and the limited 
time resolution imply that the observed variability is only a lower limit of the true one. Yet, they find that 

^ http:/ /Ihoawww. gsfc.nasa.gov/docs/gamcosray/lcgr/konus/ 
^ http:/ /www.ssl.berkeley.cdu/ipnS/ 
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Table 1 

Short vs. Long GRBs 





Short GRBs 


Long GRBs 


Section* 


General 


BATSE observed all sky rate 


Ri 170yr-i 


500 yr-i 


§2.6 


BATSE observed local rate density 


~ 10 Gpc-=* yr-i 


~ 0.5 Gpc-=* yr-i 


§2.6 


Host Galaxy types 


Early & Late 


Late 


§2.4 


Host specific SFR 


<lMo/yr/(L/L.) 


~ 10Mo/yr/(L/L.) 


§2.4 


Median observed redshift 


0.25 


R! 2.5 


§2.3 


Supernovae association 


No 


Yes (at least some) 


§2.5 


Progenitor 


NS-BH/NS-NS/? 


Massive star 


§5 


Prompt emission 


Typical BATSE duration 


Ri 0.8 s 


Ri 30 s 


§2.1.1 


Best fit spectral model'' 


Power-law + exp cutoff 


Band function 


§2.1.3 


. 

7, ISO 


10-19 _ iqSI erg 


1052 _ 1054 erg 


§2.1.4 


Lb . 

7, ISO 


1050 _ 1052 erg/s 


10^0 - 10^2 erg/s 


§2.1.4 


Average observed flux '• 


- 5 ■ 10-1° GcV/cm^/s/sr ' 


10-** GcV/cm^/s/sr 




Local energy output rate 


~ lO^'^ erg Gpc"'' yr^^ ^ 


10^^ erg Gpc--^ yr-^ 




Total energy output ° 


~ lO''^ erg 


~ 10^6 erg 





X-Ray afterglow 



X-ray dark 


Some bursts 


None 


§2.2.1 & §4.5 


Typical decay phase (~ t~^) 


Yes 


Yes 


§2.2.1 & §4.1 


Shallow decay phase (~ t-'^-^^) 


Not observed yet 


Yes 


§2.2.1 


X-ray flares 


Yes 


Yes 


§2.2.1 & §4.7 



Fireball model^ 



Lorcntz factor 


>30 


> 100 


§3.2 


Prompt emission energy source 


Internal processes^ 


Internal processes^ 


§3.4.1 


Afterglow energy source 


External shock 


External shock 


§4.1 


Synchrotron self-Compton 


Y « 1 


y >i 


§4.2 


Early afterglow & prompt emission 


Well separated 


May overlap 


§4.3 


Early afterglow (baryonic flow) 


Optically faint 


Optically bright 


§4.3 


Gamma-ray efficiency 


~ 10% 


>10% 


§4.5 


X-ray dark afterglows'' 


n < 10-^ cm-'' 




§4.5 


Average beaming 


1 < /iT' <ioo 


ft' ~ 75 


§4.6 



* The section in which the SHB table entry is discussed. Highlighted sections present new results that appear in this review 
for the first time. 

The spectral model that provides the best fit to most of the bursts (may be affected by observational selection effects) . 
b Isotropic equivalent quantities per burst. 

The total observed ~MeV 7-ray flux from GRBs, averaged over long time (3> 1 day). 

The energy output density rate in ~MeV 7-rays of the entire local GRB population. 

The total energy output in ~MeV 7-rays of the entire GRB population in the observed universe. 
' Theoretical interpretations in a framework of models where outflow of any type (e.g., baryonic, Poynting-flux-dominated) 
expands quasi-spherically and interacts with the ambient medium by driving a blast wave into it. 

s In Poynting-flux-dominated outflow models the dissipation gf the internal magnetic field that leads to the prompt emission 
may be induced by interaction with the external medium. 

b The most likely, but not the only, explanation of X-ray dark afterglows (see §4.5 for details). 



Tan [sec] 



Fig. 1. The bimodal duration distribution of GRBs. The observations (2041 bursts in the current BATSE catalog) are marked 
by the thick stairs. The decomposition of the distribution into two lognormal distributions, as determined by Horvath (2002), 
{thin solid lines) and the sum of these components {thick solid line) are superposed. 

most of the bursts in the sample exhibit variabiUty on time scales that are shorter than the bursts' durations 
(see Fig. 2 for a typical SHB light curve). More than half of the bursts in their sample show at least two 
well-separated pulses and more than a third show rapid variability in the sense that the shortest pulse is 
shorter by more than an order of magnitude than the burst duration. No correlation is found between the 
duration of a burst and the duration of its sub-pulses (given that the burst is not single pulsed). The duration 
distribution of single pulses ranges from 5 ms to 300 ms with a broad peak around 50 ms. Thus, the lower 
limit on the shortest time scale observed in these SHBs is of the order of 10 ms, and is set by the resolution 
limit. Shorter time scales are probably present, as evident from a single case in which a very bright < 1 ms 
pulse is observed in a SHB (Scargle, Norris & Bonnell, 1998, fainter pulses than this one cannot be resolved 
on ms time scale). McBreen et al. (2001) analyzed the distribution of various temporal properties of pulses 
in 100 bright BATSE SHBs. They find that the rise times, fall times, FWHM, pulse amplitudes and areas 
are all consistent with lognormal distributions and that time intervals between pulses and pulse amplitudes 
are highly correlated with each other. 

A comparison of the temporal structure of bright SHBs to the initial 2 s of a sample of long GRBs ^ shows 
similar time scales and similar distributions of pulse durations (Nakar & Piran, 2002b). This similarity is 
demonstrated in figure 2. Similarly, McBreen et al. (2001) find a great similarity between the lognormal 
distributions and correlations in the temporal structure of short and long GRBs. On the other hand, an 
examination of the temporal evolution of pulses as a function of frequency shows a different behavior in 
long and short bursts. Norris, Scargle & Bonnell (2001) compare the spectral lags of short and long GRBs. 
Spectral lag is a measurement of the spectral evolution timescale of the pulse structure, where a positive 
value indicates a hard-to-soft evolution (see Norris, Marani & Bonnell, 2000, for an exact definition). They 
find that long bursts show positive spectral lags that extend up to ~ 2 s with a core around 50 ms. SHBs, 
however, show a symmetric distribution of lags that ranges between ±30 ms. 



^ The sample includes only long GRBs that have high resolution light curve and an initial pulse that is shorter than 2 s 
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Fig. 2. Left: The light curve of GRB 910718, a bright SHB with Tgo = 0.25 s. Right: The first 0.7 s of GRB 941228, a bright 
long GRB with Tgo = 62 s. The figure demonstrates the similarity in short-time-scale structure in long and short GRBs. The 
resolution of both light curves is 5 ms. Both bursts show variability down to the resolution limit. 



Thus the temporal structure of SHBs shows both similarities and dissimilarities to that of long GRBs. 
Unfortunately additional comparisons of the temporal structure of short and long bursts were not carried 
out so far. mostly because of the difficulty in the analysis of SHB light curves. While the temporal structure 
of long GRBs was explored in detail (e.g., Norris et al. , 1996), only the general properties of the temporal 
structure of SHBs are known. 



2.1.3. Spectrum 

Comparison of the hardness ratio of SHBs to that of long GRBs shows that SHBs are on average harder. 
This result was used by Kouveliotou et al. (1993), together with the bimodal duration distribution, to 
suggest that SHBs and long GRBs are two distinctive populations. Figure 3 depicts Tgo and the hardness 
(ratio between the 50 - 100 keV and the 25 - 50 keV fiuence) of BATSE GRBs. It is evident that while 
SHBs are on average harder, the variance in the hardness ratio is such that the two populations overlap and 
the distinction between them is not well determined. 

Ghirlanda, Ghiselhni & Celotti (2004) explore the time integrated spectra of 28 bright BATSE SHBs. 
They attempt to fit the spectrum of each burst with four different spectral models — a single power-law, 
a broken power-law. a Band function (a smoothly broken power-law; Band et al. , 1993) and a power-law 
with an exponential cut-off (PLE): 

dN 

— = iVoS"cxp [^E/E^] (1) 

where N is the photon count and Ep = {2 + a)Eo is the peak of i^Fiy. They find that PLE model (Eq. 1) 
provides the best fit to the data. The values of Ep that Ghirlanda, Ghisellini & Celotti (2004) find are in 
the range 50 — 1000 keV (mean value 355 ± 30 keV), which are comparable to values observed in long GRBs. 
They find values of a between —2 and 0.5 (mean value —0.58 ± 0.1), which are significantly harder than 
those of long GRBs. These higher values of the low energy power-law index are the main driver of the higher 
hardness ratio of SHBs. 
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Fig. 3. The duration Tgo and the hardness ratio of all GRBs with available data from the BATSE catalog (dots). Swift (stars) 
and HETE-2 (plus) SHBs are marked as well. The average logarithmic values of BATSE bursts with Tgo > 2 s and Tgo < 2 
s are marked by the two squares. The hardness ratio used here is the ratio of the fluence in the 50-100 keV and 25-50 keV 
energy bands. The BATSE data is taken from the current BATSE catalog and Swift and HETE-2 data are taken from Table 
2, where the hardness ratio is calculated using the photon index. Note that the fluence ratio that is used here is different than 
the typical hardness ratio, that uses counts. 

The Konus-Wind SHB catalog (Mazets et al. , 2004) contains 140 spectra of 98 SHBs (for some SHBs 
the spectrum is time integrated and for some it is given in two or three different time intervals). Mazets 
et al. (2004) use the same spectral models as Ghirlanda, Ghiscllini & Celotti (2004) to fit the spectra 
of Konus-Wind SHBs. They fit 60% of the spectra with a PLE, 21% with a single power-law and 16% 
with a Band function. The average a value of the PLE spectra is —0.78, rather similar to the one found 
by Ghirlanda, GhiscUini & Celotti (2004) for BATSE SHBs. The average Ep of PLE spectra is 837 keV, 
significantly higher than the value found for BATSE events. This difference is probably a result of the wider 
energy window of Konus-Wind, 15 - 10000 kcV compared to 25 - 2000 keV for BATSE. The highest Ep 
value in the Konus-Wind catalog is 3.1 McV. In about 30% of the sample photons are observed up to 5 MeV 
but only in « 5% of the bursts harder photons are observed. Only a single burst shows 10 MeV photons 
(the one with a PLE spectrum with Ep = 3.1 MeV). 

An important point is that while most spectra are consistent with exponential cut-offs, the low energy 
spectrum of most SHBs is too soft to be consistent with a blackbody spectrum. Lazzati, Ghirlanda & 
Ghiselhni (2005) studied the spectra of 76 BATSE SHBs and find that in more than 75% of the bursts the 
spectrum is inconsistent with a blackbody spectrum. 

The fact that PLE provides the best fit to most SHB spectra is in contrast to long GRBs, where a Band 
function, with a rather shallow high energy power-law (^ _E"^-^^), provides the best spectral fit for most 
of the bursts (Prcece et al. , 2000; Ghirlanda, Celotti & Ghisellini, 2002). This dissimilarity may reflect a 
genuine intrinsic difference between the spectra of short and long GRBs, in which case long GRBs are actually 
harder than SHBs at high energies (^MeV). However, this difference may also be a result of observational 
selection effects. A PLE fit may be preferred over a Band function fit in cases of low signal-to-noise ratio at 
high energies, as suggested by the results of Kaneko et al. (2006) that analyze the spectra of the 17 brightest 
BATSE SHBs and find that most of them are well fitted by a Band function or a broken power-law. 



Among SHBs with known redshifts there are only two bursts with good spectral data (the spectral range 
of Swift alone is too narrow to constrain the broad band spectrum). The prompt emission time- integrated 
spectrum of HETE-2 SHB 050709 is best fitted by a PLE with Ep = 86. 5^1^ keV and a = -0.821°;^^ 
{Ea ~ 73 keV; Villasenor et al. , 2005). The Konus-Wind spectrum of SHB 051221A is fitted by a PLE with 
Ep = 402^^^ jjgY ^jjj Q, ^ _i.08i°;ii^ {Eq w 436 keV; Golenetskii et al. , 2005a). 

As I discuss below, the hardest observed non-thermal photons in the prompt emission play an important 
role in constraining the Lorentz factor of the prompt emission source (§3.2). The observations described 
above show that there are many SHB spectra that are consistent with having no non-thermal photons above 
~ 1 MeV, while other bursts show non-thermal emission at least up to ^ 5 MeV. There is no published 
evidence that SHBs emit non-thermal photons above 10 MeV. This is in contrast to long GRBs were there 
are bursts that show > 100 MeV photons (e.g., Schneid et al. , 1992; Sommer et al. , 1994; Hurley et al. , 
1994). The lack of observed > 100 MeV photons from SHBs may be in part an observational selection effect 
as well. At very high energies the background is negligible on prompt emission time scales and therefore 
the detectability depends on the fluence and not on the flux. Hence, SHBs, which have a significantly lower 
fluence than long GRBs, are much harder to detect at these energies. 

2.1.4. Isotropic equivalent energy 

In the few cases where SHB redshifts are known, the isotropic equivalent gamma-ray energy, E^ iso, and 
luminosity, L^ igo., can be calculated (i.e., assuming that the source emission is constant over the whole An 
sr solid angle). These values are listed in Table 2. The isotropic equivalent energy of SHBs ranges between 
10^^ — 10^^ erg, which is 2 — 3 orders of magnitude smaller than the energy range of long GRBs. The peak 
isotropic equivalent luminosity in most SHBs with known redshift is 10^" erg/s and in a single case (SHB 
051221 A) it is even larger than 10'"'^ erg/s. These values arc comparable to those observed in long GRBs 
(see Table 2). 

Amati et al. (2002) have shown that there is a tight relation between E~^^iso and the redshift correct peak 
spectral energy, Ep{\ + z), of long GRBs with known redshifts (see also Lloyd, Petrosian & Mallozzi, 2000; 
Lloyd-Ronning & Ramirez-Ruiz, 2002; Lamb, Donaghy & Graziani, 2004). Even before redshifts of some 
SHBs were determined, it was shown that SHBs do not follow the Amati relation^ (Ghirlanda, Ghisellini 
& Celotti, 2004; Nakar & Piran, 2005). The recent detection of SHB redshifts has shown that SHBs with 
known redshifts do not follow the Amati relation as well (Amati, 2006). 

2.2. Afterglow 

The breakthrough in the study of SHBs came in the spring-summer of 2005 with the first detections of 
X-ray, optical and radio afterglow emission (Gehrels et al. , 2005; Castro-Tirado et al. , 2005; Prochaska et 
al. , 2005; Fox et al. , 2005; Hjorth et al. , 2005a,b; Bloom et al. , 2006; Covino et al. , 2006; Bcrgcr et al. , 
2005) . This discovery was facilitated by timely and accurate localizations by Swift and HETE-2 , and came 
8 years after the first afterglow detection of a long GRB. The difficulty in detecting SHB afterglows is to 
achieve precise localization from a small number of photons (compared to long GRBs) and then to quickly 
point a sensitive X-ray instrument in order to detect the afterglow, which is significantly fainter relative to 
long GRB afterglows. 

The importance of afterglow detection is twofold. First, afterglow detection enables sub-arcsecond local- 
ization and, possibly, unambiguous determination of the host galaxy and its redshift. Second, the afterglow 
itself can teach us about the processes that take place after the explosion and to provide clues about the 
physics of the central engine and about the properties and environment of the progenitor. Here I discuss 
the observational properties of the afterglows themselves, while the environmental properties are discussed 
later. The discussion below is based on the small sample of the SHBs with observed afterglows, and even 

* Nakar & Piran (2005) and Band & Preccc (2005) have used the same method to show that many BATSE long GRBs do not 
follow the Amati relation as well. 
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within this smah sample a large variance exists. It is likely that as more SHB afterglows will be observed, 
the discussion below will have to be updated or revised. 

2.2.1. The late time afterglow 

Figures 4a & 4b depict the X-ray and optical afterglow isotropic equivalent luminosities as a function 
of the source time of all SHBs with known redshift. The observed X-ray flux as a function of observer 
time is presented if Fig. 5. The light curves roughly show a power-law decay on scales of hours to days, 
with significant additional variability in some cases. After several hours, the X-ray and optical temporal 
power-law decay appear to be rather similar in all observed SHBs with indices in the range of a ~ 1 — 1.5 
{Fi, (x t~"u~^). An exception is the X-ray afterglow of SHB 051210 (La Parola et al. , 2006), which decays as 
^^2.57 (j^j-jj^g least the first several hours (the redshift of this burst is unknown and therefore its afterglow 
is not in the figure). The optical and X-ray afterglows show power-law spectra with indices in the range 
(3 ~ 0.5 — 1.5. These general properties of late SHB afterglows are similar to those observed in long bursts. 

Whenever a dense temporal sampling of the afterglow (mostly in X-rays) exists, significant superposed 
variability is evident, with diverse characteristics. The afterglow of SHB 051221 shows an X-ray brightening 
after ~ 2 hr, with a similar power-law decay before and after the brightening. A bright X-ray bump is 
observed, superimposed on a power-law decay, between 3 hr and 5 days after SHB 050724 occurred. In both 
of these cases no high-amplitude, rapid variability {AF/F ^ 1 and At/t ^ 1) is seen. SHB 050709 does 
show rapid variability; after w 15 days, a bright X-ray spike is observed {AF/F w 10) that decays within 
2 hr {At/t w 0.01). However, the detection of this spike is based on a single Chandra observation of nine 
photons (Fox et al. , 2005). This kind of variability, if confirmed, has far-reaching implications on theoretical 
models, as discussed in §4. 

Observed SHB afterglows are significantly fainter than those of long GRBs. For example, after 1 day the 
SHB X-ray flux is fainter by more than an order of magnitude {^ 5 • lO^^"' erg/cm^/s in the four detected 
SHB late afterglows, compared to ~ 5 • 10~^^ erg/cm^/s in long GRBs^ ). Moreover, early {t ~ 100 s) X-ray 
afterglow emission was detected in every long GRB that was observed at early times by Swift , while three 
out of the 11 SHBs observed by the Swift X-ray telescope (XRT) within ^ 100 s, lack detectable X-ray 
emission. If SHB afterglow flux is related to the fiux of the prompt emission, as observed in long GRBs, 
then SHB afterglows are expected to be fainter (see §2.1.4). Fig. 4c shows that indeed this is part of the 
reason for the faintness of SHB afterglows. The figure presents SHB X-ray light curves multiplied by the time 
since the burst (illustrating the total energy emitted in the X-ray afterglow) , as well as the prompt emission 
fluence (plotted at t = 1 s). This figure clearly shows that the X-ray afterglow is strongly correlated with 
the gamma-ray fluence and that the energy emitted in the X-ray afterglow is smaller than the gamma-ray 
energy by no more than a factor of ^ 10. However, Fig. 6 shows that in SHBs in which the X-ray afterglow 
is not detected, the upper limits on the afterglows indicate that in these cases the emission is very faint, 
even compared to the gamma-ray fluence. See §4.5 for a discussion of the implications of this result. 

We conclude that late SHB afterglows (hours- weeks) show roughly similar temporal and spectral properties 
to those observed in long GRBs, while being fainter on average. The afterglows of some SHBs are considerably 
fainter than these of long GRBs, even after correcting for the lower 7-ray energy output of these SHBs. The 
small number of SHB afterglows prevents a more detailed comparison. These observations are discussed in 
the context of the external shock model for GRB afterglows in §4. 

2.2.2. Soft X-ray "tails" of the prompt emission and the early afterglow 

Several SHBs show a bright X-ray tail that follows the short prompt gamma-ray emission and lasts 
for 100 s. This X-ray component is most prominent in SHBs 050709 and 050724, where its energy is 
comparable to (SHB 050724) or even larger by a factor of '-^ 3 (SHB 050709) than the energy in the prompt 
gamma-rays (Barthelmy et al. , 2005; Villasenor et al. , 2005). This early luminous X-ray emission seems to 
be a common feature in SHBs, although there are cases in which the X-ray tail is either much weaker or does 
not exist at all (e.g. SHBs 050906 and 050925). An early hint to the existence of this component was found 



^ Values are taken from the Swift archive, http://swift.gsfc.nasa.gov/docs/swift/archivc/grb_table. 
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Prompt emission and afterglow properties 
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Various properties of SHB prompt and afterglow emission as detected by Swift and HETE-2 . All the quantities are 
in c.g.s units. 

" T90 as measured from Swift observations in the 15 — 350 keV energy band. This value is not strictly equivalent to the BATSE 
Tgo, where the 20 — 2000 keV energy band is used. 

Fluence [x lO-'^'erg/cm^] in the 15 - 350 keV energy band (except for 060801, 061201 & 061217 - 15 - 150 keV). The fluence 
in the 20-2000 keV energy band is given in parentheses for bursts that are also detected by Konus- Wind . 

Energy [xlO*^ erg] in the source frame 15 — 350 (20 — 2000) keV energy band. 

Peak isotropic equivalent luminosity [x lO^'' erg/s] in the source frame 15 — 350 (20 — 2000) keV energy band. This luminosity 
is measured over the source frame time interval that is indicated in the brackets. 

The inverse beaming factor, = ^'j/2, assuming a two sided 'top hat' jet with a half opening angle 8j. 
f The beaming corrected 7-ray energy: E^ = fbE~i^ii,o- 

9 The early X-ray flux [x lO'^erg/cm^/s] in the 0.3 - 10 keV energy band a.t t ^ 100 s (observer frame). 

Detection of afterglow in: XE - early (t ~ 100 s) X-ray, XL - late (t ~ 1 d) X-ray, O - optical, R-radio. 
* This burst was detected by HETE-2 , so Tgo, the fluence and the energy are given in the 30 — 400 keV energy band. 
** Tgo when calculated by simulating the BATSE algorithm. The maximal beaming factor and the minimal energy in the 
parenthesis are derived with the more conservative limit on the jet opening angle of SHB 050724(see §4.6). 
t The galactic latitude (0.1) and the spectrum of this burst arc consistent with an origin of a Galactic SGR. However, there is 
no known Galactic SGR at this location. 

t These bursts were observed while the manuscript was in the refereeing process and are therefore not discussed in the text. 
(?) Based on a low significance association (Ri 90% confidence). 



References - 1. Gehrels et al. (2005); Bloom et al. (2006) 2. Villasenor et al. (2005); Fox et al. (2005); Hjorth et 
al. (2005a) 3. Barthelmy et al. (2005); Berger et al. (2005); Krimm et al. (2005); Grupe et al. (2006) 4. Sato et al. (2005); 
Prochaska et al. (2005); Berger (2006) 5. Parsons et al. (2005) 6. Markwardt et al. (2005) 7. Barbier et al. (2005) 8. La 
Parola et al. (2006) 9. Cummings et al. (2005); Soderberg et al. (2006a); Burrows et al. (2006) 10. Roming et al. (2006) 
11. Bloom et al. (2007); Sato et al. (2006a) 12. Sato et al. (2006b) 13. Markwardt et al. (2006) 14. Parsons et al. (2006) 
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Fig. 4. The isotropic equivalent optical and X-ray luminosity as a function of the source frame time. Plotted are all SHB 
afterglows with a known redshift. (a) Optical luminosity (there was no optical detection of SHB 050509B). (b) X-ray luminosity, 
(c) X-ray luminosity multiplied by the time since the burst. This quantity illustrates the total energy emitted in the X-ray 
afterglow at a given time. The isotropic equivalent prompt emission fluence is plotted at t = 1 s (solid symbols) . References: 
SHB 050509B: Gehrels et al. (2005); SHB 050709: Villasenor et al. (2005); Fox et al. (2005); Hjorth et al. (2005a); Watson 
et al. (2006); SHB 050724: Barthelmy et al. (2005); Grupe et al. (2006); Berger et al. (2005); SHB 051221A: Soderberg et 
al. (2006a); Burrows et al. (2006) 

by Lazzati, Ramirez-Ruiz & Ghisellini (2001) that summed the gamma-ray hght curves of 76 BATSE SHBs 
and detected an X-ray signal starting 30 s after the trigger and lasting for ^ 100 s (see also Connaughton, 
2002; Frederiks et al. , 2004). The longer time scale and the softer spectrum of the X-ray tail, as well as 
the temporal gap from the short-hard prompt emission, clearly distinguishes it from the prompt emission. 
However, extrapolation of the late afterglow back to early times suggests that the X-ray tail is not the onset 
of the late afterglow. It is therefore unclear if the physical origin of this component is related to the prompt 
emission, to the afterglow or to a third process. In those cases where the energy in the X-ray tail is large, it 
may be detected also in gamma-rays, causing a SHB to look like a long GRB (see §2.7). Norris & Bonnell 
(2006) find eight BATSE GRBs with Tgo > 2 s that show an initial short duration emission that is followed 
by a long episode of softer emission. The spectral lags of all these bursts are short, suggesting that these are 
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Fig. 5. SHB X-ray afterglow light curves observed by the Swift X-ray telescope (XRT; Windowed Timing mode in blue and 
Photon Counting mode in red). The light curves are given in units of count/s. The conversion factor to energy flux units is 
roughly 1 count/s 5 X 10~^^crg/cm^/s, where the exact conversion factor depends on the X-ray spectrum of each afterglow. 
Courtesy of Judith Racusin and David Burrows. 



SHBs. The extended emission in some of these bursts is highly variable and its energy is comparable to that 
of the initial spike. Based on comparison between their sample and the results of Lazzati, Ramirez-Ruiz &: 
GhiscUini (2001), they suggest that the peak flux ratio between the initial spike and the extended tail can 
vary by a factor 10^ between different bursts. 

In bursts that exhibit a bright X-ray tail, it is followed by a steep temporal decay (a > 2; a = 
—dlog{Fi,)/dlog{t)) before the typical late afterglow (a ~ 1) is observed. This fast decay is clearly seen 
in SHBs 050709 and 050724. A fast decay phase is also observed in SHB 051210 (La Parola et al. , 2006). 
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Fig. 6. A histogram of the ratio between the X-ray energy flux at time t multiplied by t and the prompt gamma-ray fluence. This 
is an estimate of the ratio between the energy emitted in the late X-ray afterglow and in the prompt emission. The ratio is given 
for Swift long bursts {thin line) and for SHBs with X-ray afterglow observed after ~ 1 day [thick line). The upper limit at 4 ~ 100 
s for two Swift SHBs without detected X-ray afterglow is marked with arrow {thick line + arrow). The theoretical implications 
of this figure are discussed in §4.5. Reference: Swift archive, http://swift.gsfc.nasa.gov/docs/swift/archive/grb_table. 

Here the first X-ray observation is after 80 s and the flux is aheady faUing rapidly, until a flare is observed 
between 100 — 200 s. The energy in this X-ray flare is comparable to the prompt emission energy. After the 
flare the X-ray flux decays as t~^-^ and there is no evidence for a regular late afterglow decay. In other bursts 
the regular X-ray decay (a ~ 1) is observed starting at early times (e.g., SHB 050509B). Fig. 7 illustrates 
the main features observed in SHB X-ray afterglows. 

The overall early X-ray emission of SHBs is different from the complex early X-ray afterglows of long 
GRBs (e.g., Nousek et al. 2006). Long GRBs do not show similar X-ray tails, which are separated from the 
prompt emission ^'^ , and while there is a typical phase of fast X-ray decay in long GRBs, this decay marks 
the end of the prompt emission and is not the end of a distinctive X-ray component as it is in SHBs. The 
phase of shallow X-ray decay, which is often observed in long GRBs, was not observed so far in SHBs. Fig. 
7 presents a qualitative comparison between the afterglows of long and short GRBs. 

An early optical afterglow was detected so far only in one SHB (060313; Roming et al. 2006). In this 
burst the optical emission is constant at first (t = 100 — 1000 s) followed by a shallow decay, superimposed 
with rapid variability. The X-ray afterglow of this burst is uncorrelated with the optical emission. 

2.3. Observed redshift distribution 



It was suspected that SHBs take place at eosmological distances as soon as it was realized that they may 
be a distinct phenomenon from long GRBs. Indirect evidence for their eosmological origin were the nearly 
isotropic angular distribution in the sky (Briggs et al. , 1996; Balazs, Meszaros & Horvath, 1998) and the 
value of {V/Vmax) < 0.5 (Katz & Canel, 1996; Schmidt, 2001). These indicators also suggested that 



Long GRBs do show X-ray flares but only rarely these have an energy that is comparable to that of the prompt emission 
V/Vmax is the ratio of the volume that is enclosed within the distance in which the event is actually observed, and the one 
which is enclosed within the maximal distance to which the same event would still be detectable, assuming an Euclidian space. 
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Fig. 7. An illustration of the main features observed in SHB X- and 7-ray light curves (thick black line) and comparison to 
long GRBs (thin red line). The dashed lines in the SHB light curve represent alternative features that are observed in some 
bursts. The hard SHB prompt emission is followed in some bursts by a soft X-ray tail. The tail is well separated both in time 
and in spectrum from the prompt emission. If a 'tail' is observed, it is followed by a steep decay that in most bursts turns into 
a shallower ~ t^^ decay. Late X-ray flares can be observed at any time. So far a single SHB (051221A) has shown a late break 
in its light curve {~ t~'^) that is interpreted as a jet break. Prompt emission of long GRBs, in contrast, is typically followed by 
a fast decay that flattens into a very shallow decay, before it steepens again to a ~ decay. X-ray flares and late light curve 
breaks are observed in long GRBs as well. 

observed SHBs arc closer than observed long GRBs. Guetta & Piran (2005) find that for SHBs (V/Vmax) = 
0.39 ± 0.02. which is significantly larger than the value {V/Vmax) = 0.29 ± 0.01 they find for long GRBs. 
Magliocchctti, Ghirlanda & Celotti (2003) calculate the two-point angular correlation function for BATSE 
SHBs and find a 2cr deviation from isotropy on angular scales ~ 2° — 4° (a similar deviation is not 
found for long GRBs). They suggest that this deviation is induced by nearby large-scale structure through 
a low-redshift population of SHBs. 

The first limits on the distance to individual SHBs was obtained by examination of Interplanetary Network 
(IPN) localizations with small error boxes (Nakar et al. , 2006). The lack of bright galaxies inside the 
small IPN error boxes of six bright SHBs implied that the distances to these bursts exceed 100 Mpc (this 
limit was derived under the assumption that the SHB rate follows UV, blue or red light). The limits on 
the distances translate to energy output lower limits of order IC*^ erg, implying that these bursts can be 



Under the assumption that sources are distributed uniformly in distance, {V/Vmax) = 0.5 in Euclidian space. A smaller value, 
therefore, suggests sources at cosmological distances where the smaller the value of (V/Vmax) , the larger is the typical distance 
to the sources. 

Earlier examination of IPN error boxes was carried out by Schaofcr et al. (1998), without deriving constraints on the bursts' 
distances. See also Schacfcr (2006). 
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Fig. 8. The cumulative distribution of observed SHB redshifts (including HETE-2 SHB 050709 and not including the low 
significance redshift association of SHB 060502B) [thick line] compared to the of cumulative distribution of Swift long GRBs 
with known redshifts [thin line]. The dashed thick linos corresponds to the two possible values of the redshift of SHB 050813 
(0.72 or 1.8). 

detected by BATSE at distances greater than Gpc. This result provided another strong indication that 
SHBs constitute a bone fide cosmological population. 

The final confirmation of the cosmological origin of SHBs came with the secure identification of their host 
galaxies, following sub-arcsecond localizations of afterglows detected by Swift and HETE-2 . The measured 
redshifts of the host galaxies are listed in Table 2. All together there are 3 SHBs (050709, 050724 & 051221) 
with sub-arcsecond localizations that fall on top of galaxies which are considered to be "securely" identified 
as the hosts. The afterglow of SHB 050509B was detected only by the Swift X-ray telescope (XRT) and 
therefore was localized only up to within several arcscconds. The error circle is located in the outskirts of a 
giant elliptical. The confidence that this is indeed the host galaxy is estimated to be 3 — 4(t (Gehrels et al. , 
2005; Bloom et al. , 2006). The XRT error circle of SHB 060502B is w 20" from a bright early type galaxy 
at z=0.287. Bloom et al. (2007) suggest that this galaxy is the host of SHB 060502B and they estimate 
the confidence of this association to be w 90%. Finally, the afterglow of SHB 050813 was also detected only 
by the XRT and its location falls on the outskirts of an early type galaxy that is part of a group/cluster 
at z = 0.72 (Prochaska et al. , 2005) and on top of another, fainter, early type galaxy that is part of an 
apparent cluster at z sa 1.8 (Berger, 2006). Therefore, this burst is most likely at redshift 0.72 or 1.8 and its 
host is an early type galaxy in a cluster or a group. 

The cumulative distribution of the observed SHB redshifts (including the one detected by HETE-2 ) is 
compared in figure 8 to the distribution of Swift long GRBs. The most striking result of this comparison 
is that SHBs are detected at much lower redshifts than typical long GRBs - while the median redshift of 
SHBs is z « 0.25, the median redshift of Swift long GRBs is z w 2.5. Note that some of the observed SHBs 
are bright enough to be detected out to a much higher redshift. For example the prompt emission as well as 
the afterglow of SHB 051221 would be detected by Swift also at a redshift of z w 2. Nevertheless, caution 
is required when using the observed redshift distribution to draw quantitative conclusions. The redshifts of 
more than half of the Swift SHBs are unknown and thus the observed distribution is likely to be biased by 
selection effects that are needed to be quantified before the entire sample can be used. 

Gal- Yam et al. (2005) study the four best-localized IPN SHBs (< 10 arcmin^), searching for significant 
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luminosity overdensities inside the error boxes. They identify two putative host/cluster associations with 
SHBs. SHB 790613 is found to be significantly (3(t c.l.) associated with the rich galaxy cluster Abell 1892 
{z = 0.09) while SHB 000607 is found to be associated (at the 2cr c.l.) with a bright galaxy at z = 0.14. 
They derive lower limits on the redshifts of the two other bursts they examine of 0.25[0.06] at l[2]a c.l. 

Further clues about the redshift of SHBs are found by correlating their locations with volume limited 
samples of galaxies and clusters. Tanvir et al. (2005) find a correlation between the locations of BATSE 
SHBs and the positions of galaxies in the local Universe, suggesting that between 10 and 25 percent of 
BATSE SHBs originate at very low redshifts (z < 0.025). These SHBs may be extra-galactic versions of 
the giant flares from SGR 1806-20 (see §2.9.1) although they show a stronger correlation with galaxies with 
Hubble type earlier than Sc. So far, none of the well localized SHBs (« 15 bursts) are associated with such 
a nearby galaxy. Ghirlanda et al. (2006) find a positive 2(7 angular cross-correlation signal, on scales < 3°, 
between SHBs and galaxy clusters at z < 0.45. The significance of the signal increases when only clusters 
at z < 0.1 are considered. 

Berger et al. (2006b) explore the XRT error circles of SHBs with detected early X-ray afterglows but 
without identified hosts. They find that the galaxies within these error circles are fainter than 23-26 mag, 
and therefore are most likely at higher redshifts than the identified SHB hosts. Based on comparison to the 
redshifts of galaxies in the GOODS (Cowie et al. , 2004; Wirth et al. , 2004) and the HUDF (Coe et al. , 
2006) Berger et al. (2006b) conclude that at least 1/4 of the Swift SHBs are at z > 0.7. This conclusion is 
valid if the progenitors of these SHBs did not travel long distances from their hosts so these are still within 
the XRT error circles. 

We conclude that SHBs were shown to be genuine cosmological events which are observed, on average, 
at significantly lower redshifts than long GRBs. Indirect methods suggest that a non-negligible fraction of 
SHBs originate even closer - in the local universe, while non-detection of bright galaxies in some XRT error 
circles suggest that a significant fraction of the SHBs originate at z ~ 1. These results still require validation 
by future observations. The conclusions that can be drawn from the observed redshift distribution on the 
intrinsic redshift distribution of SHBs are discussed in §5.1.2. 

2.4. Host galaxies and cluster associations 

Various properties of known SHB host galaxies are listed in Table 3 (see Berger, 2006; Bloom & Prochaska, 
2006, for recent SHB host galaxies brief reviews). Host galaxies include both early- and late-type galaxies, as 
well as field and cluster galaxies. This diversity is in striking contrast to the properties of the host galaxies 
of long GRBs, which are typically dwarf starburst galaxies and are always (so far) found in the field. 
Moreover, the typical specific star-formation rate (SFR) in a long GRB host is w 10 Mq yr~^{L/Lt)~^ 
(Christensen, Hjorth & Gorosabel, 2004) while the specific SFR of SHB late-type hosts is significantly 
smaller, < 1 A/© yr^^(L/_L,)^^. A conservative comparison between the typical stellar ages in the hosts of 
long and short GRBs rejects the null hypothesis that the two populations are drawn from the same parent 
distribution at 99.75% confidence (Gorosabel et al. , 2006). These differences between the host galaxies of 
short and long GRBs are among the main foundations of the recent validation that long and short GRBs 
are two different physical phenomena. 

SHBs that take place in early-type galaxies are clearly associated with an old stellar population. SHBs 
in late-type star-forming galaxies, on the other hand, are not necessarily associated with young stars. The 
hosts of both SHBs 050709 and 051221 show evidence for a significant population of old 1 Gyr) stars 
(Covino et al. , 2006; Soderberg et al. , 2006a). Moreover, an examination of the exact location of SHB 
050709 within its host shows that it is not associated with the main star-forming regions in that galaxy (Fox 
et al. , 2005). 

Berger et al. (2006a) explore the association of Swift SHBs with clusters. They determine, using optical 
spectroscopy, that SHBs 050709, 050724 and 051221A are not associated with galaxy clusters. Using X-ray 
observations they find that except for SHB 050509B there is no evidence for other associations of Swift SHBs 
with clusters that are brighter than 3 • 10^^^ erg/s/cm^, or with mass Af > 5 • 10^'^ Mq, assuming a typical 
redshift z — 0.3. These results suggest that about 20% of SHBs are associated with massive clusters, which 
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Table 3 

Putative host galaxy properties 



SHB 


z 


Type" 


L/L, 


SFR 


Offset 


Offset Assoc. 


SN'' 


Ref. 










Mo/yr 


(kpc) 






limit 




050509B 


0.225 


E{c) 


3 


< 0.1 




13+? 


3 - 4cr 


-13 


1 


050709 


0.16 


Sb/c 


0.1 


0.2 


3.8 


1.8 


Secure 


-12 


2 


050724 


0.258 


E/SO 


1.5 


< 0.03 


2.6 


0.4 


Secure 




3 


050813 


0.72 or 
1.8 


E/SO(c/g) 
E/SO(c) 














4 


051221A 


0.546 


Late 


0.3 


1.4 


0.76 


0.29 


Secure 


-17.2 


5 


060502B 


0.287(?) 


Early 


1.6 


0.6 


73 ± 19 




^ 90% 




6 


IPN SHBs; 


790613 


0.09 


E/SO(c) 










3(T 




7 


000607 


0.14 


Sb 


1 


0.3 






2cr 




7 



" Hubble type of the host galaxy. (c)[(g)]: The host belong to a galaxy cluster [group]. 

^ Early-type: the rest frame K-band host luminosity in units of K-band L« . Late-type: the rest frame B-band host luminosity 
in units of B-band L* . 

The confidence level of the association between the putative host galaxy and the SHB. The association is considered to be 
secured when a sub-arcsecond localization of the afterglow coincides with the putative host. 

''■ The limit on the rest-frame absolute magnitude after 7 - 14 days (M^ for SHBs 050509B and 050709 and My for SHB 
051221) of a supernova associated with the SHB. 

References - 1. Gehrels et al. (2005); Fox et al. (2005); Bloom et al. (2006) 2. Fox et al. (2005); Prochaska et al. 
(2005); Covino et al. (2006) 3. Berger et al. (2005) 4. Prochaska et al. (2005); Berger (2006) 5. Soderberg et al. (2006a) 6. 
Bloom et al. (2007) 7. Gal- Yam et al. (2005) 



is consistent, within the uncertainties, with the fraction of stellar mass in such clusters (Fukugita, Hogan & 
Peebles, 1998). 

Finally, it is important to remember that since host galaxies were identified only for less than half of the 
Swift bursts, this sample likely suffers from selection effects. The distribution of the SHB offsets from their 
host centers (Table 3) is also expected to be biased by similar selection effects. As we discuss later (§5.2.3), 
some progenitor models predict brighter afterglows within late-type galaxies, making such hosts easier to 
detect (Belczynski et al. , 2006). Additionally, the most popular progenitor model (compact binary merger) 
predicts that some of the events will take place far from their host galaxies (> 100 kpc). This idea is also 
supported by the large possible offset of SHB 050509B (see Table 3). In such cases the afterglow is expected 
to be very dim and even when an afterglow is detected, the association with the parent host can be unclear. 



2.5. Limits on a .supernova component 

Despite the low redshifts of the first two SHB afterglows and the extensive optical follow-up, no associated 
supernova was detected. The limits are so deep (see Table 3) that an object that is ~ 10 times fainter than 
the faintest observed supernova (of any type) would have been detected. These observations are again in 
striking contrast to the observations of long GRBs. So far a supernova was detected in almost any long GRB 
with a redshift low enough to enable such a detection (see a discussion of GRBs 060505 and 060614 in §2.8). 
The lack of associated supernovae is another line of evidence that support the distinct nature of long and 
short GRBs. 
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2.6. Observed rate 

The whole-sky detection rate of BATSE was w 170 SHBs per year (Meegan et al. , 1997). This rate can 
be translated to an observed local rate (z ^ 1) using the observed Swift rcdshift distribution. The first 
three Swift SHBs with known redshift constitutes a complete sample (see discussion in §5.1.1), unbiased by 
selection effects other than the Swift sensitivity, which is roughly comparable to that of BATSE for SHBs 
(Band, 2006). Two out of these three bursts are within a distance of 1 Gpc, implying that at least 15% of 
all Swift SHBs are within this distance (at the 2a c.l.). Since the SHB redshift distributions of Swift and 
BATSE bursts are not expected to be very different (given the similar detector thresholds) the observed 
local rate is TZsHB,obs ~ 10 Gpc'^ yr~^ (Nakar, Gal- Yam & Fox, 2006). This rate is larger by about an order 
of magnitude than estimated local rate of the long GRB (Guetta, Piran & Waxman, 2005). It is also higher 
than estimates of the SHB rate that were done before the detection of SHB afterglows (Schmidt, 2001; Ando, 
2004; Guetta & Piran, 2005) by a factor of 10 - 100. 

2.7. Identifying an SHB 

An important observational challenge is to be able to distinguish between SHBs and long GRBs on a 
burst-to-burst basis . This was not so important prior to the detection of afterglows, when SHBs were 
studied statistically, using mainly prompt emission observations. However, now, when most recent scientific 
progress is based on a small sample of bursts with afterglows, this topic becomes crucial. It will probably 
remain an important issue even when the sample size of SHBs with afterglows will increase significantly if, 
as happened for long GRBs, much of the progress will be based on a few bursts with unique properties (e.g., 
very low redshifts). 

Currently, no satisfactory tests can provide a definite burst classification based on the prompt emission 
properties alone. The test which is commonly used (Tgo < 2 s) is based on the BATSE duration distribution 
and, as evident from Fig. 1, is rough and somewhat arbitrary (Tgo < 3 s would serve just as well). Clearly, 
there are long GRBs with Tgo < 2 s and SHBs that last longer than 2 s. Usage of this criterion is further 
complicated by the possible contribution of the X-ray tail, which follows some SHBs, to Tgo (Norris & 
BonncU, 2006). As demonstrated in Fig. 3, considering the hardness ratio does not significantly improve the 
ability to make a correct identification. A third prompt emission property that can be used is the spectral 
lag (sec §2.1.2; Norris, Scarglc & Bonnell, 2001), which in long GRBs is typically positive and, on average, 
longer than in SHBs. Unfortunately, the lag distributions of the two populations overlap as well. While a 
burst with a long (> 0.1 s) positive lag probably belongs to the long group, the opposite is not true, since 
there are long GRBs with short lags. Here, the lag-luminosity anti-correlation that was found for long GRBs 
may be of help (Norris, Marani & Bonnell, 2000). Long GRBs with short spectral lags (~ 10 ms) tend to be 
the brightest ones (> 10^^ erg/s). Therefore, SHBs seem to occupy a separate region in lag-luminosity space 
- short lags (< 10 ms) and lower luminosity (< 10^^ crg/s). This method, although promising, still needs to 
be validated by increasing the sample size of SHBs with known luminosity. Finally, since there are outliers 
to the lag-luminosity relation of long GRBs, even this method cannot provide absolute classification. 

Tests that assign probability to the classification using the prompt emission alone (e.g., as suggested by 
Donaghy et al. 2006) are also questionable. First, the sample that is currently used to derive such tests 
is the BATSE sample, while the tests are applied to bursts detected by other detectors (e.g.. Swift and 
HETE-2 ). Such an application may be compromised since statistical properties of an observed sample 
depend strongly on detector properties. Second, any decomposition of the observed distribution to a sum of 
underlying distributions (e.g., Horvath, 2002) depends on the assumed shape of the underlying functions, 
which is not uniquely determined. For example, the X-ray tails observed in some SHBs suggest that the 
duration distribution of SHBs alone may be bimodal - one peak corresponds to SHBs with faint X-ray tails 
and the other incudes SHBs for which the X-ray tails are bright enough to affect Tgo (according to this 



^ Here, the term SHB refers to any burst that is physically associated with the population that constitutes the majority of 
BATSE GRBs with Tgo < 2 s. Note that such a burst may not have Tgo < 2 s. 
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speculation, the second longer peak is not observed because long GRBs dominate the observed distribution 
above '--^ 2 s). 

Based on the small sample of SHBs with afterglows it seems that a much better classification can be done 
based on environmental properties (Donaghy et al. , 2006), i.e., the host galaxy type, its star-formation 
rate, the location of the burst within the host (§2.4), and on the existence of an accompanying supernova 
(§2.5). This classification is also based on physical grounds — progenitors of long GRBs are associated with 
massive stars, while SHBs occur preferentially (according to the current small sample) in places where the 
probability to find massive stars is small. Put differently, if a burst with a duration of 100 s will be detected 
in an early type galaxy and without an accompanying supernova, it is most likely related physically to SHBs 
and not to long GRBs (see discussion about GRBs 060505 and 060614 in §2.8). 

A more quantitative test can be carried out when the afterglow is imaged by the Hubble Space Telescope 
(HST). Fruchter et al. (2006) explored 42 HST images of long GRB host galaxies and found that afterglow 
locations within the hosts are more concentrated in bright blue pixels than those of core-collapse supernovae. 
This concentration is most likely a result of the very short lifetime of the progenitors of long GRBs. A 
quantitative test to determine if a burst is long can be done by analyzing its HST image in the same way 
Fruchter et al. (2006) did, and comparing the result with the known distribution of long GRB locations. 

GRB 050416A demonstrates the difficulty to use the prompt emission properties to classify an event. 
It is a soft burst with a border-line duration (Tgg = 2.4 s) and a negative spectral lag (Sakamoto et al. , 
2006). As it turns out, this specific burst took place in the brightest blue HST pixel of a star-forming galaxy 
(« 4 M0/yr/(L/L*)) and its afterglow shows the probable signature of a supernova component (Soderberg 
et al. , 2006b), implying that it is most likely a long burst. 

In this review I used the "classical" criterion, namely Tgo < 2 s, when compiling Table 2 and as the 
sample of SHBs with afterglows or significant afterglow constraints. Most of the bursts in this sample 
are short enough, Tgo <C 1 s, to confidently assume that they are bone fide SHBs. The only two SHBs that 
have border-line durations and that have significant effect on theoretical interpretation are SHBs 050724 and 
051221A. The former is observed within an early-type galaxy and therefore is not a long GRB. SHB 051221A 
has a short spectral lag compared to long GRBs, even when its high luminosity is considered (Gehrels et 
al. , 2006), supporting its classification as a SHB. Nevertheless, the possibility that GRB 051221A is not a 
genuine SHB cannot be excluded. 

2.8. Additional putative SHBs 

Given the difficulty to identify SHBs, a growing number of GRBs that do not pass the rough cut of Tgo < 2 
s were suggested, for various reasons, to physically belong to the SHB population. Table 4 lists all such Swift 
GRBs as wen as HETE-2 GRB 060121. 

The most interesting bursts in this sample are GRBs 060121, 060505 and 060614. Each one of these 
bursts affects the SHB theory if it is related physically to this population. Unfortunately, it is unlikely that 
an unambiguous classification will ever be achieved in any of the cases. I mention some of the relevant 
theoretical implications of these bursts being SHBs in the theoretical sections. 

GRB 060121 is a border line burst (Tgg = 1.97 ± 0.06 s; Donaghy et al. , 2006) with an optical afterglow 
(Levan et al. , 2006a) and a faint host with a photometric redshift probability with two peaks z = 1.7 ± 0.4 
and z = 4.6 ± 0.5 (de Ugarte Postigo et al. , 2006). The source rest-frame duration is therefore less than 1 
s. This burst also has short spectral lag, however, this cannot be used here to support the suggestion that 
it is a SHB, since short lags are expected in long bursts with similar peak flux (w 3 • 10^'^ erg/s at z=1.7). 

GRB 060505 is another border-line GRB (Tgo = 4 ± Is) which is located on top of a highly star forming 
knot in the spiral arm of a galaxy at z = 0.089 (Fynbo et al. , 2006; Ofek et al. , 2007). This burst does 
not show any evidence of supernova emission down to a very strict limit of Mb > —12.6 and therefore was 
suggested as putative SHB. 

The duration of GRB 060614 is 102 s putting it, as long as duration is concerned, safely in the long GRB 
group (Gehrels et al. , 2006). However, Gal- Yam et al. (2006) do not find any supernova emission down to 
Mb > -12 in HST images of this z=0.125 burst (see also Fynbo et al. , 2006; Delia Valle et al. , 2006). 
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Table 4 

GRBs with Tgo 2 s that were proposed to be SHBs 



SHB 



05091lt 16 
051114A 2.2 
051211A 4.2 
051227 8 
06012ltt 2 

060505 4 
060614 102 

061006 130 
061210 85 



)4g j,j,g 



Ref. 



0.165(?) 2{?) 



4.6[1.7] 24[4] X 10=* 
0.089 1.2 
0.125 80 



Tueller et al. (2005); Berger et al. (2006a) 
Sakamoto et al (2005) 
Kawai et al (2005) 
HuUingor et al (2005) 

dc Ugarte Postigo et al (2006); Levan et al (2006a) 
Fynbo et al (2006); Ofck et al (2007) 
Gal- Yam et al (2006); Fynbo et al (2006) 
Delia Vallc et al (2006); Gehrols et al (2006) 
Krimm et al (2006) 
Palmer et al (2006) 



A compilation of Swift GRBs that were appeared in the literature as putative SHBs although their duration is 2 s or longer. 

HETE-2 GRB 060121 is included as well since its redshift is constrained. 

" - Tgo as measured from Swift observations in the 15 — 350 keV energy band. 

t - The putative redshift is based on a possible association with the cluster EDCC 493 (Berger et al , 2006a). 

tt - A photometric redshift. Both ^ Ri 4.6 and z 1.7 are consistent with the observations (de Ugarte Postigo et al , 2006). 

Moreover, the star-formation rate of the host is low (< 1 MQ/yr/(L/L,)) and its lag-luminosity puts it away 
from the long GRB population and together with the other Swift SHBs. Additionally, its location in the 
HST image is on a faint pixel compare to the long GRB sample of Fruchter et al. (2006). Finally, the light 
curve of this burst is composed of an initial hard episode, which lasts ^ 5 s, that is followed by a variable soft 
emission that lasts « 100 s. Based on these observations Gal- Yam et al. (2006) and Gehrels et al. (2006) 
suggest that this might be a long duration burst that is associated physically with the SHB population. 



2.9. SHB interlopers 

It is possible that the SHB population defined according to its BATSE duration, contains sub-classes 
which are not physically associated with the dominant SHB population, or with long GRBs. Actually, the 
recent 7-ray giant flare from SGR 1806-20 implies that this is indeed the case. Such possible sub-classes are 
discuss below. 



2.9.1. Extra- galactic SGR giant flares 

Soft gamma-ray repeaters (SGRs) are compact sources of persistent X-ray emission and repeating bursts 
of soft gamma-rays (see Woods & Thompson, 2006, for a review), which are believed to be highly-magnetized 
young neutron stars (known as magnetars; e.g., Duncan & Thompson 1992; Paczyhski 1992). At infrequent 
intervals, these sources emit extreme flares of high energy radiation, releasing more than 10*** ergs in the 
form of gamma-rays alone. These "giant flares" arc characterized by a short (~ 0.1 sec), hard (~ 300 kcV) 
and very intense spike that is followed by a long (~ 300 sec) pulsating soft "tail" (not to be confused with 
the X-ray tail observed in some SHBs). The maximal observed isotropic equivalent luminosity of such giant 
flares is ^ 10**^ crg/s following the eruption of SGR 1806—20 on December 2004 (Hurley et al. , 2005; Palmer 
et al. , 2005, for as short review of this event see Taylor & Granot 2006). 

Only the initial short and hard spike of an extragalactic giant flare would have been detected by BATSE 
because its luminosity is much higher than that of the pulsating tail, in which case it would be classified 
as a SHB (Duncan, 2001; Eichler, 2002). With a luminosity of ~ 10*^ erg/s such flares can be detected by 
BATSE up to ~ 50 Mpc, suggesting that a significant fraction of the BATSE SHB sample is comprised of 
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similar extra-galactic flares (Dar, 2005a, b; Palmer et al. , 2005; Hurley et al. , 2005; Nakar et al. , 2006). This 
fraction, however, cannot be very large given that none of the current Swift SHBs took place in a nearby 
galaxy, and therefore they are all much more energetic than the giant flare from SGR 1806-20. Additional 
constraints on the rate of SGR giant flares was obtained by exploration of IPN error boxes (Nakar et al. 
, 2006; Ofek, 2007), the lack of a BATSE SHE overdensity toward the virgo cluster (Palmer et al. , 2005; 
Popov & Stern, 2006) and the small fraction of SHBs with a blackbody spectrum (Lazzati, Ghirlanda & 
Ghisellini, 2005, ; the spectrum of the SGR giant flare is consistent with a blackcody). All these methods 
yield an upper limit on the fraction of giant flares out of the BATSE SHE sample, fscR, while the single 
observed event puts a lower limit on this fraction. Having the most stringent and least assumption dependent 
result, Ofek (2007) finds that 0.01 < Jsgr < 0.14 at 95% confidence. 

Interestingly, Tanvir et al. (2005) find evidence that a local population (within sa 100 Mpc) of SHEs may 
constitute w 10% of the total SHE population. At first glance it looks as if these are the expected extra- 
galactic SGR Giant flares. However, Tanvir et al. (2005) find an even stronger correlation when galaxies 
with a Hubble type later than Sc are excluded from the sample. SGRs arc thought to be produced by core 
collapse of young stars and to be associated with star forming regions, as supported by the Galactic SGRs 
(e.g., Gaensler et al. , 2001). Therefore, the stronger correlation in the case that the latest type galaxies are 
excluded is surprising. Following these observations Levan et al. (2006b) discussed an alternative channel 
to produce SGRs, a merger of two white dwarfs, which should take place also in early-type galaxies (see also 
Thompson & Duncan, 1995; King, Pringle & Wickramasinghe, 2001). 

Distinguishing between an SGR giant flare and a SHE is not trivial, even if the location of the burst 
coincides with a nearby galaxy. An example is GRE 051103, localized by the IPN (Golenetskii et al. , 2005b) 
in the vicinity of M81 and M82 (« 3.6 Mpc from Earth). The error box includes one of the star forming 
arms of M81 (Ofek et al. , 2006; Frederiks et al. , 2007). Lacking additional observations, it is currently 
unknown whether this was an SGR giant flare in M81/82 or a background SHE. 

2.9.2. Very short GRBs 

Chne, Matthey & Otwinowski (1999, 2001) and Cline et al. (2005) suggest that SHEs with Tgo < 0.1 
s compose a sub-class that is distinct from the rest of the SHE population. The main evidence that they 
present are a value of (V/Vmax) ~ 0-5 and an anisotropic distribution on the sky of the very short GREs. 
Eoth properties suggest a nearby population and are different than those of longer SHEs. They also find 
that Konus-Wind very short GREs show on average harder spectrum above 3 MeV than longer SHEs. 
Cline, Matthey & Otwinowski (2001) and Cline et al. (2005) suggest that these events may be produced by 
primordial black hole evaporations or by nearby extra-galactic SGR giant flares. Note that these observations 
suggest that if indeed this is a subclass, very short GREs take place in the local universe. So far there are 
5 Swift SHEs with Tgo < 0.1 s, none of which are associated with a nearby galaxy. 

3. Relativistic outflows and the prompt emission 

The theory of long GRE prompt and afterglow emission was explored extensively. The generally accepted 
picture is that a stellar object undergoes a catastrophic event leading to the formation of a central engine. 
The engine rapidly releases energy in a compact region (~ 10^ — 10 cm). This energy is deposited in the 
form of radiation, heat and/or electromagnetic field and leads to the acceleration of an ultra-relativistic 
outflow. At first, this flow is optically thick and it "carries" the energy out, while cooling adiabatically, to 
large distances (> 10^'^ cm), where it becomes optically thin. From this point and on, dissipation of the 
flow's energy, most likely by internal processes (i.e., without interaction with the ambient medium), results 
in observed radiation. The dissipation heats electrons to highly relativistic velocities and these radiate in the 
presence of magnetic field and/or radiation field (e.g., synchrotron and inverse compton radiation) that is 
either advected by the flow from the source or generated by the dissipation process. This emission is observed 
as the prompt gamma-rays. At larger radii (10"'^^ — 10^^ cm) the energy that remains in the relativistic flow 
is transferred into the circum-burst medium, generating a decelerating blast wave. This blast wave, which 
propagates into the external medium, is the source of the afterglow emission. 



23 




Fig. 9. Illustration of the emission sites of four photons from a relativistic spherical shell with a width A that expands with a 
Lorentz factor T. These photons set three typical time scales (see the text). Photon 'A' is emitted from the origin at the time 
of the explosion. Photons 'BVC and 'D' are all emitted when the shell front is at radius R. 'B' is emitted from the front of 
the relativistic shell on the line-of-sight. 'C is emitted from the front of the shell at an angle 8. 'D' is emitted from the rear of 
the shell on the line-of-sight. 

This picture is based on the observed luminosity, time scales and spectral and temporal evolution of long 
GRB prompt and afterglow emission. Since most of these properties are common to long and short GRBs, 
the same model is applied also to SHBs and appears to give an adequate explanation to their observations . 

Numerous review articles describe GRB theory, focusing on the general model outlines described above 
(e.g., Piran, 1999, 2005b; Meszaros, 2006). I will review here only the main theoretical aspects of this model 
while quantitatively confronting it (in several cases for the first time) with SHB observations. 

3.1. Relativistic effects 

The ultra-relativistic velocity of the source of the prompt and afterglow emission directly affects their 
observational properties. I briefly describe here some of these effects, serving as a basis for further discussion. 

3.1.1. Time scales 

Consider a spherical shell with a width A that is expanding with a relativistic Lorentz factor F 3> 1 and 
emits photons at radius r = i? from points 'B', 'C' & 'D' (Fig. 9). Photon 'A' is emitted from the origin at 
the time of the explosion (when the shell is ejected; Fig. 9). The different arrival times of these four photons 
to the observer set three different time scales: 

Line-of-sight time: tAB is the time interval between two photons emitted by the shell front along the 
line-of-sight one when the shell is ejected, r = 0, and the second when the shell is at r = R: 

tios = (2) 

This time scale is also comparable to the time between any two photons that are emitted over the expansion 
time, i.e., the time over which the radius doubles (which is also the adiabatic cooling time). If the Lorentz 



This is true only for the physical processes that produce the prompt and afterglow emission. There are two separate branches 
of theoretical research of short and long GRBs that explore the progenitors and the processes the lead to the formation of the 
central engine (see §5). 
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factor of the emitting shell is not constant, but decays as F cx i? ^ (e.g., during SHB external shock g = 3/2) 
thcntios=R/{2[l + 2g]cr^). 

Angular time; Ibc — R9^ /2c is the time interval between two photons emitted at radius R, one along the 
line of sight and the other at an angle 6* <C 1. Light abberation from relativistic sources implies that most 
of the photons that arrive to the observer are from angles < 1/r. This sets the angular time over which 
photons that are emitted from the shell at radius R arrive to the observer: 

tang = 2^p2 ' 

Shell width time: The time scale that is set by the shell width is not affected by the relativistic motion: 

tA = A/c. (4) 

Note that if we set f = to the arrival time of photon 'A' then photons from radius R arrive to the 
observer between tios and ~ 2tios (tios ~ tang)- As I will discuss later, this property makes it difficult for a 
single spherical shell to produce high amplitude rapid variability in the light curve (§3.4.1). 

In the scenario discussed above the emitting material and the emission front are propagating at the same 
Lorentz factor F. This is not always the case. For example, according to the afterglow theory the emission is 
generated by a blast wave (i.e., emission front) that propagates at a lorentz factor that is larger by a factor 
of V2 than the fluid that is crossing the shock (i.e., emitting material). In such case the line-of-sight time 
(Eq. 2) depends on the emission front Lorentz factor while the angular time (Eq. 3) depends on the Lorentz 
factor of the emitting material. 

3.1.2. Causal connection and quasi- sphericity 

Consider a spherical shell that is expanding at a Lorentz factor F 3> 1, and a signal (e.g., a sound wave) 
that is propagating on its surface at velocity cPs (as measured in the local rest frame). Between the radius R 
and R + dR this signal propagates an angular distance dOg ~ (3sdR/{VR). Therefore, during the expansion 
time (i.e., the time over which the radius doubles) the signal propagates 9s ~ /3s /F (logarithmic terms are 
neglected). Since (is < 1 the angular size of a causally connected patch of the shell is ~ 1/F. Therefore, a 
relativistic shell is 'frozen' over angular scales that are larger than 1/F. Relativistic light abberation implies 
that most of the photons that are emitted from an angle > 1/F with respect to the line-of-sight as observed 
from the origin, do not reach the observer. The combination of the two effects implies that any event 
occurring at angles that are larger than 1 /F cannot be observed directly and are also causally disconnected 
from the region that can be seen. Thus, the observer is completely ignorant of anything outside of the 1/F 
angular region. 

While this property makes it difficult to study the angular structure of relativistic outflows, it makes the 
early stages of GRBs considerably simpler to model. As long as the angular structure of the flow varies over 
angles that are larger than 1/F, each part of the outflow behaves as if it was a part of a spherical flow with 
the local properties. Thus we can apply spherically symmetric models to the prompt emission and to the 
initial phases of the afterglow. Throughout the paper I will refer to this property as quasi-sphericity. 

3.1.3. High latitude emission 

Consider an inflnitesimally thin spherical shell that expands relativistically with a Lorentz factor F = 
(1 — /3^)~^/^ and radiates an arbitrarily short pulse at radius R. What would be the observed shape of 
this pulse? The first photons that arrive to the observer are emitted on the line-of-sight and therefore 
arrive at tios- At later times photons arrive from larger and larger angles with respect to the line-of-sight, 
0, such that any given observer time is mapped onto a given angle ('latitude'). The relativistic blueshift 
decreases with increasing latitude (and time) resulting in a decay of the energy flux. Additionally, due to 
the decreasing blueshift, a given observer frequency corresponds to an emitted rest frame frequency that 
increases with latitude. Therefore, the observed pulse shape at a given frequency window depends on the 
emitted spectrum. In the case of a power-law spectrum {Fu oc v~^) the observed flux depends on the Lorentz 
boost as Fi, cx cos 6) is the inverse of the Doppler factor. Using the relation 
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between the observer time and the emission angle, i — t(os = l—Rcos{6)/c, one obtains I? cx l + {t — tios)/tang 
in the hmit of F ^ 1, implying (Kumar & Panaitescu, 2000; Nakar & Piran, 2003a): 

F.cx 1 + . (5) 

At late times (t ^ tios and t ^ tang) the light curve decays as t~^^^^\ This relativistic high latitude 
emission is important since it predicts a prolonged observed flux, with decreasing peak frequency, even in 
the case of an arbitrary short burst of emission (e.g., a 'naked' afterglow, see §4.4; Kumar & Panaitescu, 
2000). It also constrains the amount of possible variability from a spherically symmetric source (Nakar & 
Piran, 2003a, see §4.7). 



3.2. The Lorentz factor of the outflow 



Perhaps the most prominent feature of GRBs is that they are (special) ultra-relativistic sources. This 
is a well-established result, relying on several independent indications, some of which are (almost) model 
independent. Here I discuss two methods, starting with the opacity constraint, the most robust and model 
independent method , where a lower limit on the Lorentz factor is set by requiring that the source of 
the prompt emission is optically thin. The second method provides a measurement of the Lorentz factor F 
during the early afterglow and can be applied when the onset of the afterglow is observed. It depends on the 
afterglow model and is valid only when the observed afterglow emission is produced by an external shock 
(see §4). 

These two methods (as well as several others) were used to carefully analyze only long GRB observations 
and therefore the Lorentz factor of SHBs was never properly constrained. Here I briefly present the derivation 
of the theoretical constraints and then apply them to the observations of SHBs. The main result is that, like 
long GRBs, SHBs are ultra-relativistic. However, while observations of long GRBs require F > 100 — 300, 
for most of the SHBs F > 10 — 50 is consistent with the observations. 



3.2.1. Opacity constraints 

The prompt emission of GRBs (long and short) is non-thermal (§2.1.3), implying that the source is 
optically thin to the observed photons. On the other hand, if a non-relativistic source is assumed, a calculation 
of the optical depth to Thomson scattering, rr, based on the enormous observed luminosity of MeV 7-rays, 
results in tt ~ 10^^ (Schmidt, 1978). This discrepancy was known as the compactness problem and at first 
was used to argue that GRBs are Galactic. This conflict was alleviated when it was realized that the source 
of the emission may be moving at relativistic velocities towards the observer (e.g., Guilbert, Fabian & Rees, 
1983; Piran & Shemi, 1993). The most comprehensive calculation of the opacity limit on the Lorentz factors 
of long bursts appears in Lithwick & Sari (2001). Here I carry out similar analysis, adapting it to the possibly 
different prompt emission spectra of SHBs. 

The non-thermal spectrum of GRBs implies that the source is optically thin to Thomson scattering on 
e~e"'" pairs ^® . An inevitable source for such pairs is the annihilation of photons with rest frame energy 
^p/i ™eC^, where is the electron mass. Therefore, the Thomson optical depth for a given pulse during 
the prompt emission phase is : 

orNphfie' t, > mec^) 
" ' 

where ctt is the Thomson cross-section, Nph is the total number of emitted photons within the pulse, 
/(^p/i ^ iTT-eC^) is the fraction of photons that create pairs and R is the radius of the source. Relativistic 



In the case of long GRBs the most direct evidence of a mildly relativistic motion during the afterglow phase are the resolved 
radio images of the afterglow of GRB 030329 (Taylor et al. , 2004, 2005). 
^® Opacity to 77 pair production provides less stringent constraints on SHB Lorentz factors. 

^'^Hcre I assume that ho source is moving directly toward the observer. If the source is moving at some angle with respect to 
the line of sight then the optical depth increases and so docs the lower limit on F. 
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motion of the source has two effects. First, it reduces the energy of the photons in the rest frame, thereby 
reducing /. Second, for a given observed time of the pulse, St, it increases the emission radius as i? '--^ cStT^. 
The time scales and the luminosities of individual pulses in long and short GRBs are similar (see §2.1.2), 
but the spectrum may be different. While the spectrum of most long GRBs is best described by a Band 
function (smoothly broken power-law; Preece et al. , 2000; Ghirlanda, Celotti & Ghisclhni, 2002), the best fits 
spectrum of most SHBs is a low energy power-law and an exponential cut-off (Ghirlanda, Ghisellini & Celotti, 
2004, PLE; see §2.1.3). While this spectral difference may be a result of observational selection effects, PLE 
spectrum should be used when conservatively deriving the lowest Lorentz factor that is consistent with all 
current observations. The reason is that PLE spectrum is consistent with all available data of most bursts 
and it is less constraining than spectra with high energy power-law. Note also that Gev photons, which 
support spectral high-energy power-law, were observed in several long GRBs (e.g., Schneid et al. , 1992; 
Sommer et al. , 1994; Hurley et al. , 1994) while there is no report in the literature of a photon harder than 
10 MeV that was observed from a SHB. Hopefully, GLAST observations will unambiguously determine the 
high energy spectra of SHBs, enabling more stringent limits on SHB Lorentz factors. 

Using a power-law spectrum with index a and an exponential cut-off at Eq (Eq. 1), Eq. 6 becomes: 



TT« 10 S^,-7dL^2s5t_^—-V "^exp - , , 



(7) 



where is the observed gamma-ray fluencc of the pulse, is the luminosity distance to the burst (at 
redshift z) and throughout the paper denotes A^/10^ in c.g.s units (except for Eq which is traditionally 
used in this context as the spectral cutoff energy). Requiring < 1 results in the following constraint on 
the Lorentz factor: 



^"'^'^ + (4 - a)ln(r) + In 



E. 



, , , ... , >30. (8) 

The logarithmic dependence on S'-y, 5t and is neglected in Eq. 8 (the range of the observed values of SHB 
pulses may affect the value of Eq. 8 by less than 50%). Figure 10 presents the lower limit on F as a function 
of Eq for three values of a. This lower limit is F > 15 for the majority of the bursts analyzed by Ghirlanda, 
Ghisellini & Celotti (2004), assuming that they are cosmological, while for SHBs 051221 and 050709 the 
opacity lower limits are F > 25 and F > 4 respectively. These lower limits are significantly lower than those 
obtained for long GRBs (F > 100; e.g., Lithwick & Sari, 2001). Note however that for both populations only 
lower limits on the Lorentz factor are available and, while the typical Lorentz factor of SHBs could be lower 
than that of long GRBs, a precise comparison between the real values of F is impossible. Additionally, the 
smaller lower limits on SHB Lorentz factors depend on the best-fit function of their spectra which might 
be affected by observational selection effects. Hopefully the high energy spectra of SHBs will be securely 
determined by the upcoming GLAST mission. 



3.2.2. Constraints from the onset of the afterglow 

Within the framework of all the models discussed in this review, the afterglow is produced by a blast 
wave that propagates into the circum-burst medium, which is generated by the interaction of the relativistic 
outflow with the ambient medium. Once most of the energy is deposited into the external medium the blast 
wave assumes a self-similar profile (Blandford & McKee, 1976) and the X-ray and optical emissions are 
expected to decay as power- laws (roughly as where t is the observer time since the trigger of the prompt 
emission; see §4). During the self similar phase: 

nt)^m[^f\^y\i^.f», (9) 

Where Ek.iso is the isotropic equivalent kinetic energy that remains in the relativistic outflow after the 
prompt emission phase and n is the density of the external medium, which is taken to be constant (as 
expected in the environment of SHBs; see §5). Note that the dependence on both Ek^iso and n is weak. 
Plugging the earliest observed time at which the afterglow shows a regular decay that flts the self-similar 
model into Eq. 9 provides a lower limit on the Lorentz factor of the outflow. 
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Fig. 10. The lower limit on the Lorentz factor of the prompt emission source as a function of the rest frame spectral typical 
energy -Eo(l + z) in units of rrieC^. The limit is derived by the opacity constraint (Eq. 8) using three different low energy 
power-law slopes a. The dots mark throe bursts, SHB 050709, SHB 051221A, and a typical SHB (according to Ghirlanda, 
Ghisellini & Celotti, 2004) at z=0.5. (prompt emission properties are taken from Villasonor et al. , 2005; Golenetskii et al. , 
2005a) 

Figure 4 shows three SHBs for which a regular afterglow, decaying as ~ t~^, is observed at t w 100 
s. Therefore, assuming that this emission results from an external shock, the Lorentz factor of the ejecta 
in these bursts is constrained to be F > 40 (e.g., Lee, Ramirez-Ruiz & Granot. 2005, use this method to 
constrain the Lorentz factor of SHB 050509B) . This lower limit is similar to the one obtained by the prompt 
emission opacity. Note that the Lorentz factor of SHB 050709 is not well constrained by this method, since 
the first observations of "regular" afterglow decay in this burst is more than a day after the burst. 

3.3. The composition of the relativistic outflow 

One of the main open questions about the physics of GRBs (long and short) as well as that of other 
relativistic astrophysical phenomena (e.g.. micro-quasars and active galactic nuclei) is how to launch (and 
often coUimate) a relativistic outflow. Most models of the "engine" that produces the flow are composed of 
different combinations of rotation, accretion and magnetic fields, where the energy source of the flow can 
be any one of the above. Whatever the energy source is. the outflow starts its way when a vast amount of 
energy is deposited in a relatively baryonic free compact environment. The evolution of the flow depends 
on the composition of the deposited energy (heat, magnetic field, etc.) and the environment in which the 
outflow propagates. Almost all suggested models are applicable for both long and short GRBs. where the 
main difference is the duration over which the engine is active, which is roughly comparable to the observed 
duration of the burst (§3.4.1). Another difference is that the engine of long GRBs is believed to operate 
at the center of a collapsing star, while in all current SHB progenitor models the engine is "exposed" (§5). 
Thus, relativistic ejecta from SHBs do not have to penetrate through several solar masses of surrounding 
material, and can be launched directly. 
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The relativistic outflow composition is currently unknown. The main candidates arc baryonic plasma 
(protons, electrons and likely neutrons) and magnetized plasma (at various ratios of magnetic to particle 
energy). Currently no observations of either long or short bursts conclusively point towards one of these 
cases. Since there are no studies of this topic that are specific for SHBs, and given that this topic is well 
covered elsewhere (Zhang & Meszaros, 2004; Piran, 2005b, a; Meszaros, 2006; Lyutikov, 2006), I provide here 
only a brief overview. 

3.3.1. Baryonic flow 

In this case the engine energy is deposited at the end of acceleration into an extremely small load of 
baryons, thereby accelerating the baryons to ultra-relativistic velocities. The most basic model of pure 
baryonic outflow suggests that the outflow is accelerated by radiation pressure. It was realized early on that 
releasing the observed luminosity in the form of pairs and/or radiation in a compact region (< 0.001 light 
sec) results in the relativistic expansion of a pairs-radiation fireball (Goodman, 1986; Paczynski, 1986). If 
there is a negligible load of baryons in such a fireball, most of the energy escapes as quasi-thermal radiation. 
However, if the proton load is high enough, then the accompanying electrons increase the optical depth 
and the fireball remains optically thick until almost all the the energy is converted into bulk motion kinetic 
energy (Shemi & Piran, 1990; Piran, Shemi & Narayan, 1993; Meszaros, Laguna & Rees, 1993; Grimsrud 
& Wasserman, 1998; Daigne & Mochkovitch, 2002b; Nakar, Piran & Sari, 2005; Li & Sari, 2006). The 
minimal (isotropic equivalent) proton load that results in full conversion of radiation to bulk motion energy 
corresponds to a maximal Lorentz factor: 

r,na. - ioooeI(,^r-^^^^tz^/\ (10) 

Where Rq is the radius in which the energy is released (comparable to the engine size) and T is the duration 
over which the energy, is isotropically deposited. Together with the lower limit on the Lorentz factor 
(~ 30; see §3.2) the allowed range of baryonic mass is {E ~ TMejC^): 

5 • w-HIqeII,''rI(^t'_{^ < M,, < 2 • 10-«Mq£;5o. (h) 

Note that while the mass lower limit (and the Lorentz factor upper limit) is model dependent and valid only 
in this scenario, the upper limit on the baryonic mass is generic since it is derived by the Lorentz factor 
opacity constraint (§3.2). This implies that in any model, the region in which the ejecta is accelerated must 
be almost completely clean of baryons. 

This simple model does not discuss collimation of the flow, which may be done by a hydrodynamic 
interaction with the environment or by magnetic forces. In both cases the interaction may affect the final 
Lorentz factor (as discussed below). Additionally, any variability of the outflow in this quasi-spherical model 
is generated by varying the deposition rate of the energy and/or the baryonic loading. If neutrons are 
present in the flreball as well, as expected if the engine is driven by accretion onto a compact object (e.g., 
Beloborodov, 2003b), these will be dragged by nucleon collisions along with the protons up to high Lorentz 
factor and may affect the evolution of the fireball at late times (e.g., Derishev, Kocharovsky & Kocharovsky, 
1999a,b; Pruet & Dalai, 2002; Beloborodov, 2003a,b; Vlahakis, Peng & Konigl, 2003a; Peng, Konigl & 
Granot, 2005; Rossi, Beloborodov & Rees, 2006). In this case the mass limit is roughly applicable to the 
total baryonic mass. 

An additional, purely hydrodynamic acceleration mechanism of a non-magnetized plasma flow was dis- 
cussed recently by Aloy & Rezzolla (2006) (see also Aloy, Janka & Miiller 2005). The acceleration is powered 
by rarefaction waves that propagate into a relativistic flow in the presence of large velocities tangential to a 
discontinuity. Such hydrodynamical conditions are expected to be generated by the interaction of the flow 
with a thick accretion torus that feeds the compact engine. Thus, once the flow achieves a moderate Lorentz 
factor (e.g., by radiation pressure as discussed above) additional acceleration is provided by this process. 
As a result, the maximal Lorentz factor in this model is not set by transparency considerations and can 
significantly exceed the upper limit in Eq. 11. Rarefaction waves also provide natural collimation to the flow 
and Aloy, Janka & Miiller (2005) find that for the observed SHB parameters, the resulting flow is a narrow 
jet 0.1 rad) with a core Lorentz factor of ~ 1000 and rather sharp edges. 
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3.3.2. Magnetized flow 

Many authors suggested models in which magnetic fields play an important role in the dynamics (accel- 
eration and possible coUimation) as well as in the final energetic content of the outflow (e.g., Usov, 1992; 
Levinson & Eichler, 1993; Usov, 1994; Thompson, 1994, 2006; Katz, 1997; Meszaros & Rees, 1997b; Kluzniak 
& Ruderman, 1998; Drenkhahn & Spruit, 2002; Vlahakis & Konigl, 2003a,b; Lyutikov & Blandford, 2002, 
2003, 2004; Lyutikov, 2004, 2006). These models differ in the initial ratio of magnetic to particle energies 
and in the large scale magnetic field configuration. During the evolution of the flow the energy can be 
transferred between the different components of the flow (bulk, heat and magnetic) and different models 
predict different final flow compositions (in the sense that there are different flnal energy ratios between 
these components). Some models predict that by the time that the prompt emission is generated, most, or 
comparable amount, of the energy is in a baryonic component (e.g., Katz 1997; Vlahakis & Konigl 2003a, b) 
while in other models the magnetic field is the dominant component at all times (e.g., Meszaros & Rees 
1997b; Lyutikov & Blandford 2003). 

Introducing magnetic fields enables models to achieve final Lorcntz factors that are significantly higher 
than the one obtained in pure baryonic models (e.g., ~ 10^—10^ in Meszaros & Rees 1997b) and it also relaxes 
the constraints on the minimal baryonic load, allowing for a pure magneto- leptonic fiow (e.g., Thompson, 
1994). Another example is Vlahakis, Peng & Konigl (2003a) that suggest a scenario in which the interplay 
between the magnetic and baryonic energy making it possible to increase the baryonic load in the form of 
neutrons that decouple from the flow while its Lorentz factor is still moderate. 

Regardless of the acceleration mechanism, the final composition of the flow is the relevant initial condition 
for the prompt emission phase. As we see below it mainly affects the dissipation process that gives rise to 
the prompt emission. 

3.4. The prompt emission 

The prompt emission is generated (or at least re-processed) at a radius in which the optical depth of the 
relativistic flow is low. Any heat that is generated in the flow while it is optically thick is quickly lost to 
adiabatic expansion, and therefore the outflow energy must be dissipated into internal heat near the radius 
in which the emission takes place. Dissipation processes vary according to the composition of the flow. In 
a baryonic or weakly-magnetized plasma the dissipation may result from strong shocks, while in Poynting- 
flux-dominated flows the heat is produced by dissipation of the magnetic field (e.g., via reconnection). 

3.4.1. Internal or external dissipation? 

The dissipation can be done either by processes that are internal to the fiow (e.g., internal shocks) or 
by interaction with an external medium (e.g., external shock). In the case of a baryonic flow an important 
difference between the two dissipation modes is their implications for the central engine. In internal dis- 
sipation models of baryonic outflow the duration of the burst, T, is related to the duration of the engine 
activity. Since the light-crossing time of the engine is shorter by many orders of magnitude than the burst 
duration, a long-lived engine is required. In these models the observed variability time scale, St, reflects the 
variability of the flow, induced either by the engine itself or by interactions at the base of the flow during the 
acceleration. On the other hand, the external shock model (e.g., Rees & Meszaros, 1992; Katz, 1994; Dermer 
& Mitman, 1999), which is the most popular of the external dissipation models of baryonic outflow, allows 
for an explosive source, i.e., the energy can be released on an engine dynamical time scale. If the source is 
explosive, the duration of the burst in this model is determined by the radius in which the interaction takes 
place (T ^ R/cT"^), while the variability time scale may be determined by the size of density fluctuations in 
the external medium. 

In long GRBs the consensus is that if the outflow is baryonic the prompt emission results from internal 
dissipation (see however Dermer & Mitman 2004; for a comprehensive discussion of the topic see Piran 
2005b). The main argument in favor of internal dissipation is the rapid, high amplitude variability observed 
during the prompt emission of long GRBs. Namely, the prompt emission flux varies by 100% over time 
scales, St <C T. The main obstacle to highly variable emission from external shocks is the angular smoothing 
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(Fenimore, Madras & Nayakshin, 1996) and detailed hydrodynamical analysis shows that quasi-spherical 
external shocks cannot produce the observed variability (Sari & Piran, 1997). In quasi-spherical symmetry 
the angular time (tang) is always the longest observed time scale and therefore it is comparable to the 
total duration of the burst, thereby washing out any variability that arises during the emission (for a 
definition of quasi-sphericity and relevant time scales see §3.1). Thus, viable external shock models must 
invoke spatial variability on small angular scales (^ 1/r). The main attempt to do so suggests a clumpy 
external medium (Dermer & Mitman, 1999, 2004). In the context of long GRBs this model encounters 
difficulties with achieving the required variability together with a reasonable efficiency (Sari & Piran, 1997; 
Nakar & Piran, 2002c). furthermore, in this model the width of the relativistic shell. A, sets a lower limit 
on the variability time scale [St > A/c) implying A <C R/T^ . Any variability during the acceleration of the 
flow will result in A ^ R/T^ and therefore this model requires the source of the flow to be extremely, some 
may say unphysically, steady. 

Many SHBs show variability that is similar to the one observed in long bursts (§2.1.2; Nakar & Piran, 
2002b) and therefore the arguments discussed above in favor of internal dissipation if the outflow is baryonic 
are applicable in the case of variable SHBs as well. In the case of SHBs there are additional arguments 
against external shock as the prompt emission dissipation process. First, variable external shocks require 
the ambient medium to have numerous dense clumps. Such an environment may be formed by winds from 
massive stars, the progenitors of long GRBs. However, the proposed progenitor models for SHBs predict 
a rather smooth circum-burst medium. Second, the afterglow is most likely produced by external shock. 
Therefore, if the prompt emission arises from an external shock as well, then the afterglow is expected to 
join smoothly with the prompt emission. However, the observed afterglows show a clear distinction between 
the prompt emission, the following X-ray tail and the late X-ray afterglow (Fig. 4). To produce all these 
phases by external shocks, even if the prompt emission is not variable, requires a contrived external density 
profile. Finally, as discussed in §2.2, there are SHBs with very faint X-ray afterglows. These bursts lack 
bright external shock emission following the prompt emission and therefore, for any reasonable external 
density profile, the prompt emission should arise from a different process. 

If the outflow is Poynting-flux-dominated then the prompt emission energy is generated by dissipation of 
the outflow magnetic energy. This dissipation may be affected by interaction with the external medium (e.g., 
Thompson, 2006). The observed prompt emission variability can be generated in a Poynting-flux-dominated 
flow by emission process that is not isotropic in the rest frame of the expanding outflow (e.g., Lyutikov 
& Blandford, 2003; Thompson, 2006). In such case not all the l/P angular region is observed but only 
small regions within it where the the emission happens to point towards the observer. This way the angular 
smoothing is avoided and high variability may be archived also in external dissipation process (see more 
details in §3.4.3). Note that in the two specific models of the prompt emission in a Poynting-fiux-dominated 
flow that are currently available (Lyutikov & Blandford, 2003; Thompson, 2006), the duration of the burst 
refiects that activity time of the engine. 

To conclude, if the outflow is baryonic then variable SHBs (which are the majority) are most likely 
produced by internal dissipation. In a Poynting-flux-dominated flow the dissipation process can be either 
internal of external. All viable models in both flow types relate the observed duration of the burst to the 
engine activity duration (e.g., Kobayashi, Piran & Sari, 1997; Nakar & Piran, 2002a; Janka et al. , 2006; 
Lyutikov, 2006; Thompson, 2006). SHB that shows a smooth prompt emission light curve may arise from 
an external shock in case that the afterglow joins smoothly with the prompt emission (no such cases were 
observed so far). 

3.4.2. Internal shocks 

In a baryonic or weakly magnetized flow the suggested dissipation process is internal shocks (e.g., Rees & 
Mcszaros, 1994; Narayan, Paczynski & Piran, 1992; Paczynski & Xu, 1994), namely, shocks that arise in a 
variable flow when plasma portions with different velocities collide. In this model each pulse in the light curve 
corresponds to a single episode of collision that takes place at a radius R ~ cMT'^ = 3 • 10^^ cm St -2^2- The 
total duration of the burst corresponds to the total width of the outflow, and is therefore not related to R. At 
these radii the plasma density is low enough so the flow is optically thin to Thomson scattering, as implied 
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by the observations. The low plasma density also implies that the mean free path for a binary collision 
between plasma particles (e.g., protons) is much larger than R and therefore these shocks are collisionless 
and are mediated by collective electromagnetic plasma instabilities (§4.10). Collisionless shocks arc believed 
to efBcicntly accelerate particles to relativistic velocities (§4.10) and the propagation of the accelerated 
relativistic electrons in the plasma magnetic field is thought to be the source (vie e.g., synchrotron process) 
of the observed prompt emission. 

The main advantage of internal shocks is that they are expected to arise naturally in the baryonic outflow 
and that they can easily reproduce the observed variable light curve (e.g., Kobayashi, Piran & Sari, 1997). 
In fact the observed variability reflects the variable activity of the engine and the duration of the burst 
is comparable to, or at most larger by an order of magnitude than (Janka et al. , 2006), the duration of 
engine activity. Additionally, the simplest interpretation of the emission as synchrotron radiation roughly 
predicts the correct observed frequency. The relative velocities between different plasma portions are mildly 
relativistic and therefore internal shocks heat the protons to mildly relativistic temperatures. Assuming that 
the shocks couple the electron and magnetic energies to the proton energy, the electrons have typical Lorentz 
factor 10^ — 10'^, and gyrate in a ^ 10^ G magnetic field (in the plasma rest frame) producing ~ 100 keV 
synchrotron emission (in the observer frame) . 

The main disadvantage of the internal shocks model is its efficiency. Internal shocks are expected to 
radiate only a small fraction of the total flow energy, ^ 1% (e.g., Kobayashi, Piran & Sari, 1997; Daigne & 
Mochkovitch, 1998; Kumar , 1999). Being mildly relativistic, internal shocks convert only a small fraction of 
the bulk kinetic energy into internal energy, and only a fraction of this internal energy is carried by radiating 
electrons. As a result, most of the flow energy is not radiated and remains to be dissipated by an external 
shock, resulting in a bright afterglow. In long GRBs the afterglow is not as bright as expected, suggesting 
that typically as much as 10 — 50% of the relativistic flow energy is emitted in the form of gamma-rays and 
in some cases the gamma-ray efficiency may even be higher (Freedman & Waxman, 2001; Lloyd- Ronning 
& Zhang, 2004; Granot, Konigl & Piran, 2006; Fan & Piran, 2006). As I show in §4.5 the efficiency of the 
prompt emission in SHBs is most likely high as well (see also Bloom et al. , 2006) . Therefore, efficiency is 
potentially a problem of this model in SHBs as well. The same extensions of the internal shocks model that 
were suggested in order to resolve this issue in the context of long GRBs are applicable to SHBs as well, 
e.g., very large amplitude Lorentz factor variability (Beloborodov, 2000; Guetta, Spada & Waxman, 2001) 
and repeated shock crossings (Kobayashi & Sari, 2001). 

The interpretation of the prompt emission as optically thin synchrotron radiation predicts a spectrum that 
is softer than dN/dv cx v~^/^ (Rybicki & Lightman, 1979). This prediction is in conflict with the observed 
spectrum of most SHBs (§2.1.3; Ghirlanda, Ghisellini & Celotti, 2004). This problem, raised by Preece et 
al. (1998) in the context of long GRBs, is more severe in short GRBs where the low energy spectrum is 
harder. Again, different models that were put forward to explain the hard spectrum of long GRBs are also 
applicable here. Some of these models are: random magnetic fields on very short length-scales, leading to 
jitter radiation instead of synchrotron (Medvedev, 2000), effects of synchrotron self absorption or anisotropic 
electron pitching angles (Lloyd & Petrosian, 2000; Granot, Piran & Sari , 2000), Upscatter of optical self- 
absorbed photons by inverse Compton (Panaitescu & Meszaros, 2000) and Photospheric effects (Meszaros 
& Rees, 2000). 

On the high energy end, most SHB spectra are consistent with having an exponential cut-off, while long 
GRBs show a power-law (§2.1.3; Ghirlanda, Ghisellini & Celotti, 2004). If this cut-off is intrinsic, and not 
an observational selection effect, then it most likely reflects the energy distribution of radiating electrons, 
suggesting that in SHBs, unlike in long GRBs, the electrons are not accelerated to very high energies during 
the prompt emission phase. In fact, the observed spectrum in most SHBs can be produced by a mildly 
relativistic shocks that only couple the electrons and the protons, without further accelerating electrons 



^®If the flow is accelerated by radiation pressure alone than the prompt emission light curve directly reflects the energy 
deposition rate and/or baryonic load modulation at the base of the flow (Kobayashi, Piran & Sari, 1997; Nakar & Piran, 
2002a). If the flow is accelerated continuously also at large radii (e.g., by rarefaction waves as discussed in §3.3.1) then it 
reflects the structure of the flow at the end of the acceleration phase (Janka et al. , 2006; Aloy & RezzoUa, 2006). 

Note however that this model predicts larger variance in the observed typical frequency than suggested by BA TSE observa- 
tions. 



32 



to higher (power-law) energies as required in long GRBs. In this case there are no electrons with Lorentz 
factors above 10^ — 10'^ and an exponential cut-off in the spectrum is expected at 100 keV. I speculate 
here that such difference may arise from different magnetization level of the outflow, where coUisionless 
shocks during the prompt emission of SHBs are magnetized while those in long GRBs are not. Numerical 
simulations (Spitkovsky, 2005) suggest that first order Fermi acceleration of particles may not take place 
in magnetized shocks while it has been recently suggested (Lyubarsky, 2006) that in such magnetized 
shocks the electrons can reach equipartition with the protons (with a quasi-thermal energy distribution) 
through the the synchrotron maser instability. The reason that strong magnetic fields at the shock front 
may suppress electron acceleration is that they prevent downstream thermal electrons from crossing the shock 
back into the upstream, thereby suppressing the Fermi process. On the other hand, theoretical considerations 
(Milosavljevic, Nakar & Spitkovsky, 2006) and numerical simulations (Spitkovsky, 2005) suggest that in 
unmagnetized shocks some of the thermalized downstream electrons can cross the shock back into the 
upstream and thus start the Fermi acceleration cycles (see §4.10). 

3.4.3. Prompt emission from a magnetized flow 

If the outflow energy is carried mostly by Poynting flux up to the point where the interaction with the 
external medium starts, the dissipation of the bulk energy into internal energy cannot occur through internal 
shocks^" . In this case the energy source of the radiation is most likely magnetic dissipation processes (e.g., 
reconnection) . Currently there are two models of magnetized flows treating the prompt emission in some 
detail. 

Lyutikov & Blandford (2003) suggest that electromagnetic current-driven instabilities dissipate magnetic 
energy into heat and high energy particles (see also Lyutikov, 2006). The propagation of these high energy 
particles within a strong magnetic field, naturally present in this model, is the source of the observed prompt 
7-ray emission. In this model the dissipation takes place close to the deceleration radius 10^^ — 10"'^^ cm) 
while the outflow Lorentz factor is F 100— 1000. Therefore, the observed variability time scale implies that 
the emission is not spherical (the observed variability, 6t, is much shorter than tang this case). Lyutikov & 
Blandford (2003) suggest that the reconnection of the magnetic field produces, additionally to high energy 
particles, also relativistic turbulence with a Lorentz factor 74 ^ 1, as measured in the rest frame of the 
expanding shell. This bulk relativistic motion of the emitting matter causes the radiation to be beamed into 
an angel of « 1/ (F7t) <C 1/F. This way, each observed piflse is coming from a small angular region (^ 1/r) 
on the expanding shell, allowing St ^ tang = i?/(2cF^). In this model the burst duration is greater or equal 

to tang, SO that St <^ T. 

Thompson (2006) suggests a fiow which is a combination of a strong magnetic field and pair-radiation 
plasma. In this model the radius in which the outfiow becomes optically thin is determined by the pair 
enrichment of the ambient medium through the interaction with the flow radiation (Thompson & Madau, 
2000; Meszaros, Ramirez-Ruiz & Rees, 2001; Beloborodov, 2002; Kumar & Panaitescu, 2004). The magnetic 
dissipation is triggered by the interaction with the pair loaded external medium and is regulated to take 
place around the radius where the outflow becomes optically thin. The prompt emission in this model is 
generated by inverse Compton scattering of the seed photons that are carried in the flow by pairs which 
are accelerated by the reconnecting magnetic field. Highly variable light curve is obtained by emission that 
is directed in the rest frame of the relativistic outfiow. The peak spectral energy of the prompt emission 
(observed at ^ 500 keV) reflects the temperature of these seed photons (Lorentz boosted to the observer 
frame). Therefore, it is determined by the radius at which the fiow starts expanding freely, and by the 
Lorentz factor at this radius. In long GRBs the observed spectrum corresponds to a Lorentz factor of a 
few at a radius of ~ 10^" cm which may be naturally explained as the radius of the progenitor star. It is 
unclear at this point if this model can be applied to SHBs as well, since the recent observations indicate 
that in SHBs the free expansion should also start at a large radius. For example, Thompson, Meszaros & 
Rees (2006) find that in the case of SHB 050709 the peak spectral energy corresponds to a fiow that is still 



" Note that outfiow models that arc dominated by magnetic energy at small radii but convert their energy into particle kinetic 
energy during the acceleration, can use internal shocks to produce the prompt emission. 
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mildly rclativistic and thcrmalizcd at ^ 10^" cm. Currently it is unclear how the thermalization radius can 
be so high in SHBs. 

3.5. High energy Cosmic-rays and neutrinos 

3.5.1. High energy Cosmic-rays 

Cosmic-ray acceleration is typically discussed in the context of long GRBs (for a recent review see Waxman, 
2004). Long GRBs were suggested as the source of extra-galactic ultra-high energy cosmic rays, with observed 
energies between lO^^ - 10^° eV (UHECRs; Waxman, 1995; Milgrom & Usov, 1995; Victri, 1995; Wick, 
Dermer & Atoyan, 2004; Dermer & Atoyan, 2006). This idea is supported by the similar total flux of long 
GRB gamma-rays and UHECR and by the small number of other candidate UHECR sources. Cosmic rays 
are suggested to be accelerated in long GRB internal and reverse shocks, where diffusive shock acceleration 
(a.k.s first order Fermi acceleration; §4.10) may be able to accelerate particles to 10^° eV. External shocks 
are not a promising UHECR acceleration site via diffusive shock acceleration (Gallant & Achterberg, 1999; 
Milosavljevic & Nakar, 2006a) and therefore long GRBs are expected to be UHECR sources only if the 
outflow is baryonic (§3.3), unless particles arc accelerated in the external shock by an alternative mechanism 
(e.g., second order Fermi acceleration; Dermer & Mitman, 1999; Dermer & Humi, 2001; Dermer, 2006). 

The three major factors that limit shock acceleration arc time, conflnemcnt and cooling. Assuming that the 
coherence length of the magnetic fleld in the upstream (unshockcd wind) is larger than R/T, the criterion 
that the acceleration is faster than adiabatic cooling is Rl < R/T where Rl is the Larmor radius of 
accelerated particles in the upstream field and R [F] is the radius [Lorcntz factor] of the rclativistic wind 
(note that the shock itself is mildly rclativistic in the wind rest frame) . This condition also guarantees that 
the accelerated particle is confined and that it has sufficient time to be accelerated. Writing this condition 
using the parameters of rclativistic winds one obtains (Waxman, 1995): 



where Ep^„iax is the maximal energy to which protons can be accelerated and is the wind luminosity, e^j" 
is the fraction of magnetic field energy density in the unshocked wind^^ out of the total wind energy density 
(including rest mass energy) as measured in the wind rest frame. The requirement that the acceleration is 
faster than synchrotron cooling is satisfied for R > lO^'' cm £^^20^2^^ (Waxman, 1995). If SHB winds are 
baryonic then both conditions can be satisfied during the internal shocks and the reverse shock, given that 
100 ^ F < 1000 and that the wind is mildly magnetized. Therefore, SHBs may be able to accelerate 10^° 
eV cosmic rays just like long GRBs. However, The total SHB gamma-ray energy fiux is lower by two orders 
of magnitude than that of long GRBs and the observed UHECRs (Table 1). Therefore, while SHBs may be 
sources of high energy cosmic rays, they are unlikely to be the dominant source of the observed UHECRs. 

3.5.2. High energy neutrinos 

High energy neutrinos (> TeV) are also usually discussed in the context of long GRBs (for reviews 
see Waxman, 2000; Meszaros & Razzaque, 2006). In long GRBs, four sites were suggested as primary 
sources of high energy neutrinos: (i) Acceleration of a neutron rich outfiow (Bahcall & Meszaros, 2000; 
Derishev, Kocharovsky & Kocharovsky, 1999a). (ii) Internal shocks (Waxman & Bahcall, 1997, 1999; Rachcn 
& Meszaros, 1998; Dermer & Atoyan, 2003) (iii) Reverse shocks (Waxman & Bahcall, 2000). (iv) In the course 
of a long GRB jet drilling through the envelope of its massive-star progenitor (Meszaros & Waxman, 2001). 

The last site is irrelevant for SHBs, as the progenitors are not massive stars. All other processes may take 
place in SHBs if their outflow is baryonic. If the outflow is Poyniting-flux-dominated then large neutrino 
flux is unexpected unless neutrinos are efficiently produced in the forward shock (Dermer, 2002; Li, Dai & 
Lu, 2002). This may be possible in the less likely case that ^ 10^^ eV protons are accelerated in the forward 
shock. 

Eg" is not to be confused with tlic magnetic field density in the shocked downstream, sb, which is the site of synchrotron 




(12) 



emission. 
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The energy of neutrinos produced during the acceleration of a neutron rich outflow is ^ 10 GeV and 
therefore these are hard to detect (due to neutrino detector sensitivity at this energy range). Reverse shock 
neutrinos are produced by interaction of ~ 10^" eV cosmic rays with UV photons. The UV flux in SHB 
reverse shocks is significantly fainter compare to long GRBs (§4.3) and therefore this process is inefficient 
in SHBs. Internal shock neutrinos are produced by interaction of > 10^^ eV cosmic rays with the prompt 
gamma-rays. The neutrinos are produced by the decay of pions that result from this interaction. The 
efficiency of this process depends on the efficiency of the cosmic ray acceleration, the observed gamma-ray 
flux and the variability time scale, all of whom may be similar in both phenomena. In long GRBs Waxman 
& Bahcall (1997) estimate that in optimal conditions about 10% of the burst energy is in the form of ~ 10^^ 
eV neutrinos. Assuming a similar efficiencies for long and short GRBs out of the total gamma-ray emission 
(Table 1) implies that the neutrino flux from short GRBs is: 

el'^.Aeu ~ 10^^ eV) - 5 • 10'" GeV/cmVs/sr. (13) 



4. The afterglovif - Theory 



X-ray and optical emission observed hours, days and weeks after some SHBs resemble the afterglow 
observed in long GRBs. Any GRB model that is not 100% efficient in converting energy into gamma-rays 
predicts an interaction of the remaining relativistic ejecta with the external medium, which is expected 
to produce radiation at longer wavelengths. Therefore, various models of the prompt emission, that were 
developed for long GRBs and are applicable to SHBs as well, also predict similar afterglow emission. Indeed, 
the existence of SHB afterglows was predicted long before these were observed (Panaitescu, Kumar & 
Narayan, 2001; Nakar & Piran, 2002b). Differences between long and short GRB afterglows are expected 
mostly due to different total energies, different properties of the ambient medium and, possibly, different 
outflow geometries. I give here a brief description of the standard external shock afterglow model and 
emphasize the interpretation of observed SHB afterglows in this context. I extend and adjust the basic 
model in places where the speciflc SHB parameter space requires it. The reader is referred to Piran (1999, 
2005b) and Meszaros (2006) for further details on the model . 



4.1. Standard model - synchrotron radiation from a spherically symmetric adiahatic blast wave 



The remaining relativistic flow energy, after the end of the prompt emission phase, is transferred to the 
external medium by driving a shock-wave into it - the external shock. At early times the energy content in 
the ejecta and the freshly-shocked external medium is comparable. At this time the ejecta composition has 
a strong effect on the dynamics and the emission. I discuss this early, but more complicated, phase later 
(§4.3). After all the energy is dissipated into the external medium the only memory of the initial conditions 
is the total amount of energy and its angular distribution (i.e. the initial geometry of the ffow). This phase, 
generally referred to as the late afterglow, is discussed here. 

Consider an adiabatic and spherical ^'^ blast wave with energy E that propagates with Lorentz factor F 
into an external medium composed of proton-electron plasma with a constant density n. The random motion 
Lorentz factor of the shocked plasma in the blast wave rest frame is « F and therefore the total blast wave 
energy, as measured in the observer frame, is « M{R)T^c^, where M{R) is the mass accumulated by the 
blast wave up to a radius R. Energy conservation dictates (e.g.. Sari, Piran & Narayan, 1998): 

,.,o[0''n,r«40'ij^y"\ (14) 



^^SHB afterglow models that I do not cover here are the "cannonball" model (for review see Dar & de Riijula, 2004) and a 
cylindrical external shock model (Wang, Cheng & Tam., 2005) 

The spherical approximation can be used here as long as the blast wave properties do not vary over angles that are smaller 
than 1/r (quasi-sphericity of relativistic flows). E in this model corresponds to the isotropic equivalent kinetic energy of the 
outflow. 
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where 7 is the Lorentz factor of the freshly-shocked external medium (the Lorentz factor of the shock itself 
is r = and td is the observer time in days {t « i?/47^; Waxman, 1997a; Sari, 1997). Conservation of 

mass, momentum and energy fluxes across the shock imply that the proper number and energy density of 
the freshly shocked plasma (in the plasma rest frame) are 4771 and 47^nmpC^ respectively. The profiles of 
the hydrodynamic quantities (pressure, density and velocity) behind the shock are self similar (Blandford 
& McKee, 1976) and the entire shock energy is concentrated in a very thin shell with a width A sa i?/10r^. 
The simplest model assumes that the shock couples the electrons to the proton energy and accelerates all 
the electrons to a power-law dN/dE oc starting at a minimal Lorentz factor 7e,m = ')£e{'mp/'nie)(p — 
2)/(p — 1), where p > 2 and Se is the fraction of the internal energy that is carried by the electrons. After 
crossing the shock the electron and the proton temperatures are not coupled any more, and the electrons 
cool down by radiation and pdV work. This model assumes further that the magnetic field in the shocked 
plasma carries a constant fraction £5 of the internal energy, where the magnetic field of freshly shocked 
plasma \s B = (32nmpC^£5)^/^7. Note that by definition Sb < 1. 

This model has five free parameters: E, n, Se, £b and p. These parameters can, in principle, take almost 
any value, but observation of long GRBs together with the properties of the host galaxies of SHBs suggests a 
"natural" range for these parameters. Assuming that the processes that produce the prompt emission in long 
and short GRBs are similar, we expect similar efficiency, implying that the 7-ray isotropic equivalent energy 
output is a reasonable estimator of E (Freedman & Waxman, 2001; Granot, Konigl & Piran, 2006; Fan & 
Piran, 2006), and therefore E ^ 10^" erg (see Table 2). Observed SHB host galaxies, and the location of the 
bursts within them, indicate that the progenitors are associated with an old stellar population. Moreover, 
various progenitor models suggest that the peculiar velocity of the progenitors can be high enough to escape 
from the gravitational potentials of their hosts. Therefore, possible environments of SHBs are the inter-stellar 
medium (ISM), intra-cluster medium (ICM) or the inter-galactic medium (IGM). In all of these cases the gas 
density would be rather constant on the scales relevant for SHBs, and the corresponding density values are 
10-2 - 1 cm-3 (ISM), 10-3 - 10-2 cm-3 (ICM) and ~ 10-^ cm-^ (IGM). Finally, the values of £e, £b and 
p are determined by coUisionless shock physics, which take place on microphysical scales and are therefore 
unlikely to depend on global properties such as the total energy, the density profile, etc. It is therefore 
reasonable to use the values observed in long GRBs (see also §4.10): £e = 0.05 — 0.5, Eb = 10^'^ — 0.1 and 
2 < P < 3 (e.g., Panaitescu & Kumar, 2001; Yost et al. , 2003). 

Within this parameter range the radio to X-ray afterglow is dominated by synchrotron radiation. The 
synchrotron spectrum has three characteristic frequencies (e.g., Meszaros & Rees, 1997a; Waxman, 1997b; 
Sari, Piran & Narayan, 1998; Granot, Piran & Sari, 1999b): Vm - the typical synchrotron frequency (oc ^B^l) 
that corresponds to the synchrotron frequency of an electron with Lorentz factor 7e,m; i^c - the synchrotron 
frequency that corresponds to 7e,c, which is the Lorentz factor above which radiative cooling is significant; 
and Va - the synchrotron self-absorption frequency. The observed spectrum is usually four-segment broken 
power-law with breaks at the characteristic frequencies, where the power-law indices depend on the order 
of these frequencies. The value of Va and of the peak value of F^, depend on the order of the characteristic 
frequencies as well. In the parameter range relevant for SHBs the typical frequency order would be Va < 
Vm < i^c during the relativistic phase (i.e., cooling is slow). In this case the characteristic frequencies are^^ ; 



^''Eq. 16 is generalized to include self-synchrotron Compton cooling (see §4.2) by taking Uc — > fcCl -I- ^) ^ where Y is defined 
in Eq. 22. As I show, in SHBs y < 1. 




(16) 



Fi^^m ~ 5 ^Jy (1 + 2;)£B,i2^50'^-'|rfi,,28' 



36 



where -Fl/^m is the flux at i/m- The observed flux is a broken power-law with the following spectral and 
temporal indices: 



1/ < J/a 

l^~^t~i^P~^^ V,n <V <Vc 

3p-2 



(17) 



V 2 i 



Vc<V 



A detailed discussion of the light curve in cases of other orders of the break frequencies can be found in 
Granot & Sari (2002) while the case of 1 < p < 2 is discussed in Dai & Cheng (2001). 

Once the blast wave accumulates enough mass it decelerates to non-relativistic velocities and makes 
a transition from the relativistic self-similar Blandford-Mackee solution (Blandford & McKee, 1976) to 
the Newtonian self-similar Sedov-Von Neumann- Taylor solution (Sedov, 1946; von Neumann, 1947; Taylor, 
1950). The transition between the two solutions is gradual and takes about one decade in time around: 



tNR ~ 0.5 yr 



£^50 

n-2 



1/3 



(1 + ^) 



(18) 



At tNR the blast wave velocity is w oc i "^/^ and the flux at the radio frequencies, vr, which is typically 
the only observed wavelengths at this time, behaves as: (Waxman, Kulkarni & Frail, 1998; Frail, Waxman 
& Kulkarni, 2000): 



- + - 

^ 5 ■ 



I^a, l^m < VR < Vc 



t » tNR 



(19) 



4.1.1. Comparison to the observations 

The handful of observed SHB afterglows were analyzed within the framework of the standard afterglow 
model in several papers (Hjorth et al. , 2005b; Fox et al. , 2005; Berger et al. , 2005; Lee, Ramirez-Ruiz & 
Granot, 2005; Bloom et al. , 2006; Panaitescu, 2006; Campana et al. , 2006; Burrows et al. , 2006; Grupe 
et al. , 2006; Soderberg et al. , 2006a). The simple model gives a reasonable explanation of the general 
properties of most of the observed afterglows. Moreover, most of the observed properties were predicted 
by the standard model (e.g., Panaitescu, Kumar & Narayan, 2001) and while the current observations are 
too sparse to confirm its validity, they are certainly encouraging. For example, this model predicts that the 
emitting electrons are accelerated in unmagnetized relativistic collisionless shocks to a power-law distribution 
with index p. Long GRB afterglow observations as well as theory of particle acceleration in relativistic shocks 
suggests 2 < p < 3 (§4.10). This prediction determines, with little freedom, the optical and X-ray spectral 
and temporal evolution (Eq. 17). Most SHB afterglows conform with the evolution of the afterglow as 
predicted by this range of p values. A few examples are SHBs 050709 {p « 2.5; Fox et al. , 2005; Panaitescu, 
2006)), 050724 (p=2-3; Berger et al. , 2005; Panaitescu, 2006) and 051221 {p « 2.2; Soderberg et al. , 2006a; 
Burrows et al. , 2006) where broad-band modeling (radio to X-ray) gives a reasonable description of the 
observed afterglows. 

Unfortunately, the current quality of the data is not sufficient to break all the degeneracies between the 
other four parameters of the model {E, n, Se and es), and to determine their values. Nevertheless, in most 
cases where the simple standard model fits the data, the default parameters that were predicted based 
on the observations of long GRBs (i.e., E ^ E^jso, n < 1 cm~^ and Se and Eb not much smaller than 
equipartition), give a satisfactory description of the data (e.g., Panaitescu, 2006). 

On the other hand, some of the observed details cannot be explained by the simple, spherically symmetric, 
self-similar model. The most prominent features that require deviation from the basic model are the early 
X-ray tail (see §2.2.2), the late flares observed in the X-ray afterglows of several bursts (Fig. 4) and the late 
break in the light curve of SHB 051221. Extensions of the basic model that can explain at least some of 
these features are discussed below. 
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4.2. Synchrotron self-Compton 



The same electrons that produce synchrotron radiation also upscatter synchrotron photons by the inverse 
Compton process (synchrotron sclf-Compton - SSC; e.g., Meszaros, Rees & Papathanassiou 1994; Waxman 
1997c; Wei & Lu 1998, 2000; Panaitescu & Kumar 2000; Dermer, Chiang & Mitman 2000; Sari & Esin 2001; 
Zhang & Meszaros 2001b). SSC photons themselves are upscattered again by the same relativistic electrons 
to even higher energies producing a second generation of inverse Compton photons. The process continues 
to higher and higher generations until the photon energies, as seen in the electrons' rest frame, are above 
the Klein-Nishina energy (meC^), where the scattering cross-section drops rapidly (roughly linearly) with 
photon energy. For typical long GRB parameters the first generation of the SSC is below the Klein-Nishina 
limit and is usually calculated without considering the effect of the Klein-Nishina cross-section. This first 
SSC generation is predicted to produce bright GeV emission that would be detected by sensitive GeV future 
observatories (e.g., GLAST). The second SSC generation for long GRB afterglows is well above the Klein- 
Nishina limit and is completely suppressed. As I show below for SHBs the situation is different as the first 
generation SSC component is already affected by the Klein-Nishina cross-section, though it is not entirely 
suppressed. I therefore repeat the analysis presented in Sari & Esin (2001) augmented by a simple estimate 
of the effect of the Klein-Nishina cross-section (c.f., Li & Waxman, 2006). The result is that in all but the 
brightest SHBs the SSC emission is expected to be significantly less energetic than the synchrotron emission. 

Photons that are scattered by an electron with a Lorentz factor 7e are below the Klein-Nishina limit as 
long as the observed frequency satisfies: 

r 

hi' ^ huKNile) = rUeC^ — (20) 

7e 

where h is the Planck constant and F is the plasma Lorentz factor (7e is measured in the plasma rest-frame) . 
In the slow cooling regime (i.e., v„i < &s expected for SHBs) most of the synchrotron energy is emitted 
around Vc and most of the SSC energy is emitted by ~ 7c electrons that upscatter ~ Vc photons. Therefore, 
the Klein-Nishina limit can be neglected only if Vc/vkn{1c) ^ 1: 

« UfE^,'r.-Te-^"_l{l + Y)-\ (21) 

Here I included the effect of SSC cooling on Vc using the Y parameter (defined in Eq. 22 and found to be 
small). Eq. 21 shows clearly that in SHBs the Klein-Nishina limit cannot be neglected and that the first 
SSC branch in this case is strongly affected by the Klein-Nishina cross-section^^ . 

The total energy in the first SSC generation can be estimated as follows. In the relevant parameter space 
Vm < vkn{.1c) < Vc where oc implying that for 2 < p < 3 most of the upscattered energy is 

emitted by ^ 7c electrons that upscatter ^ vkn photons. Therefore, the ratio of the SSC luminosity to the 
synchrotron luminosity is (c.f. Eq. 3.1 of Sari & Esin 2001): 

" Lsyn ^ Ub " (1 + Y)Ub " (1 + Y)eB ' ^ ^ 

where Usyn, Ub and Ue arc the energy density of the synchrotron emitted photons, the magnetic field, 
and the relativistic electrons respectively (measured in the plasma rest frame), rjrad is the fraction of the 
electrons' energy that is radiated: 



Vkn is the fraction of the radiated energy which is in frequencies smaller than VKNilc)'- 

3-p 

Vkn = I I ; Vm < vknVIc) < Vc- (24) 



In long GRBs the isotropic equivalent blast wave energy is larger by 2 — 3 orders of magnitude and the external medium 
density is usually higher. Therefore, the Klein-Nishina cross-section can usually be neglected for the first SSC generation. 
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Expressing Eqs. 23 & 24 using the five free parameters of the model and plugging them into Eq. 22 yields: 

Y = O.Ole^-'^^-'-Phl;_\e§:Xl>'nft~'^ (25) 

This equation is valid for 2 < p < 3 and as long as Vm < ^kn{1c) < J^c, as expected in SHBs. Note that 
in this range of p values the dependence on and t is very weak while by definition Eg < 1- Therefore, 
y <C 1 in all SHBs but the most energetic ones that happens to explode in relatively dense environments. 
The implications of such low Y values is that SSC cooling is inefficient and can be ignored in Eq. 16 and 
that in general SSC GeV emission from most SHB late afterglows will be dim and is unlikely to be detected 
by planed GeV observatories. 

4.3. Early afterglow and the reverse shock 

According to the internal-external dissipation model, between the prompt emission and the late afterglow 
there is an intermediate phase, during which the energy in the ejecta is dissipated into the external medium. 
The composition of the ejecta plays an important role in the physical evolution during this phase. In long 
GRBs a baryonic flow may show a clear observable signature (the reverse shock optical flash and radio flare) 
which is not expected in a Poynting flux dominated flow. As I show below for SHBs, even if the flow is 
baryonic and there is a strong reverse shock, bright early optical emission is not expected. 

Interaction of baryonic ejecta with the external medium generates two shocks. One is propagating into the 
external medium (forward shock) and the other is propagating back into the ejecta (reverse shock). Between 
these two shocks there are two shocked regions that are separated by a contact discontinuity. Similar shock 
structure can be found in many astrophysical environments such as supernova remnants and stellar winds. 
As I show below, in SHBs the reverse shock is mildly relativistic. This implies that during a single crossing 
the shock dissipates most of the bulk motion energy of the outflow into internal energy. After the shock 
crosses the outflow once and a rarefaction wave is reflected, only the forward shock remains, and is the source 
of the late afterglow. The reverse shock is short-lived and might produce an observable signature which may 
be distinguishable from the late afterglow emission. If, however, the outflow is dominated by Poynting flux, 
then the large magnetic pressure prevents the formation of a strong reverse shock. The exact prediction for 
magnetized flow depends on the magnetization level of the plasma (e.g., Lyutikov, 2006; Zhang & Kobayashi, 
2005). 

Observed signatures of reverse shocks were explored extensively in the context of long GRBs (e.g., Meszaros 
& Rees, 1997a; Sari & Piran, 1999; Kobayashi & Sari, 2000; Sari & Meszaros, 2000; Fan et al. , 2002; 
Soderberg & Ramirez-Ruiz, 2002; Zhang, Kobayashi & Meszaros, 2003; Nakar & Piran, 2004; Beloborodov, 
2005; McMahon, Kumar & Piran, 2006) and the same results are applicable, with some minor changes, 
to SHBs (Fan et al. , 2005). I summarize the predictions for reverse shock emission, making the required 
adjustments to SHBs. 

Consider a homogenous cold baryonic shell expanding relativistically into homogenous cold ambient 
medium. The problem is well deflned by the shell (isotropic equivalent) energy E, width A in the lab 
frame, initial Lorentz factor Pg and the ambient density n. As the ejecta plow through the external medium, 
a forward shock and a reverse shock are produced. The strength of the reverse shock is determined by the 
dimensionless parameter (Sari & Piran, 1995): 

e^(VA)i/2r-^/^ (26) 

where I = {3E/{ATrnmpC 

2))i/3 is the Sedov length. If ^ < 1 the reverse shock is relativistic and the shell 
decelerates, while dissipating most of its bulk motion energy, within a single shock crossing of the shell at 
w /■^/*A^/"'. In this case Ra ~ Rdec, where Rdec is the deceleration radius, namely the radius in which 
all the ejecta decelerates by a factor of two or more. For ^ ^ 1 the RS is Newtonian and many crossings are 
required to decelerate and shock dissipate a signiflcant fraction of the ejecta energy. In this case Rdec is the 
radius in which the forward shock collects an amount of mass that is equivalent to ~ E/{(?V'^). 

The Lorentz factor of the outflow is expected to vary by at least a factor of order unity, causing the shell 
to expand once it gets out to a radius Rgpread ~ ^oTq where Aq ~ cT is initial width of the shell (T is the 
burst duration): 
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A«max(Ao,i?/r2), (27) 
Correspondingly, ^ is constant up to Rspread, and it decreases at larger radii. The value of = C(^o) 
determines the evolution of the reverse shock. If < 1 then Rdec < Rspread and the relativistic reverse 
shock crosses the shell before it starts spreading. If ^ 1 the reverse shock is initially Newtonian and the 
shell begins to spread during the first crossing of the reverse shock. As the shell spreads the reverse shock 
grows stronger and at the radius where the reverse shock finally crosses the shell, ^ « 1 and the reverse 
shock is mildly relativistic (its Lorentz factor is « 1.25). 
For typical SHB parameters, 

^„ ^ 43 (^1^) r-^/^T-^/^ » 1. (28) 

Therefore, in most (and maybe all) SHBs the reverse shock is mildly relativistic 1) and does not depend 
on the initial Lorentz factor, the total energy or the burst duration. The deceleration radius in this case is: 

i?,ec« 5- 10^6 cm (29) 

Corresponding to an observer time: 

Uec « ^ (1 + z) « 90 s (1 + z) (1^) T-T. (30) 

At tdec the reverse shock dies away and its optically thin synchrotron emission (at ly > i^a) peaks and then 
starts fading rapidly (e.g.. Sari, 1997). tdec also marks the onset of the "late" afterglow - a self similar 
adiabatic solution - as the forward shock turns into a single blast wave (and indeed Eqs. 30 and 9 are 
similar). Note that in SHBs tdec does not depend on T and that tdec ^ T. This is different than long GRBs 
where Ag is significantly larger and Co ^ Ij in which case tdec ~ T. Thus, in SHBs a gap is expected between 
the prompt emission and the afterglow while in long GRBs the two are expected to overlap. 

The emission from the reverse shock is calculated as in the standard afterglow model, taking into account 
that the reverse shock is mildly relativistic while the forward shock is ultra-relativistic. The strength of the 
shock defines the typical energy per particle in the shocked region. This in turn sets the typical frequency 
of the emitted synchrotron radiation. Thus, while the forward shock is radiating in the X-ray band, the 
reverse shock is radiating at low frequencies. The reverse shock signature in long GRBs is an optical flash 
that peaks soon after, or during, the prompt GRB emission and then decays as ^ t~^ (e-g-, Sari & Piran, 
1999). R ^ 18 mag is expected from the reverse shock of a typical long GRB and in extreme cases it can 
be brighter by several orders of magnitude (Nakar & Piran, 2004). An accompanying radio flare is expected 
to peak hours to days later (e.g., Kulkarni et al. , 1999; Nakar & Piran, 2004). In SHBs, where the energy 
and typical density are lower and the reverse shock is not expected to be very strong, the predictions are 
different. For SHB "canonical" parameters (i?5o = 1 , n = 0.01 cm~'^. Eg = 0.1, = 0.01) most of the 
reverse shock emission is radiated around ~ 100 GHz and the early optical emission is faint, ~ 2 /xJy {R ~ 23 
mag) at z = 0.2 . In this case the optical emission from the forward shock is expected to be brighter than 
the reverse shock emission. Actually, the reverse shock emission from SHBs is more likely to be detected by 
early radio observations. In the above case a peak flux of ~ 0.1 mJy is expected in 8 GHz at t « 20 min. 

To conclude, early reverse shock optical flash is not expected from most SHBs (assuming canonical pa- 
rameters), even if the ejecta is baryonic. A reverse shock signature might however be detected by very early 
deep radio observations. 

4.4. "Naked" afterglow 

When the relativistic outflow that produces the prompt emission is quasi-spherical, an afterglow (late soft 
emission) is expected even if no emission from the external shock is observable (e.g., due to low external 

^® If the ejecta is mildly magnetized so the reverse shock is still mildly relativistic but eg ~ 1/3, then the optical emission from 
the reverse shock is brighter {Fu^opt is roughly linear with eg at this parameter range). Additional increase in the reverse shock 
optical luminosity is expected if only a small fraction of the electrons is accelerated, thereby significantly increasing !/„%■ 
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density). This is the high latitude emission of the prompt gamma-rays source (§3.1.3), that decays at late 
times (t ^ tios and t ^ tang] §3.1.1) as oc iy~^t~^^^^\ This type of afterglow ("naked"; Kumar & 
Panaitescu, 2000; Dermer. 2004; Page et al. , 2006) is suggested to be responsible for the first phase of long 
GRB X-ray afterglows'"^ (Nousek et al. , 2006; Zhang et al. , 2006). 

SHBs 051210 (La Parola et al. , 2006) and 050813 (Roming et al. , 2006) show X-ray afterglows that 
may be high latitude emission. The afterglow of 051210 also shows a flare at t ^ 100 s. Assuming that no 
external shock emission exists, this flare is probably a result of the engine activity, in which case tios ^ 100 
s. La Parola et al. (2006) find that assuming this value of tios, the late time decay is consistent with high 
latitude emission. In order for external shock emission to be fainter than the high latitude emission, they 
estimate that the external density is lower than n < 4 • 10^"^ cni^''. 

4.5. X-ray dark afterglows and ^-ray efficiency 

Figure 6 presents the distribution of the dimensionless ratio fx^ = Fxt/ S^, where Fx is the X-ray (0.2— 10 
keV) energy flux at time t and S~f is the prompt emission gamma-ray fluence (15 — 150 keV), of Swift long 
and short GRBs at one day (§2.2.1). Within the framework of the standard afterglow model (§4.1) and as 
long as the blast wave is quasi-spherical: 



represents the 7-ray efficiency of the prompt emission. Ek (equivalent to E in 4.1) is the kinetic energy of 
the blast wave and E^ is the energy emitted in 7-rays (all energies in this section are isotropic equivalent). 
The exact power of the parameters in Eq. 31 (e.g., £e) depends weakly on p (assuming here 2.1 < p < 3). 
For simplicity I use approximate power values. I also neglect weak dependence (power-law indices below 
1/4) on parameters, since these cannot affect the result significantly (the lack of dependence on n is exact). 
Following §4.2 the SSC cooling is neglected as well. 

Assuming that microphysics of coUisionless shocks docs not vary significantly between bursts (either long 
or short) « 0.1 and £b ~ 0.01 are adopted. Under this assumption, for long GRBs at late time (^ 1 d; 
but before the jet-break), one expects < Vx in which case fx-y is almost a direct measure of the 7-ray 
efficiency, k. As evident from Fig. 6, for long GRBs fx^i^d) = 0.1 — 10~^ implying i^lgrb ~ 0.01 — 1. 
This result is well known (Freedman & Waxman, 2001; Lloyd-Ronning & Zhang, 2004; Granot, Konigl & 
Piran, 2006; Fan & Piran, 2006). Figure 6 also shows that the values of fxf for SHBs with observed X-ray 
afterglows are comparable to those of long GRBs, ^ 0.01. 

For some SHBs the circum-burst density can be low, in which case it is not clear wether Vc is above or 
below the X-ray band at 1 day. If n > 0.01 cm"'^ then lyx ^ i^c and kshb ~ 0.1. If the density is significantly 
smaller then Vx > i^c and k decreases as n^^^. Since SHBs with observed afterglows are typically located 
within their host galaxy light, most likely n 3> 10"'^ cm~^, so k w 0.01 — 0.1. We can conclude that at least 
some SHBs (those with observed X-ray afterglow) the gamma-ray efficiency is most likely similar to that 
of long GRBs (see Bloom et al. , 2006; Lee, Ramirez-Ruiz & Granot, 2005, for a specific exploration of the 
efficincy of SHB 050509B). 

Early X-ray afterglow from long GRBs is always detected. In contrast, there are several SHBs with tight 
upper limits on any early (< 100s) X-ray emission. The values of fx-y < 5 • 10~^ for these bursts are 
exceptionally low (Fig. 6). Making the plausible assumptions that the gamma-ray efficiency of these bursts 
is typical (k ^ 0.1) as are the initial Lorentz factor and the microphysical parameters, these values of fx-y 



■^^Note that identification of high latitude emission can be tricky. At early times it requires an identification of ij^j, of the 
specific shell that dominates the late naked afterglow emission. As evident from Eq. 5, a significant underestimate of ti^s leads 
to an overestimate of the decay index at early time. 




(31) 



where, 




(32) 
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indicate that these events occurred in extremely low density environments, n < 10~^ cm~'^, typical for the 
inter galactic medium. As I will discuss below this result suggests a long-lived progenitor with high natal 
velocity (§5.2.3). An alternative would be to relax any of the above assumptions. For example, assuming 
an inter-stellar density (n > 0.01 cm"-^), the low fx^ value can be explained by ultra-efRcient gamma-rays 
production (k > 100), by unusually low electron and magnetic field energies {el^'^es < 10~^) or by low 
initial Lorentz factor Tq < 20. The latter case can explain the faint early afterglow because if Tq is low the 
deceleration time (and therefore the afterglow onset) is at t » 100 s. 

4.6. Angular structure of the outflow 

The angular structure of the relativistic outflow is of great physical interest. As discussed in §3.1.2, during 
the prompt emission phase the observer is sensitive only to a small patch of the flow (a solid angle of '-^ 1/7^ 
sr), and is ignorant to any of the flux emitted outside of this patch. This flux determines the total energy 
emitted by the burst as well as the range of the viewing angles over which a given burst can be detected. 
Therefore the angular structure is a crucial ingredient in determining the total burst energy and in extracting 
the total SHB rate out of the observations. Furthermore, the structure of the jet is most likely determined 
by the engine, and different progenitor and central engine models predict different angular structures. 

Currently, the small number of SHB afterglows prevents the determination of the exact angular structure 
of the outflow (even in long GRBs it is still an open question). Therefore, I discuss here only the simplest jet 
structure: a uniform energy over a half opening angle 9j <^ 1 and no energy outside of this angle (a.k.a 'top 
hat') ^® . This simple model turns out to be consistent with the few available observations (excluding the 
observed flares discussed in §4.7) and it is also supported by numerical simulations (Janka et al. , 2006). A 
comprehensive review of jets in GRBs can be found in Granot (2006), where the dynamics and observational 
signatures of relativistic blast waves with different angular structures are discussed. 

Two key processes govern the observational signatures of a decelerating collimated relativistic blast wave. 
The first is relativistic beaming of the radiation from an emitting element with a Lorentz factor F to an 
angle of 1/F. This implies that the radiation from a jet with half opening angle 9j is confined to an angle 
9j -\- 1/F. The second process is hydrodynamic lateral spreading of the jet, which is limited by the angular 
size of causally connected regions (see §3.1.2), w 1/F. Therefore, the half opening angle of a jet with an 
initial half opening angle Oj ^ is: Oj < Oj^, + 1/F. If the observer viewing angle (with respect to the axis of 
the jet), ^obs, is smaller than Oj^, the two effects combine to produce a light curve that is initially similar 
to that of the spherical blast wave discussed in §4.1. Only when the blast wave decelerates and F becomes 
comparable to l/%.o the jetted nature of the blast wave "reveals" itself to the observer in the form of a light 
curve break (Rhoads, 1997, 1999; Sari, Piran & Halpern, 1999). 

This jet-break in the light curve results from a combination of the two processes discussed above. When 
the radiation beaming angle becomes wider than the jet size, the "missing" emission manifests itself as a 
faster decay of the light curve. Simultaneously, when the jet center become causally connected with the jet 
edges, the jet spreads sideways and the energy per solid angle decreases resulting in a faster deceleration 
with the radius. Therefore, at: 



the light curve 'breaks' roughly simultaneously in all observed bands. The asymptotic decay at tj t t^n 
(defined in Eq. 18) is^^ (Sari, Piran & Halpern, 1999): 



° Observational signatures of the 'top hat' jet are very similar to those of a structured jet where the energy per solid angle is 
oc if the jet opening angle is replaced by the viewing angle (Lipunov, Postnov & Prokhorov, 2001; Rossi, Lazzati & Rees, 
2002; Zhang & Meszaros, 2002). Therefore, all the discussion here is valid for this case as well (replacing 9j with the viewing 
angle). Note that such a structured jet results in a larger beaming-corrected energy and lower beaming-corrected rate than top 
hat jet. 

^^The hydrodynamics of lateral spreading is a major open question. Numerical simulations (Granot et al. , 2001; Cannizzo, 
Gchrels & Vishniac, 2004) and approximate analytical treatment (Kumar & Granot, 2003) suggest that the spreading velocity 
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Comparison with the temporal decay at t ^ t;, (Eq. 17) shows that in the optical and the X-ray (where 
V > I'm) a roughly similar break is expected from a ~ 1 to a ~ 2. Therefore, when an achromatic break, 
with roughly the predicted power-law indices, is observed, it is usually interpreted as an indication that the 
outflow was initially coUimated to a half opening angle: 



To date, three SHBs have known redshift and late afterglow light curves (Fig. 4). The limits on the opening 
angle of these bursts are: 

SHB 050709: In this burst there are optical and X-ray detections at several different epoches (Fox et al. , 
2005; Hjorth et al. , 2005b). Fox et al. (2005) suggest an optical break based on three HST observations, 
but Watson et al. (2006) have shown that when additional optical data are considered, the optical data 
set cannot be fitted by a single or a broken power-law. Additionally, the last X-ray observation (t = 16 d) 
shows evidence of an intense flare. Given the small number of observations and the apparent variability, it 
is impossible to reliably determine an underlying power-law behavior that is consistent in all wavelengths. 
Therefore, unfortunately, the opening angle of this burst cannot be constrained. 

SHB 050724'- Observations include a detailed X-ray light curve (Campana et al. , 2006; Grupe et al. , 
2006) but only sparse detections in the IR, optical and radio (Berger et al. , 2005). At early times the X-rays 
are dominated by the X-ray tail (§2.2.2), and following its decay two flares are observed, one at ~ 1000 s 
and the other around « 10'* — 10^ s. Between the two flares the power-law decays as ^ t~^ and late Chandra 
observations at t w 3 weeks show that after the second flare the afterglow resumes its shallow t~^ decay. 
These observations suggest that tj > 20 d as there is no evidence for a jet-break in the data, implying 
(Grupe et al. , 2006): 6j,o > 0.4 rad (£'fc.iso,5o/"--2)~"'^^^- Possible caveats are that a jet signature may have 
been obscured by the flaring activity or, if the density is indeed as high as suggested by Panaitescu (2006), > 
0.1 cm~'^, the Newtonian transition may already affect the light curve after 3 weeks, preventing the detection 
of a jet. If instead one adopts the beginning of the second flare (~ 10"* s) as a lower limit on the jet-break 
time, the opening angle is 9j^o ^ 0.06 rad {Ek^iso,50 /n^2)~^^^ ■ Assuming 'top hat' jet structure these limits 
imply a total (beaming corrected) energy oi E^^tot = E^^.soOj/^ > 0.7[0.02]10« erg {Ek,^so.,5o/n-2)-^'^ (the 
bracketed value is for the more conservative limit on Oj^). 

SHB 051221: Observations include a detailed X-ray light curve (Burrows et al. , 2006) and continuous 
optical monitoring (Soderberg et al. , 2006a). This afterglow shows a continuous X-ray decay as between 
100 s and 4 d with a short episode of flattening around i = lhr. Atiw5d the X-ray decay becomes 
much steeper (roughly as i~^). The optical light curve decays as w t~°-^ until 4 d where it falls below the 
detection limit. The 3(t optical upper limits after four days require a simultaneous optical break. Here the 
interpretation of the break at 4 d as a jet signature seems quite suggestive given its achromatical nature and 
the regular decay before and after the break at the predicted decay rates. The corresponding jet opening 
angle is O^^q « 0.16 rad {Ek,iso,bi/n-2)~^'^ implying a total energy E-f^tot « 3 • 10''^ erg (£^fc,iso,5i/"-2)"^/''- 

Since there is only one "clean" case (SHB 051221) for which a reasonable constraint on the beaming can 
be obtained, the opening angle distribution of SHB jets is not well constrained. Current observations suggest 
that the average beaming factor (/f, = 1 — cos^j) satisfies 1 ^ (fb^) < lOOi but at this point no secure 
conclusion can be drawn. The fine details of the outflow angular structure remain obviously unconstrained. 



in the co-moving fluid frame is significantly smaller than c. The details of the lateral spreading determine the relative importance 
of each one of the two processes that cause the observed break. In any case, observational signatures do not depend strongly 
on the spreading details. Eq. 34 is derived assuming maximal lateral spreading (comparable to the light speed), but is a fairly 
good approximation for spreading at any other velocity (including no spreading at all). 




(35) 
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4.7. Afterglow variability 

Afterglow model ingredients discussed so far are regular and as such they predict smooth and regular 
light curves (smoothly broken power-law). These models are simplified and the true physical environment 
is expected to be more complicated with irregularities giving rise to light curve variability. Thus, afterglow 
variability, or the lack of it, can probe irregularities in the external shock and late engine activity. 

The external shock cannot produce large amplitude, rapid flux variability (SF/F > 1 over St/t 1 where 
F is the afterglow flux and t is the time since the burst). The reason is that during the external shock 
tang ^ ^Otios ^ lOtA (the width of the emitting shell following the external shock is A « i?/(10r^); §3.1.1). 
Emission from a radius R is observed between t = tios and t « tang + tios ~ tang and therefore spherical 
variations are smoothed over tang, implying St ^ tang ^ t. Even variations in a well localized angular region 
will be smoothed over t^ implying St > t/lO. For an external shock to produce high amplitude variability 
on shorter time scales a very small region (relative to A and R/T) should increase its luminosity so that it 
becomes much brighter than the remainder of the observed shell, and then quickly decay without affecting 
the emissivity of its vicinity. Such behavior is not expected from hydrodynamical processes. Producing 
significant variability even over time scales as short as t^ requires a unique configuration. On the other 
hand, if the engine is active at late times, it can easily produce variability with arbitrary amplitude on time 
scales as short as its dynamical time (< 1 ms), just as we observe during the prompt emission. Therefore, 
whenever a high amplitude, rapid variability is observed, it is likely a result of late engine activity. 

High latitude emission (§3.1.3) from a quasi-spherical shell sets another interesting constraint on the 
variability induced by an external shock (Nakar & Piran, 2003a). A spherical external shock cannot produce 
light curves that decay faster than Fi, oc v^^t^^'^^^^\ Rapidly decaying light curves must result either from 
late engine activity or from an aspherical blast wave on angular scales < l/P. 

Alternative explanations for afterglow temporal variability are: 
(i) Variable external density (e.g., Wang & Loeb, 2000; Lazzati et al. , 2002; Ramirez-Ruiz et al. , 2001; Dai & 
Lu, 2002; Nakar, Piran & Granot , 2003; Nakar & Piran, 2003a; Ramirez-Ruiz et al. , 2005; loka, Kobayashi 
& Zhang, 2005; Eldridge et al. , 2006; Pe'er & Wijers, 2006; Nakar & Granot, 2006): This variability source 
is discussed mostly in the context of turbulent ISM or massive stellar wind environment. In the case of 
SHBs massive stellar winds are most likely irrelevant (§5), while density inhomogeneities resulting from ISM 
turbulence are expected to have a low amplitude (of order unity). The main distinction of density induced 
variability is that it is chromatic. The light curve varies both above and below i^c, but these variations are 
uncorrelated (Nakar & Granot, 2006). The resulting fluctuations can be either 'bumps' or 'dips' (relative to 
the unperturbed light curve). A turbulent ISM can produce only rapid, low amplitude, fluctuations on top of 
a smooth power law decay (Wang & Loeb, 2000). Significant density fluctuations, which are not expected in 
the case of SHBs, can produce high amplitude light curve variability but only over a long duration {St ^ t). 
Detailed limits on the light curve fluctuations in the case of increasing density are discussed in Nakar & 
Granot (2006). 

(ii) angular variability of the ejected energy ("patchy shell" e.g., Kumar & Piran, 2000b; Nakar, Piran & 
Granot , 2003; Nakar & Oren, 2004; loka, Kobayashi & Zhang, 2005): If the outflow energy has an intrinsic 
angular structure it will produce a variable afterglow light curve. As the blast wave decelerates, the angular 
size of the observed region (~ I/7) increases. As a result, the effective (average) energy of the observed 
region, and hence the observed flux, vary. If the angular structure has a typical scale then the fluctuation 
amplitude decreases with time as more patches enter the observed region. Averaging over larger random 
structures leads to a decay of the fluctuation envelope as t~^^^ (Nakar, Piran & Granot , 2003; Nakar & 
Oren, 2004). A surprising result is that although the variability is not limited by the angular time, the 
gradual change in the visibility of the patches dictates that the variability time scale is At t (Nakar & 
Oren, 2004). An important feature in this scenario is the break of the axial symmetry and therefore the 
possible production of linear polarization. Both the amplitude and the angle of polarization are correlated 
with the light curve variations (Granot & Konigl, 2003; Nakar & Oren, 2004). 

{\\\)Energy injection (e.g., "refreshed shocks" Rees & Meszaros, 1998; Kumar & Piran, 2000a; Sari & 
Meszaros, 2000; Bjornsson et al., 2002; Granot, Nakar & Phan, 2003; Fox et al. , 2003; Nousek et al. , 
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2006; Zhang et al. , 2006; Granot & Kumar, 2006): Energy can be injected into the external blast wave by 
late-arriving shells. Such shells can be produced by late engine activity, or a slow shell, ejected during the 
initial burst, can overtake the decelerating material behind the (initially faster) afterglow shock. In both of 
these cases the blast wave energy increases monotonically and, therefore, the observed flux is expected to 
rise above the unperturbed light curve, in all wavelengths, during energy injection. When the injection stops 
the original decay is resumed. The resulting light curve in the case of discrete episodes of energy injection 
has a distinctive step-like structure. The time scale of the steps (St) depends on their timing relative to the 
jet-break. St ^ t before the brake and St < t is possible after the jet-break (Granot, Nakar & Piran, 2003). 
In any case St is limited by the time associated with the width of the shell behind the afterglow shock, 
^A, implying St > t/10. A variant of this model is the two-component jet model (e.g., Berger et al. , 2003; 
Konigl, 2004; Peng, Konigl & Granot, 2005; Granot, 2005; Wu et al. , 2005; Jin et al. , 2007) where the 
opening angel of the slow ejecta is larger than the opening angel of the fast ejecta. 

(iv) Microlensing event (Loeb & Perna. 1998; Garnavich, Loeb & Stanek, 2000; Granot & Loeb, 2001; 
Mao & Loeb. 2001; Gaudi, Granot & Loeb, 2001; Koopmans & Wambsganss, 2001; loka & Nakamura, 
2001; Baltz & Hui, 2005): The apparent radius of the afterglow (^ R/^) increases with time (as t^^^ in the 
quasi-spherical phase) and at cosmological distances its angular size after a day is ~ ^as, comparable to 
the Einstein radius of a solar mass gravitational lens at these distances (Loeb & Perna, 1998). Therefore, 
a microlensing event may induce afterglow brightening, which is expected to be seen simultaneously in 
all wavelengths. Some color dependence is expected as a result of the color dependent afterglow surface 
brightness (Granot, Piran & Sari, 1999a; Granot & Loeb, 2001). Note that the probability for an afterglow 
microlensing is small (Koopmans & Wambsganss, 2001; Baltz & Hui, 2005) and therefore only a small 
fraction of the observed afterglow variability may be attributed to microlensing. 

All well-sampled SHB X-ray afterglows show strong variability. In almost all cases St ~ t, implying that 
the source of the variability is not necessarily late engine activity (although it may be). An exception is SHB 
050709 where Fox et al. (2005) find that during the last Chandra observation (exposure time of 5 hr, taken 
at i = 16 d) most of the flux arrived during the first 6 ks (a constant emission rate is rejected at 99.9% 
confidence), suggesting a high amplitude flare, AF/F w 10, with very rapid variability St/t w 0.005. Such 
a flare, if real, requires a source that is still active 16 days after the burst. Note that the energy emitted in 
this X-ray flare was ~ 10*^ erg, about four orders of magnitude below that of the GRB itself. 

Large flares observed in other bursts (e.g., SHBs 050724, 051221 and 060313) are too bright and rapid to 
be explained by external density enhancement (Nakar & Granot, 2006) and most likely result from energy 
fluctuations of the external shock (either a patchy shell or refreshed shocks) or by internal dissipation of 
flow ejected during late engine activity. Note that SHB 051221 shows the very indicative step-like light curve 
expected by refreshed shocks (Burrows et al. , 2006; Soderberg et al. , 2006a) while flux in SHB 050724 
afterglow returns to its original level after it flares (assuming a constant underlying decay; Campana et al. 
, 2006; Grupe et al. , 2006) suggesting that the flares result from hot spots on the blast wave, or that they 
are unrelated to the external shock (i.e., late engine activity). Interestingly, similar variable afterglows are 
observed in long GRBs, such as the step- like light curve of GRB 030329 (e.g., Lipkin et al. , 2004) or the 
X-ray flares after which the original underlying flux level and decay rate are resumed in many Swift X-ray 
afterglows (e.g., Nousek et al. , 2006). 

4.8. The early X-ray "tail" 

Currently, two Smft and HETE-2 SHBs (050709 & 050724) exhibit soft X-ray tails, which last for ~ 100 s 
and are separated by ^ 10 s from the short and hard bursts of prompt gamma-rays (§2.2.2). The variability 
of such a tail is an important clue to its origin. Similarly to the prompt emission, rapid high amplitude 
variability implies continuous engine activity, while low variability emission can be explained by external 
interaction. In these two bursts no apparent rapid variability is seen, although the low signal to noise ratio 
in these two bursts may hide such variability, if it exists. Among the several BATSE bursts identifled by 
Norris & Bonnell (2006) to have a similar temporal structure, there are some in which the soft tail is highly 
variable. Similarly, the X-ray tail of the nearby GRB 060614, which may be an extreme member of the 
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SHB family (Gal- Yam et al. , 2006; Gchrels et al. , 2006), is also highly variable. Therefore, while at this 
point continuous engine activity is not absolutely necessary to explain the observed X-ray tails, there are 
suggestive hints in this direction. Note that X-ray flares observed hours and days after the bursts, might 
require late engine activity as well. 

If the X-ray tail can be explained by external dissipation then it is intriguing that its duration and timing 
are similar to the expectation from reverse shock emission (Eq. 30). According to the canonical microphysical 
assumptions, reverse shock emission is not expected to contribute to the X-ray emission (§4.3). However, it 
is possible that the composition of the outflow (e.g., its magnetization level) can affect particle acceleration 
or the magnetic field of the shocked plasma in a way that results in bright X-ray emission. For example, 
a large e b together with effective acceleration of only a small fraction of the electrons (resulting in a high 
value of 7m) can produce soft X-ray emission that begins ~ 10 s after the prompt emission, lasts for ^ 100 
s and has a comparable energy to the one observed in the prompt emission. 

MacFadyen, Ramirez-Ruiz & Zhang (2005) suggest a binary SHB progenitor that is composed of a neutron 
star and a non-compact companion (§5.3.1). In this model the SHB takes place once the neutron star collapses 
to a black hole (driven by accretion). The X-ray tail is produced by interaction of rclativistic ejecta with 
the non-compact companion. The time scale as well as the typical emission frequency and the total radiated 
energy from this interaction can fit observed X-ray tails. This model, however, can explain only a single 
smooth episode. 

If however X-ray tails are produced by extended engine activity then there are several different schemes 
that were suggested in order to extend the duration of hyperaccretion of ~ 0.1 onto a black hole (the most 
popular central engine model) beyond the natural timescale of < 1 s (§5.2.1). Examples include magnetic 
barrier that chokes the accretion (van Putten & Ostriker, 2001; Proga & Zhang, 2006), disk fragmentation 
(Perna, Armitage & Zhang, 2006) and a complex evolution of the neutron star disruption in a NS-BH 
binary (Faber et al. , 2006a). Alternatively, viscosity that is lower then expected can prolong the accretion 
time (§5.2.1). Fan, Zhang & Proga (2005) argue that if the X-ray tail results from a low accretion rate, 
the energy of the accretion cannot be extracted by neutrino-anti neutrino annihilation, since the efficiency 
of this process falls rapidly when the accretion rate is lower than 0.01 Mq/s. For a detailed discussion of 
accretion disk lifetime and the efficiency of different energy extraction mechanisms see §5.2.1. A different 
solution may be that the central engine is a hyper-magnetized neutron star (e.g., Usov, 1992, ;§5.3.2 and 
§5.2.11). The activity duration of this type of engine may be significantly longer than several seconds (e.g. 
Gao & Fan, 2006; Fan & Xu, 2006; Dai et al. , 2006). 

4.9. Macronova 

An interesting possible source of late time emission is radioactive decay of matter that was ejected sub- 
relativistically during the burst. There are several possible scenarios in which such ejecta arc launched along 
with the rclativistic outffow that produce the SHB. For example, during the coalescence of a neutron star 
with another compact object (neutron star or a black hole) considerable amount of mass, of the order of 
O.OIM0, may be ejected from the system (e.g., Rosswog et al. , 1999; Ruffert & Janka, 2001). During the 
merger the density and temperature of the ejected mass can be high enough to produce heavy neutron 
rich elements, most of them radioactive (e.g., Lattimer & Schramm, 1974, 1976; Symbalisty & Schramm, 
1982; Eichler et al. , 1989; Rosswog et al. , 1999; Ruffert & Janka, 2001). Another example of a source of 
radioactive sub-relativistic ejecta may be nucleosynthesis in the wind that is driven by a neutrino-dominated 
hyper-accreting disk, which is the leading central engine model (e.g., MacFadyen, 2003). 

Sub-relativistic radioactive ejecta also power the emission from supernova explosions. The higher expected 
velocities in the SHB case, together with the lower amount of mass, lead to a fainter and faster evolving, 
yet supernova reminiscent, emission. Radioactive energy dominates the observed emission of sub-relativistic 
ejecta since initial thermal energy is quickly lost to expansion. The properties of the observed emission 
are determined by the interplay between the typical radioactivity lifetime (which depends on the nuclear 
composition of the outflow), and the maximal time (and radius) at which the ejecta are still dense enough 
to reprocess this energy into a blackbody radiation. If a significant amount of radioactive heat is generated 
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while the ejeeta are dense enough we may detect the emitted blackbody emission. If most of the energy is 
released at later times the radioactively generated photons may be observed directly. 

Li & Paczyhski (1998) were the first to discuss this possibility. Having the unknown composition ejected 
matter during a compact binary merger in mind, they postulate a radioactive source which decays as a 
power-law in time. Assuming that 0.01 Mq is ejected at a velocity of c/3, and that 0.1% of the rest mass 
energy is converted to heat, they find that after one day the bolometric luminosity can be as high as lO''^^ 
erg/s (Mfeoz ~ —21). Such a bright emission is already ruled out by tight constraints on the optical emission 
from SHB 050509B (Hjorth et al. , 2005a). 

Kulkarni (2005), introducing the term "Macronova" , presents detailed light curve calculations for two 
specific nuclear heat sources, free neutron decay and radioactive Ni 56. He finds that neutron-heated ejeeta 
may generate an optical emission that peaks hours after the SHB with a luminosity of ~ 10'*"'^ erg/s (M ~ 
— 14). Such emission can be detected from z ~ 0.2 by a 10-meter-elass telescope. Ni 56 heated ejeeta, on the 
other hand, are very hard to detect in optical observations, unless the velocity of the ejeeta is unreasonably 
slow (much below the escape velocity from a neutron star surface), otherwise the ejeeta becomes optically 
thin before most of the Ni 56 decay. Most of the radioactive energy in this case would escape as hard X-ray 
photons. 

To conclude, macronova emission may be detected if SHBs eject a significant amount of radioactive mass 
with a favorable decay time, and even then only by rapid (less than 1 d post burst) and deep optical 
observations. The characterizing signatures of macronova optical emission are blackbody colors and non- 
correlated optical and X-ray emission. 

4.10. Relativistic collisionless shocks 

The microphysics of collisionless shocks is a major open question in GRB physics (for a recent review 
see Waxman, 2006). The mean free path of a proton to undergo a binary collision of any type during GRB 
external or internal shocks is larger by many orders of magnitude than the size of the system. Therefore, any 
interaction between the particles must be due to collective processes. Such electromagnetic interactions are 
known to take place in anisotropic plasma and they give rise to collisionless shocks in GRBs. The external 
shock is ultra-relativistic and the "upstream" (e.g., the unshocked ISM) is unmagnetized. Internal shocks, 
if they take place, are mildly relativistic and the upstream may be magnetized by magnetic field that is 
advected all the way from the central engine. 

Afterglow observations require collisionless shocks to generate the two major ingredients of the observed 
emission: (i) The shock should generate a near equipartition magnetic field while a fiG upstream field is 
some eight orders of magnitude below equipartition (the fraction of the generated magnetic energy out of the 
total energy is parameterized by et)- The generated field must be sustained along the downstream, (ii) The 
shock should accelerate electrons to Lorentz factors that are much higher than the thermal temperature of 
the "downstream" (the shocked ambient medium). The outcome of the acceleration process is parameterized 
by £e and p (see §4.1). 

There is currently no complete "first principles" theory that explains how collisionless shocks accelerate 
particles or generate and sustain the magnetic field, e^, Eg and p are simply phenomenological (almost 
certainly oversimplified) model parameterizations, which are constrained by the observations. In long GRBs 
these parameters are reasonably constrained (e.g., Panaitescu & Kumar, 2001; Yost et al. , 2003) and so is 
the requirement that the magnetic field is sustained long after the shock (~ 10* A^, where is the plasma 
skin depth; Rossi & Rees, 2003). In SHBs the current observations are not detailed enough to constrain 
these parameters, but the similarities to long GRB afterglows suggests that the microphysics is similar as 
well. This is not very surprising given that the microphysical processes are most likely affected only by local 
properties such as the shock velocity and the upstream magnetization and not by global properties such as 
the total outfiow energy or large scale gradients in the external density. 

Unmagnetized collisionless shocks are not well understood. Moiseev & Sagdeev (1963) suggested that 
unmagnetized Newtonian shocks are generated by the filamentation mode of the Weibel instability (Weibel, 
1959), a kinetic instability that arises when the velocity distribution of plasma particles is anisotropic. The 
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typical scale of this instability is A.,. Medvedev & Locb (1999) and Gruzinov & Waxman (1999) suggested 
that the same process may take place in relativistic GRB shocks (see however Lyubarsky & Eichler, 2006). 
The fillamentation modes of this instability (e.g.. Fried, 1959; Silva et al. , 2002; Bret, Firpo & Deutsch, 2004, 
2005) result in the breakdown of the plasma into current filaments and in the generation of a magnetic field 
in the plane of the shock. In this picture this magnetic field, which grows to equipartition levels, provides the 
effective coUisionality that produces the shock. Medvedev & Loeb (1999) suggested that this magnetic field 
is also responsible for the observed synchrotron emission. The generated field coherence scale is Xg and basic 
considerations suggest that such a field should decay over a similar scale after the shock (Gruzinov, 2001; 
Milosavljevic & Nakar, 2006b). The field coherence scale must grow quickly enough if it is to survive along 
the downstream, as suggested by Silva et al. (2003) and Medvedev et al. (2005). It is currently unclear 
whether the coherence length of the field is growing fast enough and therefore the fate of the generated field 
is unknown. 

An alternative idea is that the field is generated already in the upstream by interaction between shock ac- 
celerated protons and unperturbed upstream (Gruzinov, 2001; Milosavljevic & Nakar, 2006a). Interestingly, 
Li & Waxman (2006) used the requirement that the acceleration time of the radiating electrons (assuming 
Fermi acceleration) is shorter than the inverse Compton cooling time of these electrons in the upstream, in 
order to put a lower limit on the upstream field in front of the shock (a larger field reduces the acceleration 
time). They find that in several long GRBs the field must be larger than 0.2nQ^* mG, where n is the external 
density, suggesting that the upstream field is amplified. Repeating the same analysis for SHBs shows that no 
meaningful constraints can be derived, and a /xG field is consistent with the observations. The reason is that 
for typical SHB parameters inverse Compton cooling in the upstream is suppressed by the Klein-Nishina 
cross-section (§4.2). 

Acceleration of non-thermal particles in collisionless shocks is believed to be done through diffusive shock 
acceleration (a.k.a. first order Fermi accleration ; for reviews see Drury, 1983; Blandford & Eichler, 1987). 
Acceleration in this process takes place when a particle is scattered back and forth across the shock by 
interactions with plasma waves on each side. With each shock crossing the particle gains energy. GRB 
observations require that the acceleration process will produce electrons with a high energy power-law 
distribution with an index 2 < p < 3. Theoretical predictions of relativistic diffusive shock acceleration 
typically suggest p « 2.2 (e.g., Aehterberg et al. , 2001; Ellison & Double, 2002; Lemoine & Pelletier, 
2003; Keshet & Waxman, 2005), which fits the observations nicely'^" . However, there are some simplifying 
assumptions that go into these calculations that are not necessarily valid, and the full acceleration process 
is far from being understood. 

During the last several years there has been a large effort to simulate relativistic collisionless shocks in 
three dimensions, using the particle-in-cell (PIC) method (Silva et al. , 2003; Fonseca et al. , 2003; Frederiksen 
et al. , 2004; Jaroschek, Lesch & Treumann, 2004, 2005; Nishikawa et al. , 2005; Spitkovsky, 2005). These 
simulations show, as expected, that when two unmagnetized plasma shells collide relativistically, shocks are 
formed by the Weibel instability, building skin depth scale equipartition magnetic field in the plane of the 
shock. None of the current simulations is large enough (in time and space) to achieve a steady state shock 
or to explore the decay of the generated magnetic field in the downstream, due to computational resource 
limitations. None of the simulations show a clear indication of diffusive shock acceleration. This is not very 
surprising, given the simulation size and the fact that the upstream field is practically zero. The reason is 
that in these conditions there is no upstream scattering agent that can reflect high energy particles back 
into the shock and start the acceleration process. In reality the pre-existing upstream (e.g., ISM) magnetic 
field is necessary to enable reflection of particles back into the shock over scales that are much larger than 
the skin depth. 



^"The exact value of p depend on the details of the interaction between the accelerated particles and the magnetic turbulence 
(e.g., Ellison & Double, 2004; Lemoine & Revenu, 2006; Katz, Keshet & Waxman, 2007; Keshet, 2006). 
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5. Progenitors and the central engine 

In my mind, the most interesting open SHB question today is the nature of the progenitor and the 
processes leading to the formation of the central engine (for a recent review see Lee & Ramirez-Ruiz , 2007) . 
This branch of research has been active for many years and some of the models were explored in great detail. 
Interest in this field was boosted following the recent detection of SHB afterglows, as we now have an idea 
about the typical environment and redshift in which SHB progenitors dwell. An immediate implication of 
this information is the understanding that SHB progenitors are different systems than the progenitors of long 
GRBs, providing the proof that long and short GRBs are two distinct physical phenomena. Nevertheless, 
the identity of the progenitors of SHBs is still unknown. 

Coalescence of a compact binary, either a double neutron star or a neutron star (NS) and a black hole 
(BH) binary, is currently the leading progenitor candidate. It is also the most extensively explored model 
and at least NS-NS systems are observationally known to exist. Therefore, this model will be the focus 
of this chapter (§5.2). This progenitor candidate is of great interest since NS-NS and NS-BH coalescence 
events are the most promising sources of gravitational-wave (GW) signals accessible to ground-based GW 
observatories. The evidence supporting the merger model are indirect, and some aspects of this model seem 
to be 'uncomfortable' with the observations. Therefore, it is important at this point to keep an open mind 
to the possibility that SHBs may have a different origin than compact binary mergers. Other models are 
discussed in §5.3. 

The debate about the origin of long GRBs ended only with the detection of SN spectral signatures in 
optical long GRB afterglow (Stanek et al. , 2003; Hjorth et al. , 2003). Lacking similar direct evidence 
for SHBs, the next best thing is an examination of indirect environmental properties. The hope is that this 
information, together with the requirement that the progenitor will be, in principle, able to produce an SHB, 
will single out viable progenitor systems. Therefore, before discussing specific models and their predictions, 
I review the new constraints that the detection of SHB afterglows pose on progenitor properties. 

5.1. Progenitors lifetime and intrinsic rate 

SHB progenitor lifetime is constrained using two independent methods - (i) based on the observed redshift- 
luminositiy distribution, and (ii) based on the host galaxy type distribution and the fractional association of 
host galaxies with galaxy clusters. The results from both methods are similar - SHB progenitor population 
is dominated by long-lived (^ several Gyr) systems. Unfortunately, the current sample available for this 
kind of analysis is small, and it is unlikely that its size will significantly grow in the near future (see below). 
Therefore, the conclusions that are based on this sample should be considered with care. 

5.1.1. The sample 

The constraints on SHB progenitor lifetime and on the intrinsic SHB rate are statistical by nature. 
Therefore, the sample that is used should be chosen carefully, and selection effects should be well understood. 
When the first analysis exploring this topic was carried out (the end of the summer of 2005; Nakar, Gal- Yam 
& Fox, 2006; Guetta & Piran, 2006) all four SHBs localized with arcsecond precision were associated with 
host galaxies (SHBs 050509, 050709, 050724) or with a galaxy structure (SHB 050813). Therefore, these 
four bursts constitute a complete sample in the sense that no events were left out due to possible selection 
effects. Nakar, Gal- Yam & Fox (2006) also carried out an analysis based on an extended sample of SHBs, 
including four IPN localized bursts (§2.3; Gal- Yam et al. , 2005). These constitute a complete sample of IPN 
bursts selected based on the size of their error boxes and locations far from the Galactic plane. 

Unfortunately, the hosts and the rcdshifts of SHBs discovered after these first well-localized four bursts, 
were not identified (in some cases even an early X-ray afterglow was not detected), in part as a a result 
of diminished observational effort. As a result, any additional SHB for which the redshift and/or host is 
determined (e.g., SHB 051221) cannot be currently included in an unbiased sample. The reason is that 
selection effects that prevent easy identification of the hosts of these bursts may strongly bias any sample 
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that docs not compensate for the undetected bursts. HopefuUy, it wiU be possible in the future to assemble 
a larger unbiased sample in order to obtain more definite conclusions. 



5.1.2. Constraints from SHB redshift distribution 

The observed redshift-luminosity two-dimensional (2D) distribution is determined by the intrinsic redshift 
and luminosity distributions, modified usually by the detector sensitivity. Thus, the goal of the analysis is 
to use the observed distribution to constrain the intrinsic redshift distribution (Piran, 1992; Gal- Yam & 
Maoz, 2004; Ando, 2004; Guetta & Piran, 2005, 2006; Nakar, Gal- Yam & Fox, 2006). Then, the lifetime 
of the progenitor can be constrained by making the most probable assumption that the progenitor birth 
rate follows the cosmic star-formation history, and that the intrinsic redshift distribution is a convolution 
of the cosmic birth rate with the lifetime distributions. The main disadvantage of this method is that for 
small observed samples the intrinsic distributions are under-constrained. It is only possible to assume several 
physically motivated parametric functional forms of the intrinsic distributions and constrain their parameter 
values. 

The two-dimensional observed redshift, z, and luminosity, L, distribution is derived from the intrinsic 
distributions via: 

where Nobs is the observed SHB rate and 0(i) is the intrinsic peak luminosity function, which is assumed 
to be independent of z. < S{P) < 1 is the probability for detection, including redshift determination when 
redshift information is needed, of a burst with a peak photon flux P, which in turn depends on L and z 
as well as on the spectrum of the bursts. Rshb{z) is the intrinsic SHB rate per unit comoving volume and 
comoving time. Since SHB progenitors are most likely of a stellar origin it is expected that: 



/■°° dt 
Rshb{z)oc J SFR{z')f{t{z)-t{z'))—dz', 



(37) 



where SFR{z) is the star formation rate at redshift z (per unit comoving volume and comoving time), t(z) 
is the age of the universe at redshift z, and f{T)dT is the fraction of SHB progenitors that are born with a 
lifetimes between r and t -|- dr. 

S{P) can describe a single detector or a combination of several detectors, each weighted by its field of 
view and operational time. In principle if S{P) is well known and if the observed sample is large enough 
then the intrinsic distributions can be extracted from Eq. 36. In reality we have to work with a limited 
sample as well as poorly understood S{P). In the case of GRBs (long and short) there is a large sample of 
bursts, observed by BATSE , for which only the peak flux distribution is available while the redshift (and 
thus luminosity) is unknown. The BATSE sample can constrain the intrinsic distributions by considering 
the observed flux distribution which is an integration of Eq. 36: 

' f^dzT dL^^, (38) 



dP dP Jq jl(z,p) dLdz 
where 

L{z,P) = Andlk{z)P (39) 

is the luminosity for which a burst at a redshift z has a peak flux P. dL{z) is the luminosity distance and 
k{z) depends on the spectrum of the bursts and includes the k-correction as well as the conversion from 
energy flux to photon flux. k{z) is assumed to be a function of the redshift only. 

If </>(i) is a single power-law, 4'{L) = (j)oL~^ , with no upper or lower cutoff (within a luminosity range 
that is discussed below) then the integral over z in Eq. 38 does not depend on P and thus the observed peak 
flux distribution docs not depend on Rshb arid simply satisfles: 

oc P-PS{P). (40) 

Note that since Eq. 40 is independent of the burst redshifts S{P) here is the probability for detection alone. 
Similarly, the integral over L in eq. 36 results in: 
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^ = (4.difc(z))-^0o^^ / P-^S{P)dP, (41) 
az I + z az J 

thereby eliminating the dependence on S{P) up to a constant normahzation factor. Naturally, for f3 < 2 
an upper cutoff must exist while for /3 > 1 a lower limit is necessary. However, if the lower cutoff is low 
enough so that it affects only a negligible volume, and if the upper cutoff is high enough so it affects only the 
detection at high rcdshift, then Eqs. (40) and (41) are applicable (these cut-offs also prevent the integral over 
P in Eq. 41 from diverging). Therefore, if the observed peak flux distribution can be fitted by Eq. 40 then 
the luminosity function can be a single power-law. In this case data sets for which S{P) is not well known 
can be readily used and Eq. 41 enables a comparison of the one-dimensional observed redshift distribution 
with model predictions. Unfortunately the observed luminosity distribution depends on S{P) even when 
the luminosity function is a single power-law. If S{P) is well known, a better constraint on the intrinsic 
distributions can be obtained by a comparison with the two-dimensional luminosity-rcdshift distribution 
(Eq. 36). 

The observed BATSE peak flux distribution, ^^^^^p^, can be fitted successfully by a single power-law 
with an index f3 = 2 ± 0.1 (Nakar, Gal- Yam & Fox, 2006). Therefore, it is consistent to assume a single 
power-law luminosity function and to use Eq. 41 on a sample of bursts from different detectors with unknown 
S{P) - altogether 8 bursts (4 Swift / HETE-2 bursts and 4 IPN bursts). Figure 11 (from Nakar, Gal- Yam 
& Fox, 2006) presents the observed cumulative redshift distribution comparing to the predictions of several 
models (Eqs. 37 & 41) taking the star formation history formula SF2 from Porciani & Madau (2001) and 
using different functional forms of life time distribution: (i) lognormal, f{T)dT ~ (TcrA/27r)^^ea;p[— (Zn(T) — 
Inir^))'^ j1cP'\dT with various values of r^, and narrow (a = 0.3) or wide (cr = 1) dispersions; (ii) power-law 
distributions /(t) oc with a lower cutoff at 20 Myr and an upper cutoff that is larger than the Hubble 
time. This figure vividly illustrates that models with a typical delay are consistent with the data only if this 
delay is long (> 4 Gyr). Models with no typical lifetime (power-law distributions) must have a birthrate of 
progenitors per unit logarithmic lifetime that increases significantly as a function of lifetime (i.e., 77 > —0.5). 
The reason that models dominated by short-lived systems do not fit the data is that they under-predict the 
fraction of bursts at low redshift (z < 0.3). Since the star formation rate is higher at large redshifts (> 1), 
while SHBs arc found preferentially at low redshift (w 0.3), although they arc bright enough to be detected 
at high redshfit, implies that a long delay is required for the progenitors so they arc born when the universe 
is young and produce SHBs only when it is much older. Quantitatively, Nakar, Gal- Yam & Fox (2006) find 
that for lognormal lifetime distributions (narrow and wide) r* > 4[1] Gyr in 95 [> 99.9]% confidence and 
for a power-law lifetime distribution the most probable power-law index is 77 = 0.6 and 77 > — 0.5[— 1] at 
95[99.5]% confidence '^^ . These results are similar to those obtained by Guetta & Piran (2006). 

The fact that the observed luminosity function fits a single power-law does not necessarily imply that 
the intrinsic luminosity function is a single power law as well. When an intrinsic luminosity function that is 
not a single power-law is considered, the threshold of the detector must be known (at least partially) and 
therefore the sample that can be used is reduced only to the first three Swift SHBs. However, in such case 
the full two-dimensional L — z distribution can be considered. Nakar, Gal- Yam & Fox (2006) find that when 
a 'knee shaped' broken power law (i.e., 4>{L) is steeper above the break) lifetime distribution is assumed 
the results are very similar to the case of a single power-law. These results suggest that the observations 
prefer long-lived progenitors (Gyrs) for unimodal lifetime distributions. If a bimodal luminosity function 
is considered then there are solutions in which the progenitor lifetime is dominant by short-lived systems. 
However, such luminosity functions are not expected unless SHBs are composed of two separate populations. 

This analysis does not include the recent indication found by Berger et al. (2006b) that at least 1/4 of 
the SHBs are at z > 0.7 (assuming that we are not fulled by SHBs that take place at large distances, ^ 100 
kpc, from there hosts). If confirmed, this result indicates that the lifetime of a non-negligible fraction of the 
SHB progenitors is rather short (Gyr or less), implying that the lifetime distribution is wide (includes both 

These values are calculated assuming that SHB 050813 is at z=0.72 (Gladdcrs et al. , 2005; Berger, 2005; Prochaska et al. 
, 2005). If a redshift of 1.8 is taken instead (Berger, 2006) then a narrow lognormal lifetime distribution is ruled out (cannot 
explain both the low and high redshift bursts). For a wide lognormal lifetime distribution t, > 2[1] Gyr in 95[> 99]% confidence 
and for a power-law lifetime distribution the most probable power-law index is = and rj > — 0.8[— 1] at 95[98]% confidence. 
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Fig. 11. The cumulative observed redshift distribution as predicted by various lifetime distributions when the luminosity 
function is <f>(L) oc and the star formation history is given by SF2 formula in Porciani & Madau (2001). The cumulative 
redshift distribution of the observed data (shaded area) is bracketed between the lower solid line, which is the cumulative 
redshift distribution of the six bursts with known redshifts, and the upper solid line, which includes also the contribution of the 
two bursts with upper limits (see Nakar, Gal- Yam & Fox (2006) for details). The figure demonstrates that only models which 
arc dominated by long-lived progenitors provide a good fit to the data while models with short-lived progenitors under-predict 
the fraction of bursts at low redshift z < 0.3. [taken from Nakar, Gal- Yam & Fox (2006)] 

old and young progenitors). It rules out a narrow lognormal {a = 0.3) lifetime distribution or a power-law 
distribution with > 0. A wide lognormal (cr ~ 1) lifetime distribution with 4 < < 8 Gyr and a power-law 
distribution with —0.5 ^ ^7 ^ arc consistent with all current observations. 

5.1.3. Based on host galaxy types 

An alternative method to estimate the progenitor lifetime is based on the spectral types of the host galaxies 
(Gal- Yam et al. , 2005). The stellar populations in early type galaxies (E & SO) is entirely dominated by 
old systems (several Gyr) while the population in galaxies of later types is composed of a mix of old and 
young systems. Therefore, the fraction of SHBs in different types of galaxies is an independent measure of 
the progenitor lifetime. Gal- Yam et al. (2005) compared the distribution of SHB host types to that of type 
la supernovae (SNe) from Mannucci et al. (2005), building on the extensive effort done in undertaken to 
constrain SNe la progenitor lifetime (Fig. 12). They find that a larger fraction of SHBs, compared to SNe 
la, take place in early type galaxies (at 93% confidence), implying that it is most likely that SHBs are older, 
on average, than SNe la. The typical time delay of SNe la, r/a is still under debate, but all different models 
that assume a unimodal lifetime distribution agree that r/a > 1 Gyr (Tonry et al. , 2003; Barris & Tonry, 
2006; Strolger et al. , 2004; Gal- Yam & Maoz, 2004; Maoz & Gal- Yam, 2004; Mannucci et al. , 2005) and 
therefore SHBs, which are most likely older, have a lifetime of several Gyr. It was recently suggested that the 
population of SNe la is composed of two sub-groups - a long-lived component which is responsible for the 
entire SN la population in early type galaxies, and a short-lived component, which follows the star formation 
rate (Mannucci et al. , 2005, 2006; Scannapieco & Bildsten, 2005). If this is the case, the distribution of 
SHB host types is consistent with that of the long-lived SN la component which is several Gyr old. 
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Fig. 12. A comparison between the host galaxy types (E/SO, Sa/Sb, Sbc/Sc and Sd/I) of SHBs (from Table 2) and SNe la 
(Mannucci et al. , 2005). The fraction of SHBs in early type galaxies is significantly larger than the fraction of SNe la observed 
in such galaxies in the nearby Universe, indicating that the progenitor systems of SHBs are probably longer-lived than those 
of SNe la. [form Gal- Yam et al. (2005)] 

Zheng & Ramirez-Ruiz (2006) constrain the hfctime of SHE progenitors directly by decomposing the total 
star formation history to the contribution from early- and late-type galaxies. The convolution of the split 
star formation history and a given lifetime distribution predicts the fraction of SHBs in each host type. 
Zheng & Ramirez- Ruiz (2006) find that if the lifetime distribution is a single power law, /(t) oc , then the 
observations suggest ?/ > 3/2 implying a typical lifetime of ^ 10 Gyr (note that this result is inconsistent 
with the recent observations of Bcrgcr et al. (2006b) that suggest that tj <0). 

Shin & Berger (2006) follow a similar route to that of Zheng & Ramirez- Ruiz (2006), but rather than using 
the relative SHB rate in early- and late- type galaxies, they use the relative rates in cluster and field early-type 
galaxies. The advantage of this method over the early/late ratio is that the star formation histories of early- 
type galaxies are simpler. The disadvantage is that this method does not sample the lifetime distribution 
below ~ 1 Gyr since short-lived SHBs would not take place in early-type galaxies. Shin & Berger (2006) find 
that the current ratio of cluster to field early-type galaxies, 2, corresponds to < 77 < 1, which is consistent 
with earlier results from redshift distribution and early- type/late- type ratio analysis. 



5.1.4. Intrinsic local rate 

The observed BATSE local rate is TZsHB.obs ~ 10 Gpc'^ yr"^ (§2-6). This is clearly a lower limit on the 
true intrinsic rate. An upper limit on the SHB rate can be derived under the most probable assumption that 
the progenitors of SHBs are a product of at least one core-collapse supernova (e.g., a neutron star) and that 
the SHB itself is a catastrophic non-repeating event. In this case the rate of SHBs is bound by the rate of 
core-collapse SNe. However, since the typical lifetime of SHB progenitors seems to be several Gyr, the local 
SHB rate corresponds to the rate of core-collapse SNe at a redshift « 0.7, which is ~ 5 x 10^ Gpc^^ yr^^ 
(Dahlen et al. , 2004). Therefore, the intrinsic local rate, TZshb = Rshb{z = 0), must be within the range: 

10 < UsHB < 5 X 10^ Gpc-3 yr-i. (42) 
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Constraining the local rate within this range requires an estimate of the fraction of undetected bursts. 
These may be events that are coUimatcd and point away from the observer, or bursts that are simply too 
dim to be detected. The number of bursts that are beamed away from us can be estimated by the beaming 
factor (i.e.; the fraction of the total solid angle into which the prompt gamma-rays are emitted), /^^^. As 
I discuss in §4.6, the value of the beaming factor is reasonably constrained only in two bursts. In one case 
(SHB 051221A) it seems to be w 80 (assuming a 'top hat' jet) and in the other (SHB 050724) it is probably 
much smaller. Thus it is impossible at this stage to confidently estimate the beaming factor, but current 
observations suggest that 1 ^ /(, ^ 100. The correction for undetected dim bursts depends on the luminosity 
function and most strongly on its lower cutoff, Lj-nin (defined so (t>{L) = for L < Lmm)- Assuming the 
single power-law luminosity function that is found to be consistent with the data {4>{L) oa L^^ L > Lmin) 
the local rate is (Nakar, Gal- Yam & Fox, 2006): 

5.2. Coalescence of a compact binary 

The coalescence of two neutron stars was recognized as a potential GRB progenitor already two decades 
ago. This possibility was briefly mentioned in Blinnikov et al. (1984), Paczynski (1986), Goodman (1986) 
and Goodman, Dar & Nussinov (1987), and discussed in detail for the first time by Eichler et al. (1989) 
and later by Narayan, Paczynski & Piran (1992) and many others. Paczynski (1991), Narayan, Paczynski & 
Piran (1992) and Mochkovitch et al. (1993) discussed similar coalescence models of a neutron star (NS) and 
a stellar mass black hole (BH) as a possible progenitor. Binary mergers take place because of orbital energy 
and angular momentum loss to gravitational-wave radiation, as was observationally confirmed (Taylor & 
Weisberg, 1982). These models are natural candidates since such events must take place at a reasonable 
rate (see below) and the amount of gravitational energy that is liberated during the coalescence (~ 10^^ 
erg) is large enough, so a small fraction of it can generate the GRB. Additionally, in the 'standard' version 
of the merger model, the outcome of the coalescence is the same central engine as in long GRB models - a 
disk accreting onto a black hole. This similarity between the engines naturally explains the great similarity 
between the observational properties of the two phenomena. In both scenarios the duration of the burst is 
determined by the lifetime of the disk. The difference is that in binary mergers the disk is expected to be 
consumed within a fraction of a second, while in the 'coUapsar' model for long GRBs, in-falling matter from 
the collapsing star feeds the disk for much longer, enabling the production of a long duration GRB. 

5.2.1. The "central engine" 

I. An accretion driven engine 
In the most popular central engine scenario a binary merger results in the formation of a disk around a 
black hole (newly born in the NS-NS case). Accreting black holes are known to produce relativistic jets in 
other systems such as active galactic nuclei and microquasars, however the accretion rates that are required 
in order to power GRBs are higher by many orders of magnitude, implying different physical conditions. 
The models of the formation and operation of this type of SHB engine are discussed below. 

/./ Forming an accretion disk in a NS-NS merger 
An extensive effort was dedicated to numerical simulations of NS-NS mergers in the context of the GRB 
central engine. Starting from 3-dimensional Newtonian simulations with a polytropic equation of state (EOS) 
and no neutrino effects (Davies et al. , 1994), the simulations evolved to include increasingly more realistic 
EOS, neutrino modeling and asymmetric scenarios. These simulations were carried out mainly by two groups 
using different numerical schemes - Lagrangian Smooth Particle Hydrodynamics (SPH) (Rosswog et al. , 
1999, 2000; Rosswog & Davies, 2002) and Eulerian Piecewise Parabolic Method (PPM) (Ruffert, Janka & 
Schaefer, 1996; Ruffert et al. , 1997; Ruffert & Janka, 1998, 1999, 2001). The outcomes of these simulations 
are rather similar. A central quasi-axissymmetric object is formed within a few orbital periods after the 
initial contact. This central object is surrounded by a thick disk of material at a radius of tens of km with 
densities of ^ 10^^ — 10^^ gr cm~^ and temperatures of ^ 1 — 10 MeV. The mass of the disk depends 
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somewhat on the initial conditions (e.g., spin, mass ratio) but is typically in the range of 0.03 — O.SA/q for 
a system with an initial mass of w 3Mq. A narrow funnel along the axis of symmetry is relatively baryonic 
free and may serve as a potential site for launching a GRB jet (although it is unclear if this funnel is clean 
enough; Ruffert & Janka, 1998, 1999). 

In recent years several groups carried out approximate (e.g., Oechslin, Rosswog & Thielemann, 2002) and 
fully (e.g., Shibata & Uryii, 2000) general relativistic (GR) simulations of double neutron star merger. In the 
context of SHBs these simulations show that the general relativistic picture may be more complicated than 
the Newtonian one. Oechslin & Janka (2006) use a relativistic SPH code with conformally flat approximation 
of the Einstein field equations and a realistic EOS. For a range of initial conditions they find results that are 
rather similar to those of Newtonian simulations - a central object with a surrounding disk at a radius of 
tens of km with a mass 0.06 — 0.26^^0. Shibata and collaborators (Shibata, Taniguchi & Uryu, 2005; Shibata 
& Taniguchi, 2006) carried out a full general relativity simulation with a hybrid EOS that mimics a realistic 
stiff nuclear EOS. They find that depending on the initial mass of the system, the merger can proceed in 
two ways: (i) If the total mass is below some threshold (< 2.6Mo depending on the EOS) a hypermassive 
neutron star is formed at first (a NS that is supported by differential rotation; Baumgarte, Shapiro & 
Shibata, 2000). If this hypermassive neutron star collapses to a black hole after a sufficient delay so it can 
transport angular momentum to the surrounding matter, then the mass of the remaining disk is > O.OIMq, 
independent of the initial NS mass ratio. Adding magnetic fields, Shibata et al. (2006) carry out a full 
GR magneto-hydrodynamic simulation of the final stages of the merger starting with a hypermassive NS 
as initial conditions (see also Duez et al. , 2006). They find that magnetic breaking and magnetorotational 
instability transfer angular momentum very efficiently and as a result a disk of ^ 0.05Mq is formed by 
the time that the NS collapses to a BH. (ii) On the other hand, if the total mass of the system is above 
the threshold, the central object collapses promptly (in less than 1 ms) into a black hole, taking with it 
most of the mass and leaving a very small disk. A massive enough disk to be a possible SHB central engine 
(~ O.OIA/0) is left only in cases where the mass of the two progenitor neutron stars is significantly different 
(mass ratio < 0.8), contrary to the results of Oechslin & Janka (2006). Shibata & Taniguchi (2006) suggest 
that the most likely reason for this difference is the approximations made by Oechslin & Janka (2006) , that 
do not take into account gravitational waves radiation reaction. In the simulation by Shibata & Taniguchi 
(2006), GW radiation is one of the main channels through which the central object loses angular momentum 
and ultimately collapses. 

/.// Forming an accretion disk in a NS-BH merger 
The outcome of a NS-BH merger is less clear. The most important initial condition affecting the coalescence 
evolution is the mass ratio q = Mns/Mbh- If 1 ^ 0.1 the tidal disruption radius of the NS is within 
the innermost stable circular orbit, for any BH spin, and the NS plunges into the BH without leaving any 
residual disk (e.g., Vallisneri, 2000; Miller, 2005; Rosswog, 2005; Faber et al. , 2006a). As I describe below, 
the fate of systems with q ^ O.l seems more promising as SHB progenitor, but the outcome of this case is 
not entirely known yet. 

Numerous numerical simulations of NS-BH coalescence were carried out using Newtonian potentials. The 
main result of these simulations is that the formation of the disk depends strongly on the NS equation of state 
(EOS). In the case of soft equation of state (polytropes with F ^ 2.5 or Lattimcr and Swesty EOS; Lattimer 
& Swesty 1991) the NS is tidally disrupted during its first approach, forming a massive disk (~ O.3M0; 
Janka et al. , 1999; Lee & Kluzniak, 1999b; Lee, 2000, 2001). A more complex evolution is found in the 
case of a stiff EOS (Shen EOS [Shen et al. 1998] or a polytrope with F = 3; Kluzniak & Lee 1998; Lee & 
Kluzniak 1999a; Lee 2000; Rosswog, Speith & Wynn 2004), as seems to be required by the highest measured 
pulsar mass M = 2.1±O.2M0 (Nice et al. , 2005). These simulations show that a "stiff" NS is not disrupted 
during the first episode of mass transfer. Instead, in some cases, it bounces back into an eccentric orbit and 
the mass transfer resumes only after additional angular momentum is lost to gravitational wave radiation. 
The result is an episodic mass transfer process that lasts much longer than the duration of the simulations 
(massive disks are not formed before the simulations end). Based on these results, Davies, Levan & King 
(2005) developed a semi-analytic model of intermittent mass transfer, finding that the NS is disrupted after 
~ 1 s, leaving a ~ O.lAf0 disk. 

While newtonian simulations seem to give a qualitatively correct picture of the dynamics in at least some 
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of the NS-NS merger scenarios, it is not clear at all that they are applicable also in the case of NS-BH 
coalescence. An important general rclativistic (GR) feature which is not treated by Newtonian simulations 
is the innermost stable circular orbit (ISCO). Miller (2005) argues that the NS tidal radius may be within 
the ISCO not only in the case of low q. but also in the case of comparable NS and BH masses. The reason 
is that for q > 0.25 the mass of the NS itself begins to affect, increasing the the ISCO radius (compared to 
an isolated BH; e.g., Buonanno & Damour, 1999). For q > 0.5, even in the case of a maximally spinning 
BH, the tidal radius may be within the ISCO. Moreover, Miller (2005) suggests that even in cases where 
the tidal radius is slightly outside of the ISCO, GR effects can cause the NS to plunge directly into the BH 
without leaving a disk. 

The most advanced numerical simulation of a NS-BH merger to date were presented by Faber et al. 
(2006a, b), that use the conformal flatness approximation for GR for a Schwarzschild "^^ BH and a polytropic 
EOS with r = 2 for the NS (no simulations were carried out for stiffer EOS). First, they carry out a 
simulation with q ~ 0.1 and find that indeed the NS is 'swallowed' by the BH without being disrupted. Then 
they perform a simulation in which the tidal radius is comparable to the ISCO (equivalent to q=0.24 for their 
NS modeling ■^^ ). Here they find that instead of a direct plunge, angular momentum transfer during tidal 
disruption deposits 1/4 of the NS mass back outside the ISCO. Half of this mass becomes unbound and leaves 
the system, while the rest, « 0.15Mq, remains bound. Interestingly, only a fraction of the bound matter 
forms a hot and dense disk that is ready to perform as a SHB central engine. The remaining bounded mass, 
sa O.O5M0, is ejected to eccentric orbits and should fall back toward the BH after more than a second. At 
the end of the simulation the polar axis is not significantly polluted by baryons. These results are promising, 
however, BH spin is expected to play an important role in the dynamics of systems with 0.1 ^ q 0.5, so 
the final answer about their fate will require a full GR simulation. 

/./// The accretion rate and the lifetime of the disk 
Regardless of the initial progenitor (NS-NS or NS-BH binary) the SHB central engine in this model is a hot 
and dense torus of 0.01 — 0.3Mq that is accreted onto a stellar mass BH. The engine is active as long as 
efficient accretion takes place. Assuming that the burst is powered by accretion, the accretion rate can be 
estimated using energy requirements. The observed isotropic equivalent luminosity of SHBs is 10^*^ — 10'^^ 
erg/s. Correcting for beaming, the true luminosity is probably smaller by about one, or at most two, orders 
of magnitude. Taking a reasonable efficiency of accretion energy conversion into gamma-rays of ~ 0.01 — 1% 
(see below) implies hyper- accretion rates of 0.01 — 10A'/q/s. Are such accretion rates expected for a typical 
disk that is formed by a compact binary merger? 

The accretion rate and the lifetime of the disk in this scenario was first estimated analytically and semi- 
analytically (Popham, Woosley & Fryer, 1999; Narayan, Piran & Kumar, 2001; Di Matteo, Perna & Narayan, 
2002) . The amount of energy emitted during the burst implies that the disk must be accreted efficiently (i.e, 
most of the disk mass must be accreted and not expelled as a wind). Efficient accretion requires efficient 
cooling, which at these densities can be achived only by neutrinos. For this to occur, the temperatures and 
the densities in the disk should be high enough, which in turn implies that the radius of the disk should be 
small enough. Narayan, Piran & Kumar (2001) present an analytic approximation that is valid assuming 
that the disk is optically thin to neutrinos. They find that in order for efficient neutrino cooling to take 
place, the radius in which the disk is deposited must satisfy: 



where Rs is the Schwarzschild radius of a central BH with mass Mbh, Md is the disk mass and a is 
dimensionless viscosity parameter (Shakura & Sunyacv, 1973). As discussed above, the disk that is apparently 
formed in NS-NS and NS-BH mergers satisfies this criterion ■^^ . For a disk which is within this radius Narayan, 
Piran & Kumar (2001) find an accretion rate of: 

This approximation is exact for spherically symmetric potentials and is less appropriate for a Kerr BH. 
^^As a result of the conformal flatness approximation Faber et al. (2006a) use q=0.1 while changing the neutron star com- 
pactness (Mffg/Rffs) in order to imitate a q=0.24 case. 

''''Mergers of a BH with less compact objects such as white dwarfs or helium stars do not satisfy this criterion and therefore 
are unlikely to produce long or short GRBs (Narayan, Piran & Kumar, 2001). 




(44) 
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Mace = 0.6 (^) 

and a corresponding disk lifetime of 

/a 

^acc — 0.2 (^Q^ 

Note that the accretion time is set by the viscous time scale and is linear with l/a, which is the least 
constrained parameter in the problem. These equations assume that the disk is optically thin to its own 
neutrino emission, which for the expected accretion rates is valid for R > lORg (Di Matteo, Perna & Narayan, 
2002). A semi-analytic extension of this model to the optically thick regime was carried out by Di Matteo, 
Perna & Narayan (2002) and a further semi-analytic calculations of the the disk properties evolution in time 
is presented by Janiuk et al. (2004). Considerations of the strong magnetic fields that might be present 
in hyper-accreting disk as well as observable quantum electrodynamical effects arising from supercritical 
fields are discussed by Kohri & Mineshige (2002). Finally, General Relativistic effects of Kerr metric on the 
structure of the accretion disk are calculated by Chen & Beloborodov (2007). 

A number of detailed numerical simulations in 2D (Lee & Ramirez-Ruiz, 2002; Lee, Ramirez-Ruiz & Page, 
2004, 2005) and 3D (Setiawan, Ruffert & Janka, 2004) follow the evolution of the disk during its accretion. 
The latest simulations include detailed equation of state, approximate neutrino treatment and a viscosity 
described by an a-law. These simulations provide insights about the neutrino luminosity of the disk (see 
below) and they are in rough agreement with analytic calculations of the accretion rate and duration. 

As discussed in §3.4, the burst duration is most likely comparable to, (e.g., Kobayashi, Piran & Sari, 
1997), or at most larger by an order of magnitude than (Janka et al. , 2006; Aloy & Rezzolla, 2006), the 
duration of the engine activity. Hence, for a typical value of a = 0.1 the accretion rate as well as the burst 
duration nicely agree with those required for SHBs. However, late engine activity, which is suggested by the 
observations (see §4.7 and §4.8), is not naturally explained in this model. 

Several processes have been suggested to explain prolonged accretion that lasts longer than the viscous 
time scale, van Putten & Ostriker 2001 (see also van Putten & Levinson 2001, 2002, 2003; van Putten 2005) 
suggest a configuration in which the accretion is temporarily suspended due to a "magnetic wall" around a 
rapidly spinning black hole. This model was originally proposed for long GRBs, but it may be viable during 
the post-merger conditions of a compact binary, and explain long-lasting emission. In this model, the jet is 
energized by the spin of the BH and the accretion is suspended as long as the BH is rapidly spinning (tens 
of seconds for a 2.5M0 BH and a O.IA/q disk). Proga & Zhang (2006) suggest a different scenario in which 
magnetic flux can be accumulated near the black hole and behave as an alternating barrier, that allows for 
intermittent accretion, thereby increasing the accretion time. Perna, Armitagc & Zhang (2006) suggest that 
instabilities in the disk may cause fragmentation of the outer parts of the disk, producing separate blobs 
which slowly inspiral toward the black hole, producing late engine activity. Yet another possibility is that 
during the time that the disk is formed, some of the material is ejected into an eccentric, but bound, orbit 
and falls back into the BH after the disk was accreted. Such behavior is seen in the NS-BH simulation of 
Faber et al. (2006a) and may take place also during NS-NS coalescence, where some simulations show that 
^ 1O~^M0 are ejected from the system (e.g., Rosswog et al. , 1999; Ruffert & Janka, 2001). 

I. IV Launching the jet 

Two classes of processes were suggested to extract the energy of an accretion disk-black hole system and 
launch a relativistic jet: neutrino-anti neutrino annihilation and magnetically driven mechanisms. 
I. IV. I Neutrino driven jet 

In neutrino-cooled accretion disks, a significant fraction of the disk gravitational energy is converted during 
the accretion into neutrino flux. This energy is available to launch a jet through i/9 annihilation and subse- 
quent pair production. The possible importance of neutrino annihilation as the energy source of GRB jets was 
recognized by many authors (e.g., Goodman, Dar & Nussinov, 1987; Eichler et al. , 1989; Narayan, Paczyn- 
ski & Piran, 1992; Meszaros & Rees, 1992; Mochkovitch et al. , 1993, 1995; Witt et al. , 1994; Jaroszynski, 
1993, 1996). While neutrinos are generated in various sites during the merger (e.g., Salmonson & Wilson, 
2001), the most promising process for energy extraction is by neutrinos emitted from the cooling torus. This 
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emission lasts long enough to power the jet and is originates far enough from the black hole so it is not 
strongly affected by gravitational capture or redshift of the black hole (see Asano & Fukuyama. 2000, 2001, 
for a discussion of gravitational effects). Moreover, the geometry of the torus implies that the most efficient 
annihilation sites are along the rotation axis, which also suffer the least from baryon pollution. Therefore, 
most numerical simulations of these systems (Ruffert et al. , 1997; Ruffert & Janka, 1998, 1999; Fryer et 
al. , 1999; Rosswog & Ramirez-Ruiz, 2002; Rosswog & Liebendorfer, 2003; Rosswog & Ramirez-Ruiz, 2003; 
Rosswog, Ramirez-Ruiz & Davies, 2003; Lee & Ramirez-Ruiz, 2002; Lee, Ramirez-Ruiz & Page, 2004, 2005; 
Setiawan, Ruffert & Janka, 2004), as well as semi-analytic calculations (Popham, Woosley & Fryer, 1999; 
Di Matteo, Perna & Narayan, 2002; Janiuk et al. , 2004), were used to estimate the energy deposited by 
neutrinos emitted from the disk. 

The energy extraction efhciency is defined as the fraction of the rest mass energy of accreted matter that 
is converted into pairs by vv annihilation e^p = Li,jj/Mc^ (most of the annihilation takes place along the 
rotation axis, where baryonic pollution is low). This efficiency is a combination of the fraction of the rest 
mass energy that is emitted by neutrinos, = L^/ Mc? , and the fraction of neutrino energy that 

annihilates into e^e~ pairs, fLu^Lvu = T^vvl^'v The total efficiency depends on the accretion rate. At 
low accretion rates (M < 0.1Mq/s) the neutrino luminosity drops fast and the annihilation rate (which 
is proportional to the luminosity square) drops even faster. Therefore, at these accretion rates the total 
efficiency is very low, <C 10^** (Popham, Woosley & Fryer, 1999; Setiawan, Ruffert & Janka, 2004). On 
the other hand, disks with accretion rates that are too high (Af ^ O.IMq/s) are optically thick, and most 
of the generated neutrinos are adverted into the BH, thereby reducing f^^^j^^ while fLv^Luu remains 
roughly constant (Di Matteo, Perna & Narayan, 2002). Maximal efficiency is obtained at M ~ 1Mq/s were 
fhi^Lu ^ 0-01 - 0-1 and fLu^Luu ^ 10"^ - 10"^ (e.g., Di Matteo, Perna & Narayan, 2002). Therefore, 
the total efhciency of neutrino annihilation under optimal conditions is e^i; ~ 10^^. In the context of SHB 
central engine this optimal accretion rate is also the expected one, implying that the total energy that can 
be expected from this process is ~ lO^^(Md/O.lM0) erg where Md is the disk mass. This energy is enough to 
power luminous SHBs (e.g., SHB 051221) only if their flow is narrowly coUimated '^^ and less luminous SHBs 
(e.g., SHB 050509B) also if they are not narrowly beamed. Current limits on the total energy of SHBs (Table 
2) are still compatible with this model, but future observations may require higher energy. Ramirez-Ruiz & 
Socrates (2005) suggest the existence a hot corana that modifies the spectrum of the emitted neutrinos as 
a way to significantly increase fhv-tLiyu- 

I. IV. II Magnetically driven jet 
An alternative method to launch relativistic jets from the BH-accretion disk system is via strong elec- 
tromagnetic fields (e.g., Narayan, Paczynski & Piran, 1992; Levinson & Eichler, 1993; Thompson, 1994; 
Ruffert et al. , 1997; Meszaros & Rees, 1997b; Lee, Wijers & Brown, 2000; Brown et al. , 2000b; Rosswog 
& Ramirez-Ruiz, 2002; Rosswog, Ramirez-Ruiz & Davies, 2003; Daigne & Mochkovitch, 2002a; Lyutikov, 
2006). Even if the initial magnetic field is low, it is expected to play an important role due to amplification 
by the magneto-rotational instability (for review of MHD accretion disks see Balbus & Hawley, 1998). An 
extensive numerical effort is invested in systems of magnetized accretion onto BH, with rapid progress from 
simulations that use newtonian and pseudo-Newtonian potentials (e.g., Hawley, 2000; Hawley & Krolik, 
2002; Hawley & Balbus, 2002; Armitage, Reynolds & Chiang, 2001; Armitage & Reynolds, 2003; Proga, 
2003; Proga & Begelman, 2003; Machida & Matsumoto, 2003; Kato, Mineshige & Shibata, 2004) to GR 
simulations using a Kerr metric in two dimensions (McKinney & Gammie, 2004; McKinney, 2005, 2006) 
and in three dimensions (De Villiers, Hawley & Krolik, 2003; De Villiers et al. , 2005; Hirose et al. , 2004; 
Krolik, Hawley & Hirose, 2005; Hawley & Krolik, 2006). Currently, these simulations do not include neutrino 
physics and use simple equation of state (a constant adiabatic index of 4/3 or 5/3 is assumed) and therefore 
the results may not be directly applicable to the case of SHBs. Nevertheless, magnetic fields most likely play 
important, and maybe the major, role in jet launching. Different GRMHD simulations show a qualitatively 
similar picture. An accretion flow is generated within the disk by magnetic viscosity that results from am- 
plification of the magnetic field by the magnetorotational instability (Balbus & Hawley, 1991) even if it is 



Possible mechanical coUimation processes are discussed and simulated by Rosswog & Ramirez-Ruiz (2002) and by Janka et 
al. (2006) 
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initially weak. The accretion results in a strong wind of magnetized plasma which is ejected from the system 
along the boundaries of the ccntrifugally evacuated funnel, and a strong Poynting flux jet is observed within 
the funnel. The baryonic load within the Poynting-flux-dominated funnel is low enough to allow a terminal 
Lorentz factor > 100 (McKinney, 2006). 

Recent GRMHD simulations find that the energy of Poynting flux jets is directly related the spin of the 
BH and is in general agreement with the predictions of Blandford & Znajek 1977'^^ (McKinney & Gammie, 
2004; McKinney, 2005; Hawley & Krohk, 2006, see also Komissarov 2004, 2005). Energy output in Poynting 
flux jets increases sharply with the spin of the BH, and this process can be far more efficient than a neutrino- 
driven jet. Hawley & Krolik (2006) find efficiencies Cem = 3 • 10~*, 6 • 10~^, 0.04 & 0.2 for spin parameters 
a/M = 0,0.5,0.9 & 0.99 correspondingly [eem = Lem/M(? where Lem is the luminosity of the Poynting 
flux jet and a/AI = Jc/GM^fj where J is the BH angular momentum and G is the gravitational constant). 
McKinney (2006) simulates an accretion on a BH with a spinning parameter a/M = 0.9375 and finds an 
outfiow efficiency of eem ~ 0.05 where 10% of the jet luminosity is spread roughly uniformly within an 
opening angle of 0.1 rad around the polar axis. Taking an accretion disk of O.IA/q, this e^m corresponds to 
a total emitted energy of « 10^^ erg out of which 10^^ erg are within 0.1 rad of the jet axis. An observer 
with a viewing angle within 0.1 rad of the jet axis will infer a total isotropic equivalent energy of sa 10^^ erg. 
Therefore, this process, if applicable, can easily account for the energetics and the Lorentz factors observed 
in SHBs. 

//. Hypermagnetized neutron stars 
An interesting alternative GRB engine is a compact object with an ultra-high magnetic field (^ 10^^ G). 
Here, it is the energy of the magnetic field and/or the star rotation that drives a relativistic outflow, and there 
is no need for an accompanying accretion disk. Such a magnetized object can be a transient hypermassive 
NS that is formed by the coalescence of two neutron stars (Kluzniak & Lee, 1998; Rosswog & Davies, 2002; 
Rosswog & Ramirez-Ruiz, 2002; Rosswog, Ramirez-Ruiz & Davies, 2003; Shibata et al. , 2006; Duez et al. , 
2006) . Price & Rosswog (2006) carry out an MHD simulation of the merger of two neutron stars with initial 
magnetic fields of 10^^ G. They find that on very short timescales 1 ms) a magnetic field > 10^^ G is 
generated by the Kelvin-Helmholtz instability along the shear interface that forms where the two neutron 
stars come into contact. This field strength is clearly limited by the resolution of the simulations and therefore 
the actual generated field is most likely much higher. An upper limit on the fleld strength is the equipartition 
level ^ 10^^ G. If the hypermagnetized hypermassive neutron star is stable for ~ 0.1 s, as suggested by the 
results of Shibata & Taniguchi (2006), then the energy in the magnetic field can drive a relativistic outfiow 
in various ways. Price & Rosswog (2006) suggest that bubbles of equipartition magnetic field can become 
buoyant and fioat up through the star surface producing a relativistic flow (Kluzniak & Lee, 1998; Lyutikov, 
2003; Dai et al. , 2006). If a large scale field is built up then an alternative launching mechanism is the is 
the magneto-centrifugal slinging (Usov, 1992; Thompson, Chang & Quataert, 2004; Thompson, 2005). The 
emission time scale is the spin-down time of the star (assuming that it does not collapse first), which is ~ 0.1 
s in the case of a 10^^ G magnetic field (the initial spin of the hypermassive NS is 1 ms). Since the energy 
source of this engine is the magnetic field and/or the rotation of the transient neutron star (^^ 10^'^ erg), it 
can in principle produce very energetic bursts with more than 10^^ erg. Rapidly rotating hyper-magnetized 
neutron star (i.e., msec magnetar) may be produced also by other progenitor systems and it is discussed in 
some more detail in §5.3.2 

5.2.2. The lifetime of compact binaries and their merger rate 

The merger rate and typical lifetime of NS-NS binaries are estimated in two ways, based on observed 
systems in our Galaxy and using theoretical population synthesis. In the case of NS-BH binaries only the 
latter method is used, as no observed systems of this kind are currently known "^^ . These aspects of compact 



^®For a discussion of the similarities and the differences between the simulations and the BZ mechanism see Hawley & Krolik 
(2006). 

''^ The fact that no NS-BH systems are observed is not surprising given that binary detection requires a recycled pulsar, and 
that the birth rate of BH-recycled pulsar binaries is expected to be very low (Pfahl, Podsiadlowski & Rappaport, 2005). 
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binaries arc covered in detail by other papers in this volume and therefore I will briefly review here only the 
main results of each method. 

I. Constraints from observed systems in the Milky Way 

This method is based on counting observed NS-NS binaries in our Galaxy, and then deducing the Galactic 
merger rate by correcting for the completeness of the surveys and the estimated lifetime of each system (from 
the second supernova till the coalescence). The cosmological local rate is then derived by extrapolating the 
Galactic rate to a cosmological volume. The main advantage of this methods is that it is based on solid 
observations. As such it provides a robust lower limit on the total merger rate. The main disadvantage 
of this method is that it is possible that there are binary populations that cannot be detected by current 
surveys. For example currently detection of NS-NS binaries depends on at least one of the companions being 
a recycled pulsar. Therefore, it is insensitive to formation channels that do not recycle any of the pulsars, if 
such exist. The extrapolation from the Galactic rate to the cosmological one is also not trivial since it is not 
clear what is the right indicator which should be used for such an extrapolation (e.g., blue light?). Finally, 
the sample is small — there are only three NS-NS binaries that are used in order to estimate the Galactic 
rate - PSR B1913-hl6, PSR B1534-H2, PSR J0737-303938 . 

The first calculations of the Galactic merger rate (Phinney, 1991; Narayan, Piran & Shemi, 1991) were 
of the order of 10~^ — lO^^yr"^ , where the uncertainty was dominated by uncertain beaming correction 
and survey selection effects. More detailed calculations that followed obtained similar values (Curran & 
Lorimer, 1995; van den Heuvel & Lorimer, 1996; Arzoumanian, Cordes & Wasserman, 1999). Kalogera et 
al. (2001) carried out detailed investigation of the uncertainties in these rate estimates and found that the 
uncertainty in the correction for undetected faint systems can span over two orders of magnitude. With 
two NS-NS systems observed , PSR B1913+16 and PSR B1534+12, they estimate a rate of - lO"'^ yr'^ 
when only the beaming corrections is included. Considering the correction for faint undetected pulsars, they 
estimate that the rate can be as high as 5 x IQ-'^yr-^ Later, Kalo gera, Kim & Lorimer (2003) developed 
a method to assign statistical significance to these estimates. They find a merger rate of Slg x 10"^ yr"-"^ 
(68% confidence). The subsequent discovery of the relativistic binary pulsar PSR J0737— 3039 (Burgay et 
al. , 2003) increased the number of systems used for rate estimation to three. The corresponding range of 
Galactic merger rates becomes between 1.7 x 10^^ to 2.9 x 10^"* yr^^ at 95% confidence (Kalogera et al. , 
2004), increasing the most likely estimated rate by a factor of 6-7. This significant revision in the rate is due 
to the unique properties of PSR J0737— 3039 (short lifetime, 2.4hr period orbit and different beam profile) 
that make it difficult to detect. 

The extrapolation of the galactic rate to the local cosmological rate is usually done assuming that the 
merger rate is proportional to the blue stellar luminosity. Phinney (1991) obtains a conversion factor of 
10^ Gpc""^. Using a similar conversion factor the most updated estimate of the local universe merger rate is 
200 - 3000 Gpc-3yr-i (Kalogera et al. , 2004). 

The merger rate derived based on the observed Galactic sources is dominated by PSR J0737— 3039, which 
is short-lived (« 100 Myr). Therefore, this method predicts that the lifetime distribution of NS-NS binaries is 
dominated by rather short lived-systems (^ 1 Gyr). The fact that there are observed systems with lifetimes 
that are comparable the Hubble time or longer, implies that the observed systems lifetime distribution is 
broad. 

II. Constraints from population syntheses 

An alternative approach to evaluate the NS-NS merger rate - and the only current mean to evaluate the 
NS-BH merger rate - is via population synthesis (e.g., Tutukov & Yungelson, 1994; Brown, 1995; Lipunov 
et al. , 1995; Portegies Zwart & Spreeuw, 1996; Portegies Zwart & Yungelson, 1998; Bethe & Brown, 
1998; Fryer, Burrows & Benz, 1998; Fryer, Woosley & Hartmann, 1999; Bloom, Sigurdsson & Pols, 1999; 
Belczynski & Bulik, 1999; Brown et al. , 2000a; Belczyhski & Kalogera, 2001; Belczynski, Bulik & Kalogera, 
2002; Belczynski, Bulik & Rudak, 2002; Belczynski, Kalogera & Buhk, 2002; Perna & Belczynski, 2002; 
O'Shaughnessy, Kalogera & Belczynski, 2005; O'Shaughnessy et al. , 2005; de Freitas Pacheco et al. , 2006; 



Other observed NS-NS systems will not merge within the Hubble time and therefore are irrelevant for this method. The only 
exception, PSR B2127+11C, is excluded because of its association with a globular cluster (Anderson et al. , 1990; Prince et al. 
, 1991). 
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Dcwi, Podsiadlowski & Sena. 2006). The evolution of a binary is followed as a function of the initial binary 
properties (e.g., the zero age main sequence [ZAMS] mass of the two companions and the initial orbital 
separation). The advantage of this method is that it can be used to explore binary systems that cannot 
be observed, such as non-recycled NS-NS or NS-BH binaries. The main disadvantage of the method is that 
the uncertainties involved are substantial (e.g., the supernova kick distribution, common envelope evolution, 
etc.). As a result, the rate estimates span over two orders of magnitude, which include the range inferred 
from observed Galactic binaries. 

An interesting outcome of this research is the suggestion by Belczyhski & Kalogera (2001) that a formation 
channel of very short-lived, practically unobservable, NS-NS binaries exists. Including this channel in their 
most updated population synthesis code, Belezynski et al. (2005, 2006) find that the distribution of the time 
that the binary spends as two neutron stars, from the second supernova until coalescence, is bimodal with 
one peak at ~ 10^ yr and another at ~ 10^" yr. When the time between the formation of main sequence stars 
and the two supernovae is added, the lifetime distribution (from star-formation till death), /(r), becomes 
flat in log space (/(r) oc in the notation of §5.1) with a cutoff below 10 Myr. They obtain a rather 
similar lifetime distribution also for NS-BH binaries. 

An alternative formation channel of NS-NS binaries is by exchange interactions in globular clusters (as 
in the case of PSR B2127+11C that is associated with M15). The number of NS-NS binaries that are 
formed by this process is expected to be significantly smaller than the number of binaries that are formed 
through binary evolution (e.g., Phinney, 1991). Grindlay, Portegies Zwart & McMillan (2006) estimate a 
Galactic merger rate in globular clusters of ~ 4 • 10~^ yr~^which corresponds to a cosmological local rate of 
^ 2 Gpc~'^ yr~^. This rate is indeed much smaller than the estimates based on the observed systems that 
are not associated with globular clusters. Hopman et al. (2006) evaluate the lifetime distribution function 
of binaries that form in globular clusters and find that it is dominated by old systems with an average life 
time of ~ 6 Gyr. 

Estimates of the rate of NS-BH mergers vary by orders of magnitude, and so does the ratio between the 
NS-NS merger rate and the NS-BH merger rate. Two major open questions that determine if a NS collapses 
into BH during the binary evolution, and thus strongly affect the ratio between NS-NS and NS-BH binaries, 
is the accretion rate of a NS during common envelope phase and the maximal mass of a NS. The interesting 
work by Bethe & Brown (1998) suggests that the mergers of low-mass BH and NS should be more common 
than NS-NS mergers (see also Chevalier, 1993; Brown, 1995; Wettig & Brown, 1996; Bethe, Brown & Lee, 
2005a,b), since the first neutron star cannot survive a common envelope phase and it collapses into a BH 
following hypercritical accretion from the companion envelope. They suggest that the only way to avoid 
this outcome is by having two companions with similar masses (ZAMS masses no more than 4% apart) so 
they get into the helium burning phase at a similar time, thereby avoiding vigorous mass accretion during 
a common envelope phase. Using a Salpeter initial mass function (Salpeter, 1955) they find that the rate 
of NS-BH mergers is 20 times the rate of NS-NS mergers. This ratio is reduced to a factor of 5 when a flat 
initial mass function (constant dn/dM) is considered (Lee, Brown & Park, 2006). The mean lifetime of a 
NS-BH binary in this scenario is w 5 Gyr and the predicted rate (normalized to the supernova rate at a 
redshift of « 0.7) is ~ lO^* Gpc^^ yj.-i (Bethe, Brown & Lee, 2007, this volume). 

Finally, population synthesis can be used to put an upper limit on the local rate of NS-NS and NS-BH 
mergers, by evaluating the fraction of binaries that survive both supernovae. This fraction depends mostly 
on the distribution of the supernova kick velocity, which is not well constrained. Several works show that 
this fraction of surviving binaries is unlikely to exceed 2% (e.g., Pfahl et al. , 2002; Lipunov, Postnov & 
Prokhorov, 1997). Taking this fraction from the local rate of core-collapse SNe (Cappellaro, Evans & Turatto, 
1999) implies an upper limit of ^ 1000 Gpc~'^ yr~^ on the merger rate of short-lived binaries (t < 1 Gyr). 
An upper limit of lO'* Gpc^'^ yr^^ on the local merger rate of long-lived binaries(r > 4 Gyr) is obtained 
by relating it to the rate of core-collapse SNe at redshift 0.7 (Dahlen et al. , 2004). Therefore, if SHBs are 
mergers of compact binaries, then from Eq. 42, the SHB-merger rate must be in the range: 



10 < nsHB^rrrerger < W GpC-^ yf-l. (47) 
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5.2.3. Offsets from host galaxies and external medium densities 

An important clue about the nature of the progenitors is the location of the bursts with respect to their 
host galaxies. The offset of the merger site of a NS-NS or NS-BH binary from its birth place was explored 
by several groups (Fryer, Woosley & Hartmann, 1999; Bloom, Sigurdsson & Pols, 1999; Perna & Belczynski, 
2002; Belczynski et al. , 2006). This location depends on the natal kick of the progenitor system and the 
time between the kick and the merger, as well as the galactic gravitational potential. In the case of compact 
binaries the natal kick is expected to be between several tens to several hundreds km/s (e.g., Hobbs et al. , 
2005) and different models of binary formation channels predict a wide duration range between the second 
SN and the merger (e.g., Belczynski et al. , 2006). A binary with an inspiral time of several Gyr born with 
a kick of several hundred km/s, can merge after traveling 1 Mpc, while a binary that merges within 1 Myr 
will always end up near its birth place. Similarly, binaries with moderate kick velocity are expected to be 
bound to large galaxies but not to small ones, and therefore offsets from larger galaxies are expected to be 
smaller. Thus, in principle, mergers can take place at any distance (< 1 Mpc) from their hosts, and the 
predicted merger distribution depends on the specific binary model considered. Belczynski et al. (2006) 
find that for their model of bimodal inspiral time distribution, the offset from highly star-forming galaxies 
is very small (since the merger rate is dominated by the short inspiral time population). The offset from 
elliptical galaxies, in which only long-lived binaries merge, is expected to be large. The typical offset from a 
small [giant] elliptical is about 100 [10] kpc. 

The circum-burst gas density depends of course on the offset. If the merger site is outside of a field host, 
in the inter-galactic medium (IGM), the density is expected to be ~ 10~^ cm~'^. If, on the other hand, the 
merger takes place within its host galaxy or in the intra-cluster medium of a galaxy cluster (ICM) the density 
is expected to be > 10^'^ cm^'^ . Since the brightness of the afterglow depends on the external density (§4), 
afterglows of bursts with large offsets should be fainter. For example, Belczynski et al. (2006) present the 
distribution of external densities of merger sites predicted by their model for different types of hosts. They 
find that the lowest densities are expected for binaries that are born in small early-type galaxies, while the 
highest densities are found for those that are born in large star- forming galaxies. 

If compact binary mergers lead to SHBs then this behavior predicts a strong selection effect that biases 
the observed redshift and host-type distributions. The reason is that it is very hard to detect the afterglow of 
a burst that takes place at a large distance from its host and even if an afterglow of such a burst is detected, 
it is unlikely that a host association will be secure, and therefore the SHB redshift, can be determined. As a 
result, a bias is expected favoring bursts with short-lived progenitors, implying higher redshifts and late- type 
hosts. Similarly, a bias is expected in favor early-type cluster galaxies over early-type field galaxies. 

5.2.4. Comparison with the observations 

The most rudimentary predictions of NS-NS and NS-BH mergers fit the SHB observations well: (i) These 
mergers must take place at a rate that is comparable to the SHB rate (Eq. 42). (ii) The gravitational binding 
energy that is liberated in these mergers is more than enough to power the observed prompt emission and 
afterglow. Moreover, extensive numerical and theoretical work suggests several very plausible processes to 
channel enough of this energy into a relativistic outflow, (iii) The duration of SHBs is explained in a rather 
natural way in these models, while the engine size is compact enough to allow for the observed rapid 
variability during the prompt emission. Following recent observations, several papers have shown that the 
energy and the timescales of specific events with known redshifts can be explained within the framework 
of the merger model (Lee, Ramirez-Ruiz & Granot, 2005; Oechslin & Janka, 2006). (iv) The most popular 
engine model in this scenario, hyper-accretion onto a black hole, is similar to the engine suggested for 
long GRBs, naturally explaining the similarity between the two phenomena, (v) SHBs do not trace the star 
formation, and they take place in both early- and late-type galaxies, as expected for compact binary mergers, 
(vi) Many SHBs are observed at low redshifts (^ 0.2) while there are indications of a population of SHBs 
at higher redshift (~ 1). This result indicates on a wide lifetime distribution, as some papers predict for 
compact binary mergers (e.g., Belczynski et al. , 2006). (vii) Observations suggest that some of the bursts 
take place in very low densities and at large offsets from their hosts, as predicted for a binary population 
with strong natal kicks and long lifetimes. 
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A more detailed comparison of this model to the observations reveals that there are also some conflicts 
and issues that should still be resolved between the model and the observations. The growing amount of 
evidence suggesting that the central engine is active for a duration that is much longer than the duration of 
the initial hard prompt gamma-ray emission is a major challenge to the merger model (e.g., the late flare 
in SHB 050709). Here, the natural engine activity time scale that the merger model predicts (< 1 s) is an 
obstacle. Nevertheless, several extensions and modifications of the model that may explain such activity have 
already been suggested. An additional discrepancy between the NS-NS merger model and the observations 
is the inferred lifetime distribution of the merging systems. While the current observations point toward a 
lifetime distribution that is dominated by several Gyr old progenitors, the observed NS-NS systems in our 
Galaxy represent a population with a typical lifetime of ~ 100 Myr. Note that the poorly constrained lifetime 
distribution of NS-BH binaries is consistent with the observations of SHBs. Both the observed SHB sample 
and the observed NS-NS sample are very small, and therefore this discrepancy may be resolved with any of 
the two lifetime estimates, or both, being inaccurate. However, if supported by future SHB observations and 
assuming that Galactic binaries are representative of the cosmological population, lifetime considerations 
may rule out NS-NS coalescence and leave only NS-BH mergers as a viable SHB progenitor system. 

To conclude, given the very promising rudimentary comparison of compact binary merger model predic- 
tions and observations, this remains the leading progenitor model. However, it is important to remember 
that all supporting arguments are indirect and none of them are even close to being as conclusive as the 
detection of SN spectra in the afterglows of nearby long GRBs. Therefore, it is important, especially given 
the fact that there are some apparent discrepancies between this model and SHB observations, to explore 
the possibility that SHBs are produced by entirely different systems. 

5.3. Other progenitor models 

5.3.1. Accretion induced collapse 

An accretion induced collapse (AIC) of a rapidly rotating neutron star to a black hole was suggested by 
several authors as the source of GRBs (Vietri & Stella, 1998, 1999; MacFadyen, Ramirez-Ruiz & Zhang, 
2005; Dermer & Atoyan, 2006). The GRB central engine in this scenario is similar to the one expected in 
a compact binary merger, the rapid accretion of a disk onto the newly formed BH. Here, the source of the 
disk is the NS material with the highest angular momentum, on the equator, while the BH is formed by 
material from the collapsing NS poles. 

MacFadyen, Ramirez-Ruiz & Zhang (2005) propose this progenitor model to explain the ^ 100 s X-ray tail 
observed in some SHBs (see §2.2.2). They suggest that the collapse is initiated by accretion from a close and 
less compact companion and that the late X-ray emission is produced by the interaction of relativistic ejecta 
with the companion. Note that this scenario predicts a single smooth pulse and cannot explain variable 
X-ray tails (which might have been already observed; Norris & Bonnell, 2006) or several afterglow flares, as 
observed following SHB 050724. Dermer & Atoyan (2006) further suggest that neutrino radiation from the 
NS collapse may heat and disrupt the companion, sending some of its mass toward the young BH, leading 
to late-time flares. They also find that the rate of such events can fit the SHB rate. 

The main uncertainty in this model is the mass of the formed disk. Semi-analytic calculations and general 
relativistic simulations of the collapse into a black hole of a uniformly rotating star spinning at the mass- 
shedding limit, show that the mass of the remaining disk depends strongly on the equation of state (Cook, 
Shapiro & Teukolsky, 1994a, b; Shibata, Baumgarte & Shapiro, 2000; Shibata & Shapiro, 2002; Shibata, 
2003; Shapiro, 2004). These calculations use a polytropic equation of state and they show that assuming 
a very soft EOS, T — 4/3 <C 1, a massive disk may form (Shibata & Shapiro, 2002). However, a star with 
r > 1.5 collapses directly into a BH leaving practically no disk, Md < 10"'^ M© (Shibata, 2003; Shapiro, 
2004) . The reason is that the angular momentum of the material on the equator is not high enough in order 
to support it from collapsing. It is possible that if there are strong magnetic fields that help with angular 
momentum transfer towards the equator the outcome is different. Since neutron stars are expected to be 
much stiffer than a polytrope with F = 1.5, the current results suggest that unless strong magnetic field 
significantly alter the collapse outcome, AIC of a NS cannot be a GRB progenitor. 
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5.3.2. Magnetars 

Central engine models that involve highly-magnetized neutron stars (magnetars) with magnetic fields that 
are similar to or larger than those observed in our Galaxy 3^ (Duncan & Thompson, 1992; Paczynski, 1992; 
Kouveliotou et al. , 1998), come in two flavors. One is a newly born millisecond magnetar and the second 
is a version of SGR giant flares, in which case the magnetar already span-down significantly and is rotating 
slowly. A millisecond magnetar may be formed during the merger of two neutron stars as a transient, which 
collapses into a black hole during its spin-down, (as discussed in §5.2). A span down magnetars in old stellar 
population (as required by §5.1) must be therefore formed by a different progenitor system. A milliseconed 
magnetar that span-down without collapsing may be formed for example by a merger of two white dwarfs 
(e.g., Saio & Nomoto, 1985; Levan et al. , 2006b), or an accretion-induced collapse of a white dwarf (e.g., 
Nomoto & Kondo, 1991). 

A millisecond magnetar was first suggested as the central engine of GRBs by Usov (1992, see also Thomp- 
son, Chang & Quataert 2004; Thompson 2005). The energy source in this model is rotational energy, ~ 
5 • 10^^ erg, which is channeled into a relativistic outflow by magnetic luminosity. Assuming magnetic dipole 
emission, and that gravitational radiation can be neglected, the luminosity is (Usov, 1992): 

L^d ~ 2.5 X lO'^^Bf^nt erg/s, (48) 

where B is the magnetic field on the neutron star surface and CI is the magnetar angular velocity. The spin 
down time is: 

T^d ~ 2B^^nfs. (49) 

Therefore, SHB energy considerations require B > 10^^ G while time considerations require either B > 10^^ 
G or termination of the magnetic emission before the entire rotational energy is emitted (e.g., by collapse 
of the magnetar into a BH). 

A millisecond magnetar with a lower magnetic field (~ lO^"' — 10^^) is unlikely to produce the SHB itself, 
but it may affect the afterglow (Dai & Lu, 1998a, b; Wang & Dai, 2001; Zhang & Meszaros, 2001a; Dai, 
2004) . If somehow the end-product after the SHB is launched is such a magnetar then the rotational energy 
powers an outflow over 10^ — 10** s and this energy is injected into the external shock that produces the 
afterglow. This energy injection was suggested as an explanation of late flares and light curve flattening 
observed in SHBs 050724 and 051221 A (Gao & Fan, 2006; Fan & Xu, 2006). 

The second scenario is of an older magnetar that is producing the SHB in the same way as it produces 
SGR giant flares. Here the energy source is the magnetic field. If the field of the magnetar is similar to 
those inferred for Galactic objects 10^'^) then it cannot produce a burst with more than ^ 10"*^ erg. Dar 
(2005a) suggests that giant flares, such as the one observed from SGR 1806-20, are beamed in such a way 
that for some observers the isotropic-equivalent energy appears to be ^ 10^" erg, thereby explaining the 
whole SHB population as similar extragalactic giant flares. While this model is very economic, it fails to 
produce the observed SHB afterglows, which are brighter by many orders of magnitude than the unbeamed 
afterglow of the giant flare from SGR 1806-20 (Cameron et al. , 2005; Gaensler et al. , 2005). An alternative 
intriguing possibility is that there are old magnetars with external magnetic fields of 10^^ G (T. Thompson, 
private communication). Giant flares from such magnetars can be energetic enough to produce SHBs. Such 
magnetars may be the end product of a merger of two white dwarfs (Levan et al. , 2006b), and their birthrate 
and lifetime may be such that it is not surprising that none are observed in our Galaxy. Note that in this 
model is currently the only one where the same source may produce more than one SHB, enabling, in 
principle, repeated detection of SHBs from the same source. 

5.3.3. Quark stars 

Quark stars, if they exist, may be an entirely different energy source for both long and short GRBs. There 
are many models in which GRB progenitors involve a quark star (e.g., Schramm & Olinto, 1992; Ma & Xie, 
1996; Cheng & Dai, 1996; Fryer & Woosley, 1998; Dai & Lu, 1998b; Dar, 1999; Bombaci & Datta, 2000; 
Ouyed & Sannino, 2002; Ouyed, Dey & Dey, 2002; Berczhiani et al. , 2003; Bombaci, Parenti & Vidaha, 



see Woods & Thompson (2006) for a review of Galactic magnetars 
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2004; Bombaci, 2005; Paczynski & Hacnscl, 2005; Ouycd, Rapp & Vogt, 2005; Staff, Ouycd & Bagchi, 2006; 
Menezes et al. , 2006; Drago, Lavagno & Pagliara, 2006). In these models the burst is initiated by a nuclear 
phase transition in the composition of some, or all, of the star material, usually the conversion of hadrons 
to quarks. This phase transition is also, usually, the main source of the burst energy. The duration of the 
bursts in different models is determined by different mechanisms. For example, Menezes et al. (2006) suggest 
that the duration in which an entire neutron star is converted into a strange quark star, triggered by the 
formation of seed strange quark matter, is between 1 ms and 1 s, and may lead to a SHB. Ouyed & Sannino 
(2002) suggest that both long and short GRBs are produced by hot quark stars where the difference between 
the two types of bursts is the progenitor mass. They find that the structure of the star is different above and 
below a critical mass and that the phase transition in hot quark stars below this critical mass takes ~ 1 s 
while for a more massive star it takes ^ 80 s. Ouyed, Dey & Dey (2002) suggest a model in which the engine 
is activated by the accretion onto a quark star. In this scenario, similar to typical GRB accretion models, 
the difference between long and short GRBs is in the duration of the accretion. Dar (2005a) suggest that 
SHBs, as well as soft gamma-ray repeaters and their giant flares, originate from hyper-stars, i.e. neutron 
stars where a considerable fraction of their neutrons have converted to hyperons. The energy source is a 
continuous episodic phase transition from neutron matter to hyper matter in the inner star layers. 

The outcome of most quark star scenarios is the release of a large amount of energy in the form of 
radiation within a short time on the surface of the star. At this point these models merge with the loaded 
pair-radiation fireball model that leads to a relativistic outflow. 

5.3.4. Type la SN 

Dar & de Rujula (2003) suggested that SHBs are associated with Type la SNe. The prompt emission and 
the afterglow in this model are produced by "cannon balls" (for a review of the cannon ball model see Dar 
& de Rujula, 2004). The stringent limits on SN emission, which are a factor of ~ 1000 below that of a typical 
SN la, together with the fact that the variance on SN la peak magnitudes is small, disfavor this model. 

6. Gravitational waves from SHBs 

If SHBs are generated by NS-NS or NS-BH mergers then they are the electromagnetic counterparts of 
the most accessible gravitational-wave (GW) sources for ground-based GW observatories. As such, SHBs 
received a lot of attention from the GW community. The local SHB rate might have direct implications for 
the detection rate of GW telescopes, and more importantly searching for GW signals in association with 
SHBs may increase the sensitivity, and thus the detection rate, of these telescopes. Even if SHBs have a 
different origin they still may be powerful GW sources, although in this case the prospects for detecting 
their GW signals are not as promising. 

6.1. NS-NS or NS-BH coalescence 

Mergers of compact binaries are expected to be the first cosmological source of GW detected by ground- 
based observatories. The exact nature of the signal expected and the prospects for its detection are discussed 
in numerous papers (see Cutler & Thorne, 2002, for a comprehensive review). In principle the mergers are 
expected to have three different phases of GW emission: the chirp inspiral signal, the coalescence signal, 
and the signal from the ringdown of the remaining BH. The frequency of the typical signal during the final 
inspiral stages of NS-NS and NS-stellar mass BH systems passes conveniently through the most sensitive 
frequency range of GW observatories such as LIGO and VIRGO (~ 100 — 500) Hz, making it the main 
target for detection. At its current sensitivity LIGO-I can detect a NS-NS merger up to an average distance 
of about 15 Mpc, and a merger of a NS with a IOMq BH up to about 30 Mpc. Beginning at the end of 2005 
LIGO-I started its S5 science run, in which a full year of data is expected to be collected at this sensitivity 



The Laser Interferometer Gravitational- Wave Observatory; http://www.ligo.caltech.edu/ 
http: / /wwwcascina. virgo.infn.it/ 
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(which is comparable to LIGO design sensitivity over ahuost the entire frequency range) . The VIRGO design 
sensitivity is similar, and it is expected to be reached in the near future. 

If SHBs are generated by compact binary mergers then their local rate is given by Eq. 47, implying a 
detection rate by current observatories of 1[10] x lO""* < Tlaw < 0.1[1] yr^^ for NS-NS [NS-10MqBH] 
mergers. Therefore, in the most optimistic scenario, where SHBs are narrowly beamed and/or there are 
many dim undetected SHBs, LIGO-I might detect GW signals from undetected SHBs. In all other cases 
such a detection is unlikely. This GW detection rate is for a blind search within the data for signal from 
mergers. If the search is preformed only during some short time window (say 1 min) around the time in which 
electromagnetic emission from SHBs is observed, the sensitivity increases (a search for a significant signal 
within data with longer duration requires higher signal to noise ratio). Kochanek & Piran (1993) estimate 
that a search preformed only when SHBs are detected can increase the range of LIGO-I and VIRGO by 
a factor of ~ 1.5. If SHBs arc preferentially beamed perpendicular to the binary orbital plane the range 
is increased by an additional factor of^^ k, 1.5. Therefore, a SHB-coincident search for GW signals may 
increase the detection range to ~ 30 [60] Mpc. The detection rate of Swift is about 10 SHBs per year and that 
of HETE-2 and other IPN spacecraft is similar. If several percent of these bursts are within 50 — 100 Mpc, 
as suggested by Tanvir et al. (2005), then a coincident detection may be obtained with current facilities. 
Note that if SHBs are beamed so their GW-SHB coincident detection range is significantly increased then 
their true rate (after beaming correction) within 15 Mpc is at least comparable to the observed SHB rate 
within 30 Mpc. However, given that even the most optimistic scenarios predict a rate that is smaller than 
one merger event per year at these distances, a coincident search for SHB and a GW signal may make the 
difference between a detection and a non-detection by current GW observatories. 

Simultaneous detection of the inspiral GW signal from a compact merger and a SHB will provide conclusive 
evidence that SHBs originate from compact mergers and would improve our understanding of both merger 
physics and SHBs significantly. GWs can provide a unique view of the formation, and possibly the operation, 
of the inner engine powering the burst, which are difficult to observe via any other method. A coincident 
detection is also most valuable in order to determine the cosmological parameters (Dalai et al. , 2006), as 
the GW signal of a binary merger enables an accurate determination of the luminosity distance, while the 
electromagnetic signal may provide an accurate measurement of the redshift. If GW signal from a merger 
is detected without a coincident SHB, an association may still be secured by the detection of an orphan 
afterglow - afterglows that are not associated with prompt emission (e.g., Rhoads, 1997; Granot et al. , 2002; 
Nakar, Piran & Granot , 2002; Totani & Panaitescu , 2002; Nakar & Piran, 2003b; Levinson et al. , 2002). 
However, orphan afterglows should be searched for quickly, within days or weeks after the detection of the 
GW signal. The reason is, that the nearby SHBs, which might be detected by GW observatories, arc most 
likely at the low end of the luminosity function (< lO^^erg), simply because these are the most frequent, 
and therefore their afterglow is not expected to be detectable many days after the burst (even though they 
are nearby). Moreover, the transition to Newtonian blast wave and quasi spherical emission takes place 
days-weeks after the burst if its energy is low. 

Next generation observatories, which are planned to become operational in the first half of the next decade, 
are expected to be ten times more sensitive. For example advanced-LIGO (LIGO-II) is designed to detect 
NS-NS [NS - IOMq BH] up to about 300[650] Mpc. As discussed above, simultaneous detections will increase 
LIGO-II range by a factor of 1.5-2.5 to « 0.6[1.3] Gpc . So far, HETE-2 observed one burst at a distance 
< 700 Mpc while Swift detected at least 2 additional SHBs at a distance < 1 Gpc. The GBM detector on 
GLAST^^ is expected to have a threshold that is similar to BATSE and more than half-sky field of view, 
and thus it is expected to detect at least 5 SHBs within a distance of 500 Mpc every year. Therefore, if SHBs 
are compact binary mergers, and assuming that an efficient GRB detector will be operational simultaneously 
with LIGO-II, a coincident electromagnetic and GW detection is guaranteed. The benefit of the simultaneous 
operation of LIGO-II and an efficient GRB detector goes beyond the high likelihood to observe simultaneous 



The average range is obtained by averaging over all the possible orientations. The signal from a merger that is observed 
face-on is larger by a factor of ^5/2. 
http: / /f64. nsstc.nasa.gov/gbm/ 
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SHBs and mergers if they are associated - it will also enable to disprove the assoeiation if no simultaneous 
detections are observed. 

6.2. Other processes that radiate gravitational waves 

SHBs may emit GWs also if the progenitor is not a compact binary merger. Any progenitor that involves 
the collapse of a rotating compact object to a black hole (e.g., the collapse of a rotating neutron star 
triggered by accretion; §5.3.1) will produce gravitational waves (e.g., Stark & Piran, 1985). The amplitude 
of these waves is highly uncertain. Moreover, the absence of an accurate signal templates will reduce thier 
detectability. Such GW signals would most likely not be detected even by LIGO-II at distances much greater 
than 10 Mpc (Kokkotas & Stergioulas, 2005, , and references therein). 

If the central engine is an accretion disk - BH system then disk instabilities may break the axissymmetry 
and lead to GW emission (e.g., van Putten & Levinson, 2001, 2002, 2003; Kobayashi & Meszaros, 2003; 
Araya-Gochez, 2004; van Putten, 2005; Bromberg, Levinson & van Putten, 2006). The strength of this 
signal is uncertain and it depends on the evolution of the disk, which is different from one model to the 
other. Optimistic estimates suggest that it might be detectable by LIGO-II to a distance of ^ 100 Mpc 
(Kobayashi & Meszaros, 2003; Bromberg, Levinson & van Putten, 2006). 

An additional source of GW is the acceleration of relativistic jets (e.g., Segalis & Ori, 2001; Sago et al. , 
2004; Piran, 2005b). Unlike electromagnetic emission, GW radiation in this case is not coUimated. However, 
it is expected to be too weak for detection by LIGO-I and even with LIGO-II the detection range is not 
expected to exceed 10 Mpc by much. 
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